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FOREWORD 


This report describes the work performed by the AiResearch Manufacturing 
Company of Los Angeles, California under Tasks I and 2 of Contract NASI -8859, 
"Development of Regenerative CC^ -Removal -System Design Techniques," issued by 
NASA Langley Research Center. The purpose of Tasks I and 2 of the contract 
was to produce basic design data and design tools (in the form of computer 
programs) to facilitate the design process of, and to aid efforts in improving 
spacecraft regenerable CO^-removal systems which employ molecular-sieve and 

silica-gel inorganic sorbents. The time period of the analysis and experimen- 
tation that this report represents was from November 1968 to May 1972. 

Two CO^-removal -system performance-prediction computer programs were pro- 
duced under this contract. These programs are discussed only briefly in this 
report; full documentation of these programs is contained in a separate report 
"A Transient Performance Prediction Method for C0 ? . Removal with Regenerable 
Adsorbents", NASA CR-I 12098 (Reference l). * 

NASA Manned Spacecraft Center provided a portion of the financial support 
for the dynamic coadsorption mass-transfer tests; these tests were monitored at 
NASA MSC by Mr. Will Ellis. Dr. Sam Davis served as a consultant to NASA MSC 
in this area of interest. 

At AiResearch, Mr. J. P. Byrne was the initial program manager; later, 

Mr. C. W. Browning was program manager. Initially, Mr. J. M. Ruder was 
principal investigator in charge of all equilibrium data taken in Tasks I 
and 2. Later, Dr. R. M. Wright was principal investigator responsible for 
dynamic mass-transfer experimentation. Dr. K. C. Hwang was responsible for 
the performance-prediction computer programs developed under the contract. 

Dr. Jack Winnick of the University of Missouri was a consultant to AiResearch 
during the early phases of the work. 

Experimental equipment design was handled by Mr. K. Walther, Mr.' M. Gee, 
and Mr. W. Dunn. Laboratory technicians were Mr. R. Forrey, Mr. D. Falzett, 
and Mr. R. Cl i fton. 
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SUMMARY 

This report describes the results of analytical and experimental research 
studies which were directed toward the acquisition of basic design information 
on molecular sieve sorbents for C0 2 removal systems. Equilibrium and dynamic 
mass-transfer data are presented for a number of molecular sieve sorbents 
over a range of conditions appropriate to future manned spacecraft C0 2 - 
removal systems. 

A large amount of equilibrium data for molecular sieve and silica gel 
sorbents was obtained and is presented here in graphical form. Some of the 
equilibrium data is for coadsorbed pairs of sorbents, with water and CC^ 
coadsorption of primary interest. 

The dynamic mass-transfer behavior of special prototype beds was studied 
under control led condi tions. This involved studies with single sorbates in 
the carrier gas, and tests where the coadsorpt i on of CO^ and water was studied. 

Many observations resulted from the latter series of tests. From the single- 
sorbate tests, mass-transfer coefficients and intraparticle d i ffus i vi ti es 
were evaluated and are presented herein. Comparisons were made of different 
size plate-fin heat-transfer surfaces placed within the packed sorbent beds. 

Packed-bed densities and process-gas flow pressure drops were measured 
for various bed heat-transfer-surface configurations. 

Two large-scale computer programs were produced to predict the transient 
performance of spacecraft, regenerable, COg- removal systems which employ molec- 
ular sieve sorbents. One program considers water-dump/CO^-dump systems which 
use composite sorbent beds. The second, more-general program allows for two- 
or four-bed systems, with water-save and CO^-save provisions. It considers 

the second-order effects due to the coadsorption of CO^, water, oxygen, and 

nitrogen. Several auxiliary programs useful in molecular-sieve CO^-removal- 
system design work were produced. 

Along with presentation of graphical and numerical data, pertinent obser- 
vations and conclusions have been made. Suggestions for further design and 
research activity are made. 
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SYMBOLS 


Specific external sorbent area exposed to the process gas stream, 
ft 2 /ft 3 of bed. 

2 

Superficial bed cross-section area for flow, ft 
Concentration, mass of C0 2 or HgO/mass of N 2 

o 

Intraparticle diffusivity for sorbate k, ft /hr 

Bed void fraction 
Flow length, ft 

Differential heat of adsorption Btu/lb of sorbate 
Constant in pressure drop equation 

2 

Gas-phase mass-transfer coefficient, lb-mole/hr-ft -mm Hg 

Mass of sorbent contained in the test bed, lb 

Molecular weight of sorbate k; molecular weight of process gas 
stream 

Partial pressure of sorbate k, mm Hg 

Total pressure, mm Hg 
Pressure drop, psi or in. H 2 0 

Radial distance from sorbent pellet center, ft 
Universal gas constant, 1.986 Btu/lb-mole-°R 
Time, adsorption-breakthrough time, hr 
Temperature, °F or °R 
Interstitial gas velocity, ft/hr 
Superci f icial gas velocity, ft/sec 



w 


Sorbate loading on the sorbent, mass sorbate/mass sorbent; or 
as a percent loading, mass sorbate/100 mass units of sorbent 

W Process gas flow rate, lb/hr 

x Axial distance from front of bed, ft 

0 Perfect-bed breakthrough time, hr. Defined by Equations 4-3 

and 4-4. 

P Density, lb/ft^ 

a Ratio of actual density in lb/ft 3 to that at standard conditions 

cr = p/0.07176 

t Dimensionless time, r = t/9 

i(f Dimensionless concentration, defined by Equation 4-6 

Subscri pts 

g Process gas 

k Sorbate k; for example COg or HgO 

s Sorbent 

sb Sorbent bed 

sg sorbent-to-gas 

Superscr ? pt 

* Refers to equilibrium condition 
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SECTION I 
INTRODUCTION 


In all Mercury, Gemini, and Apollo manned spacef l ights, carbon dioxide 
has been removed from the spacecraft cabin by means of absorption by 
lithium hydroxide. This method of C0 2 removal is extremely reliable and 
efficient, and in general, the design of LiOH systems is relatively easy. 

The main objection to LiOH is that it is not regenerable and as an expend- 
able it can represent a large weight penalty to the vehicle, especially 
for long duration missions. 

Considering the normal CO 2 output of a man to be 2.2 lb/day, and using a 
92-percent utilization efficiency without contingencies, 2.60 lb of LiOH are 
required per man-day. For first generation manned-spacecraft, the penalties 
associated with LiOH were reasonably small. A one-day Mercury flight required 
about 5.4 lb LiOH, contained in a single bed. For the two-man Gemini vehicle, 
one large LiOH bed was employed which could hold 80.5 lb LiOH for a 14-day 
mission. The Apollo vehicle employs a dual-bed LiOH C02-removal unit; each 
bed contains a removable LiOH charge of 3.979 lb. Each bed Is charged by the 
crew every 24 hours; bed recharging is alternated to provide one charge every 
12 hours. For a 12-day mission, some 24 LiOH charges (95.5 lb) are carried. 

For long duration missions with large crews, the total weight and volume 
penalty due to LiOH could be substantial. Also, the lower cabin levels of 
CO 2 and higher total pressures desired in future missions makes the design of 
LiOH systems more complex. This is due to the accumulation in the bed of water 
produced by the chemical reaction with LiOH, and the subsequent slow-down in 
the CO 2 removal process. The water accumulation is a result of excessive 
cooling of the bed due to the large flow of cabin gases. In previous missions, 
at 7 mm Hg CO 2 partial pressure and 5 psia total pressure, the mass ratio of 
CO 2 to cabin gases is reasonably large, and the heat of reaction of COp and 
LiOH is sufficient to keep the bed warm and to drive off water formed By the 
reaction. For future mission conditions, the mass ratio of CO 2 to cabin gases 
is much lower and the bed would run much cooler. To make the system work, 
thermal conservation by means of regenerative heat exchangers is necessary. 

The system and its design are thereby complicated. 

As a consequence of the penalties and added complexities of LiOH 
systems for long-term flights, there has been a great deal of emphasis 
to develop low-penalty, regenerable (^-removal systems. One of the most 
promising concepts involves the use of inorganic molecular-sieve sorbents. The 
concept has been used by Ai Research in a system for Skylab missions. This 
system is capable of maintaining the partial pressure of CO 2 in Skylab at 5 mm 
Hg or below for the three-man crew. The entire system weighs 195 lb; two entire 
systems are included in the vehicle for redundancy, although only one is re- 
quired for the CO 2 control. The normal energy requirement (exclusive of fans) 
is 87 watt-hr per day. The flow through the system is 15 lb per hr; cycle time 
is 30 minutes. Each system contains two beds. Each bed contains 
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7.5 lb of molecular sieve 1 3X and 9.9 lb of molecular, sieve 5A. In addition, 

10 lb of activated charcoal is included in* each RCRS package. The Skylab 
vehicle is to be manned for a total of I 12 days. If LiOH were used for COo 
removal, 874 lb of the material would be required; hot to mention the weight 
of the attendant hardware. 

At the time that this program was initiated (during the early development 
period of the RCRS), it was obvious that there were many areas in the design of 
molecular-sieve CO 2 removal systems where data quantity and quality were mar- 
ginal or insufficient, and that many areas existed where improvements could be 
made in design techniques. Thus, this program was directed at providing infor- 
mation and design tools where there were known deficiencies. .This report pre- 
sents the results of these efforts. Some of the results of the program were 
immediately used in the RCRS development. 


With any learning and development process, as more information becomes 
available and as desired applications change, new areas emerge which warrant 
study. This has occurred with respect to molecular-sieve CC^-removal systems. 
Throughout this program, where it has been possible, new areas of concern 
were incorporated into the study investigations; however, it was not possible 
to cover all such areas. In some cases, the need for further study has become 
apparent only recently; some of these areas are discussed in the last section of 
this report. In summary, it is considered that the concept of CO 2 removal by 
molecular-sieve sorbents is sound and practical for manned spacecraft. Further, 
the concept has considerable growth potential and thus should be considered 
as a leading candidate for CO 2 removal on future spacecraft. 

Program results are presented in the remainder of this report. The 
contents of the report are summarized below: 


e Section 1 Introducti on 

© Section 2 Basic Scope of the Study Program — Background informa- 

tion on spacecraft molecular-sieve CC^-removal systems 
and a review of the scope of the program. 


© Section 3 Equi 1 ib.rium Data --Equi 1 ibrium data for COg, .water, O 2 , 

and N 2 on various molecular sieve types, including 
coadsorption of sorbate pairs on molecular sieves. 


o Section 4 Single Sorbate Dynamic Mass-Transfer Data- -Si ngl e- 

sorbate dynamic breakthrough tests and comparisons of 
packed-bed heat-transfer surfaces. 

© Section 5 Evaluation of Mass-Transfer-Coefficients and Intra- 

particle Pi ffusi vi ties — •This evaluation based on the 
tests given in Section 4. 


© Section 6 Coadsorption Dynamic Mass-Transfer Tests- -Coadsorption 

dynamic breakthrough tests, involving CO 2 and water 
simultaneously adsorbed by molecular sieves. 
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o 


Section 7 

Section 8 
Section 9 


Performance-Prediction Computer Prog rams --Computer 
programs developed under this contract. 

Pressure Drop and Packed-Bed Densities — Data for 
molecular sieves in internally finned beds. 

Conclusions and Recommendations 
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SECTION 2 

BASIC SCOPE OF THE STUDY PROGRAM 


GENERAL 

This section presents background information on molecular sieve adsorbents, 
and on features of spacecraft C02-removal systems using molecular sieves. On 
the basis of the background information, this section reviews the scope of the 
research and analyses performed during this program. 


REGENERABLE MOLECULAR-SIEVE COg-REMOVAL SYSTEMS FOR SPACECRAFT 


Basic Characteristics of Molecular Sieves 

Molecular sieves belong to the class of materials called zeolites which 
are generally crystalline, hydrated metal alumino Silicates. Many types of 
zeolites are known. The molecular-sieve type of zeolite has a special crystal- 
line structure which contains many interconnecting cavaties, all of uniform size. 
The cavities, or "cages," are interconnected by narrow openings, sometimes 
called windows. When formed, the cages are filled with water. Upon heating, 
the water can be driven out leaving a tremendous internal surface area capable 
of adsorbing many different gases and liquids, especially polar molecules like 
CO 2 and water. 

When gases are adsorbed by molecular sieves, heat is given off in the 
process. Simply speaking, the heat release is approximately equal to the 
latent heat of vaporization of the gas at the adsorption temperature, plus a 
heat of adsorption. For water, the total heat release is about 1400 Btu per 
lb of adsorbed water. It is from this rather large heat effect associated with 
water that the name zeolite is derived: from the Greek, "zeo" meaning boil, 
and " 1 i thos" meaning stone. Early observers noted that the natural zeolite 
materials seemed to melt and boil when heated--as water was given off. Also, 
in the baked-dry condition, rather violent heat releases can occur, simulating 
boiling, when zeolites are dumped into water. Some zeolites literally explode 
when plunged into water. 

The unique property of molecular sieves--the property for which they are 
named--is that the openings between cavities are very uniform in size and 
small enough (3 to 10 angstroms, depending on the crystal type) to prohibit 
the entry and subsequent adsorption of large molecules. Thus, in many 
instances separations can be made strictly on the basis of the size of the 
molecules. 

The chemical process industry today uses millions of pounds of molecular 
sieves in many diverse applications. The great majority of molecular sieves 
are made synthetically. There are three manufacturers of molecular-sieve or 
zeolite sorbents in the United States: Linde Division of Union Carbide Corp., 
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Davison Chemical Division of W.R. Grace and Co., and Norton Company. Of major 
interest are the synthetic molecular sieves of the type A crystal structure 
and the type X crystal structure. With each of these types, the building 
blocks of the crystal are alumina and silica tetrahedra. These tetrahedra are 
each formed by four oxygen atoms surrounding an aluminum or silicon atom. 
Considering that neighboring tetrahedra share oxygen atoms (i.e., each oxygen 
atom is shared among four tetrahedra), silica tetrahedra, with silicon in the 
+4 valence state, are electrically neutral. However, with aluminum in the +3 
valence state, the alumina tetrahedra are negatively charged. This requires 
the presence of a cation in the crystal structure. The cations are exchange- 
able and are located in the vicinity of the windows. By altering the size 
and valence of the cations, the window size may be made larger or smaller. 

For example, with sodium cations, molecular sieve type 4A has 4 . 2-angstrom 
windows. This is the parent A-type molecular sieve; it has the basic formula: 

Na |2 [(AI0 2 ) |2 (Si0 2 ) |2 j x H 2 0 

1 , 

By exchanging calcium ions for sodium ions (in the ratio of two sodiums for 
each calcium), type 5A molecular sieve is produced with approximately 
5-angstrom wi ndows. With potassium ions exchanged for the sma 1 1 er sodium 
ions, the openings are closed to approximately 3 angstroms; the molecular 
sieve is called type 3A. Similarly, with the X-type crystal, types I OX and 
I 3X are available. Type I 3X has the chemical formula: 

Na 86 [ CAI0 2 ) a6 (:Si0 2 > IO^ * H 2° 

It has windows of about 9 angstroms. 

Because of their ' inorganic crystalline nature, molecular sieves are 
extremely inert and stable. For example, with normal commercial regeneration 
temperatures of about 600° to 700°F, an indefinite number of regenerations of 
the basic crystal are possible without affecting adsorption performance. Excel- 
lent repeatability of performance has been experienced with molecular sieves. 

Most molecular sieves are first manufactured in the form of fine 
crystallites. For 4A and 5A molecular sieves manufactured by Linde, the 
crystallites are cubic with crystal sides measuring from I to 5 p,, with the 
average being about 3 p, (Reference 2). This material (usually referred to 
as a powder) is not practical for packed-bed operations (pressure drop in 
beds would be too large); therefore, the base powder is usual ly combined with 
an inert clay binder and formed into pellets or spheres of practical dimen- 
sions (e.g., 1/16 to 1/8 in.). 

( 

The final molecular-sieve pellet or sphere is about 20 percent clay 
binder. It is desirable that the binder be as inert as possible, both 
chemically and with respect to adsorption; it should be strong, to allow long 
sorbent life; and it should possess an open-pore structure to allow unhindered 
diffusional access to the molecular-sieve crystallites. Each manufacturer 
claims to have produced a molecular-sieve sorbent which is near ideal with 
respect to these properties. In addition, bi nderl ess sorbent particles have 
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been produced. Early binderless parti cles were reputed to be fragile and to 
have a rather restrictive pore structure. Newer binderless sorbents are 
advertised as being much stronger and equal in dynamic performance to the 
standard sorbent particles with clay binder. 


Characteristics of Spacecraft COg-Removal Systems using Molecular Sieves 

Require ment for Water Predrier Beds . - Molecular sieve 5A has an equilib- 
rium loading capacity ‘for CO^ of 0.068 gm CO^/gm 5A at 70°F and a C0„ partial 
pressure of 7 mm Hg. At the same temperature and at a water partial pressure 
of 9.2 mm Hg (50 U F dew point), molecular sieve 5A has an equi librium capacity 
for water of about 0.18 gm H 2 <Vgm 5A. However, the material will not simul- 
taneously adsorb both CO 2 and water to these capacities/ In fact, when both 
are present, molecular sieves unfortunately have a strong preference for 
water and will absorb much less C02* Further, if adequate water vapor is 
available, CO 2 will be stripped off molecular sieves, resulting in no net COg 
adsorption. 

This hydrophilic nature of molecular sieve sorbents has an important impact 
on spacecraft C 02 -removal systems: each system must contain a predrier bed 

upstream of the molecular-sieve CC^-removal bed. This predrier (or desiccant) 
bed must be capable of trapping almost all of the water in the process air 
stream, allowing very little to reach the CC^-removal bed. A removal efficiency 
of 99.5 percent is considered as a design goal for the predrier. When a 
molecular-sieve bed has adsorbed water to the point that its C02-removal per- 
formance is impaired, it is said to have been "poisoned." This poisoning is 
reversible in nature — like almost all adsorptions with molecular sieves. By 
heating the bed sufficiently to drive off the water, the C02-removal perfor- 
mance is totally restored to the bed. 

two sorbents have been proposed for the predrier bed: molecular sieve 

I3X and silica gel. Because of this use of silica gel, it has been included 
as a subject in this study. Molecular sieve I 3X has less overall capacity for 
water at spacecraft cabjn conditions than silica gel, but it has substantially 
greater capacity in the low partial -pressure (low concentration) range. The 
choice of sorbent to be used for a particular application depends upon design- 
analysis tradeoffs. 1 For the Skylab RCRS, molecular sieve I3X is used. 

Cyclic Adsorption/Desorption of Beds . - Molecular-sieve CCL-removal sys- 
tems operate in a cyclic fashion with two or, in rare cases, three sets of 
adsorption beds. While one set is adsorbing water and €0^, the other set is 
undergoing a desorption process. When the respective adsorption processes 
have reached the desired limits, the process gas is switched to the freshly 
desorbed bed system, and the water- and COg-laden beds undergo the desorption 
process. 

For adsorption, cabin gas is passed through the two beds (predrier and 
CO -removal) and returned to the cabin. Low head fans are sufficient for this 
application. As the adsorption process occurs, heat is given off in both beds. 
The design of the system must fully account for this heat release, as the 
molecular sieve adsorption capacities for CO^ and water are inversely related 
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to temperature. As the beds warm, their equilibrium capacities for CO2 and 
water become less. 

There are two approaches to providing thermal control of adsorbing beds. 
First, in an active manner, cooling can be supplied to the beds from the space 
craft heat-transport system. This requires that cooling coils or passages be 
incorporated within the packed beds. The more efficient the cooling process 
(with the ultimate being an isothermal bed) and the lower the bed temperature, 
the more efficient will be the adsorption and the smaller will be the allowabl 
bed size. The major disadvantage to internal bed cooling is that an interface 
is created with, and a dependency placed upon, the spacecraft coolant system. 
This latter system also sustains a penalty due to thi s appl i cation. 

The second approach i s to forego all active cooling and allow the adsorp- 
tion beds to operate adiabatical ly. Basically, a larger bed mass is provided 
so that the heat released in adsorption does not raise the bed temperature to 
an intolerable level. The bed temperature floats somewhat above the inlet 
temperature of the process gas stream (usually cabin temperature). With this 
approach, a great deal of simplicity is allowed in the system hardware design; 
for example, internal coolant passages are dispensed with. However, in order 
to assure proper performance, very accurate and encompassing analyses are 
required. The RCRS for Skylab employs adiabatic beds, thus making the system 
independent of the spacecraft heat-transport system. 

In general terms, desorption is accomplished by subjecting the bed to 
conditions where the equilibrium part ia 1 -pressure of the sorbate on the bed 
is greater than the actual partial pressure above the bed. This can be 
accomplished by subjecting the bed to an envi rontnent where the sorbate is in 
low concentration; it can also be accomplished by increasing the temperature 
of the bed, thus increasing the equilibrium partial pressure of the bed; 
or both methods may be used. If desorption is accomplished by subjecting the 
bed to low pressure (and thus a low partial pressure of the sorbate), the 
bed is said to be operating in a "pressure swing" mode. That is, the bed pres 
sure swings from the adsorption level, usually cabin pressure, to the low leve 
employed for desorption. If the low desorption pressure is provided by space 
vacuum, the bed is said to be "vacuum dumped." If desorption is carried out 
by heating the bed, the operation is called "thermal swing." If desorption 
is carried out by flowing through the bed a gas which has a very small concen- 
tration of sorbate, the process is called "gas stripping.;" Combinations of 
the desorption methods are possible; for example, a bed might undergo a com- 
bined pressure/thermal swing if the bed is heated as it is exposed to vacuum. 
Heating (thermal swing) can also be used with gas stripping. 

Types of Spacecraft COp-Removai Systems . - The desorption method to be 
applied to a particular application depends greatly upon the ultimate desired 
dispos i tion of the separated sorbates, CO2 and water. If there is no need to 
save the adsorbed CO2 and water, then simple vacuum dumping of both predrier 
and CC>2-removal beds is probably the best available desorption method. In 
such a water-dump/C 02 “dump system, it is possible to pack both sorbents in a 
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Figure 2-1. Composite-Bed, Water-Dump/C02-Dump CO 2 Removal 

System; After the Regenerative C02~Removal System 
(RCRS) for the Skylab Vehicle 

single bed container. Figure 2-1 depicts this type of a composite-bed, water- 
dump/C02 - dump C : 02-removal system, which is often called a two-bed system. 

The Skylab RCRS is precisely this type of system. 

With composite bed systems, it is almost mandatory that the desorption 
take place in such a manner that the water leaving the predrier portion of the 
bed does not travel through, and have the opportunity to adsorb on, the 002 - 
removal bed. In effect, the desorption is counterflow to the adsorption. As 
a somewha.t more complicated alternative, a void space may be provided between 
predrier and C 02 -removal portions of the bed. Upon desorption, this void 
area can be opened to vacuum. By also opening the predrier inlet to vacuum, 
double-end desorption of the predrier bed is afforded. This arrangement is 
shown in Figure 2-2. 

If the adsorbed water, but not the adsorbed CO 2 , must be saved, another 
system arrangement is required. It may be necessary to save water either to 
maintain cabin humidity conditions or to reclaim and use all water aboard the 
spacecraft. Such a system would be known as water-save/C02-dump system. 

For this system, it is necessary to physically isolate the predrier beds from 
their respective C02-removal beds and to provide a different means of desorbing 
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Figure 2-2. Two-Point Desorption Scheme for Regeneration of Composite 
Sorbent Beds 


the predrier beds. This can be accomplished by gas stripping. By proper 
ducting and valving, it is possible to introduce the very dry, C 02 -lean gas 
exiting from the adsorbing CC^-remova 1 bed back into the desorbing predrier 
bed, as shown schematically in Figure 2-3. That is, the desorpti on of water 
from the predrier beds is carried out by a gas-stripping arrangement, using 
all or a portion of the same process gas stream that just passed through the 
adsorbing beds. Thus, as can be seen from Figure 2-3, at any particular time 
the process gas stream flows through three beds: the predrier bed, where water 

is removed from the processed gas stream; the adsorbing CO 2 bed, where CO 2 is 
removed; and finally back through the previously water-loaded predrier bed. 

In this arrangement, care must be taken that the gas-stripping operation is as 
efficient as the adsorption operation. Meanwhile, the desorbing CC^-removal 
bed is opened to vacuum for a more or less simple vacuum dump. Upon comple- 
tion of the cycle, the bed positions are essentially reversed. 
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Figure 2-3. Four-Bed Water-Save/C02-Dump CO 2 Removal System. 


With somewhat greater system complexity, it is possible to produce a semi - 
water-save/C02“dump system. This might be desired where there is an overbun- 
dance of water in the cabin, and it is desired to dump water to vacuum for a 
fraction of the desorption time. In such a system, a gas-stripping operation 
of the predrier bed might proceed for a fraction of the desorption period, and 
then the bed would be opened to vacuum to finalize the desorption and to remove 
excess water from the spacecraft cabin. This, however, requires a rather 
complex system arrangement, which is probably undesirable. In most instances, 
the same result can be obtained by using a composite two-bed system and altering 
k the process gas flow rate to fit the production rates of CO 2 and water. This 
type of system is based on the assumption that the production of CO 2 and water 
are linked, and that by controlling cabin humidity by the modulation of the 
process gas flow rate, CO 2 removal is also properly controlled. 
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For extremely long duration flights, it appears desirable to save both 
water and CC^- A water-save/C02-save system differs from the water-save/C02- 
dump system only in the manner of CO 2 desorption. Water is again removed from 
a predrier bed by gas stripping; however, the CO 2 is removed from the CO 2 - 
removal bed by means of heating and vacuum pumping, the CO 2 being delivered to 
an accumulator which would be upstream of an oxygen reclamation unit. 


SCOPE OF RESEARCH AND ANALYSES 

The research and analyses conducted under this contract for the purpose of 
generating basic design data and providing design tools for spacecraft regener- 
ate C0 2 -removal systems was in three areas: (l) experimentally defining basic 

equilibrium data for molecular-sieve sorbents, (2) evaluating the dynamic mass- 
transfer performance of prototype beds under representative conditions, and 
(3) producing general performance-prediction computer programs for spacecraft 
C02-removal systems. 


Equi 1 i bri urn Testi ng 

Equilibrium sorption capacity refers to the ultimate capacity of a sorbent 
to adsorb another material (sorbate) under given conditions of temperature and 
sorbate partial pressure. Equilibrium sorption capacity data is usually taken 
by exposing a sorbent sample to a controlled environment of temperature and 
sorbate pressure. After a sufficiently long period of time, when equilibrium 
between gaseous and adsorbed sorbate has been established, the weight increase 
of the sorbent material is noted. The sorbent capacity or sorbent loading is 
usually expressed in terms of the mass of sorbate adsorbed on a unit mass of 
sorbent; percentage loadings (mass of sorbate per 100 mass units of sorbent) 
are sometimes used. Because sorbent loading is a function of temperature and 
sorbate partial pressure, equilibrium data is usually displayed as a map of 
loading vs partial pressure, with temperature used as a parameter. The lines 
drawn for constant temperature are called adsorption isotherms. 

Equilibrium data may be applied in a number of ways. First, by its 
definition (expressing maximum sorption capacities), it is useful in initial 
comparisons of sorbents or sorp'tion conditions. Because equilibrium data 
says nothing about rates, such comparisons can be used only as guidelines. 
Properly used with rate information, the comparisons become more definitive. 

Equilibrium data is also used in the prediction of the dynamic perfor- 
mance of beds, as the mass-transfer driving force is usually considered to be 
the difference between the actual sorbate partial pressure in the gas stream 
(p^) and the equilibrium partial pressure ( p^*, evaluated from the equilibrium 
map) with respect to the instantaneous bed temperature and loading. Thus, 
equilibrium data is necessary for the use of performance-prediction computer 
programs. Another use of equilibrium data is that the shapes of the adsorption 
isotherms can give insight into the relative ease or difficulty of adsorption 
processes. 
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In spite of the basic nature of equilibrium data, at the onset of this 
program much of the required equilibrium data was scarce, was missing entirely, 
or was questionable. A substantial part of the reasoning why such data was not 
available lies in the range of partial pressures of interest for spacecraft 
systems as contrasted to those of interest in commercial practice. For example, 
CO 2 removal in industry is not usually accompl i shed wi th the use of solid 
adsorbents. Even when CO 2 adsorption by molecular sieves occurs commeri cal ly, 
CO 2 partial pressures are usually much higher than the levels experienced in 
spacecraft. 

With regard to water equilibirum on molecular sieves, it should be noted 
that even at very low partial pressures these sorbents can hold large quanti- 
ties of water. However, it is only recently that high-quality instrumentation 
has become available (l) to accurately measure water vapor pressures in the 
range of interest (down to below Ip,), and (2) to measure, within the apparatus 
used for equilibria determination, small weight changes of sorbent samples. 

Thus, much of the previously available data is questionable. To clarify the 
interest in low water-vapor pressure, it is noted that Linde 5A has an equili- 
brium capacity of about 8.5 percent at 70°F and at a water vapor pressure of 
only 0.006 mm Hg; that is, at -80°F frost point. Further, at 8.5 percent water 
loading, molecular sieve 5A exhibits a drastic poisoning effect with respect 
to CO 2 adsorption. 

For many new sorbents (for example, sorbents with improved binders) there 
was little equilibrium data by which to judge their potential usefulness. 
Furthermore, there was essentially no data at all on the equilibrium capacities 
of molecular sieves when both water and CO 2 were present. As discussed earlier, 
because of the poisoning nature of water, it is most desirable to have quanti- 
tative information on the equilibrium capacity for CO 2 in the presence of water 
vapor; therefore, a good deal of the equilibrium data taken in this contract 
has been for CO 2 and water simultaneously on molecular sieves. There was also 
interest in sorbate pairs other than CO 2 and water. Thus, data was also taken 
for CO 2 paired with oxygen, CO 2 paired with nitrogen, water paired with oxygen 
and water paired with nitrogen. The data and explanations of the equilibrium 
tests that were conducted are presented in Section 3. 


Si ngl e-Sorbate Dynamic Testing 

Equilibrium data is probably the most basic of al 1 categories of mass- 
transfer information. .However, the knowledge of the equilibrium sorption 
capacity of the particular sorbent is not sufficient to declare that it is a 
good sorbent, and that it will perform well under dynamic conditions. One 
of the goals of this program was to characterize molecular-sieve sorbents 
for water and CO 2 under representative dynamic conditions. To accomplish this, 
special beds were produced which incorporated the latest concepts in internal- 
coolant heat-transfer surfaces. Using these beds, si ngl e-sorbate mass-transfer 
experiments were carried out under isothermal conditions. The reason for using 
isothermal adsorptions was to eliminate heat-transfer considerations from 
evaluation of mass-transfer parameters. Therefore, to the greatest extent 
possible, the beds were kept isothermal during the entire adsorption runs. 
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The sorbent/sorbate combinations studied in the dynamic tests were: 
water adsorbing on silica gel; water adsorbing on Linde I3X molecular sieve; 

CO 2 adsorbing on Linde 5A molecular sieve; CO 2 adsorbing on Davison 5A molecular 
sieve; and CO 2 adsorbing on Linde I3X molecular sieve. All of the single- 
sorbate adsorption tests were initiated by subjecting the bed to a prolonged 
vacuum-bakeout to remove essentially all traces of CO 2 and water. All adsorp- 
tions were carried out with nitrogen as the carrier gas for the CO 2 or water 
vapor. The nitrogen/C02 or nitrogen/water gas stream was passed through the 
sorbate bed until the outlet concentration from the bed finally rose to that 
of the inlet concentration; that is, until complete breakthrough of the 
sorbent was obtained. For some cases, gas stripping desorptions using dry 
nitrogen were carried out after adsorption tests. The data from the break- 
through runs and gas strips, and the discussion of the results of these tests, 
are contained in Section 4. 


Evaluation of Mass-Transfer Coefficients 

Using the single-sorbate breakthrough curves as basic data, and with the 
help of the performance-prediction computer programs produced under this program, 
mass-transfer parameters were evaluated. For the sorbent systems in question, 
the mass-transfer parameters of interest are the mass-transfer coefficient Kg 
and the effective intraparticle diffusivity, D- The mass-transfer coefficient 
relates to the conductance of the convective mass-transfer process as the 
sorbate moves from the free stream of gas within the bed to the sorbent-pellet 
boundary. The effective sorbate diffusivity relates to the diffusional con- 
ductance as the sorbate travels through the pore structure of the sorbent 
pellet. The means of determining the mass-transfer coefficients and diffusi- 
vities from the breakthrough data, and the results of these evaluations are 
contained in Section 5. 


Coadsorption Dynamic Testing 

In order to obtain a greater understanding of the water poisoning effect 
on C02-removal beds, a series of coadsorption tests were conducted in the 
dynamic mass-transfer apparatus. A few of these runs involved the simultaneous 
adsorption of CO 2 and water onto a sorbent bed. However, most of the coadsorp- 
tion runs were conducted in somewhat different manner: After an initial 

vacuum-bakeout, the bed was allowed to acquire a preload of water, up to a 
designated loading, by means of adsorption from a water-bearing nitrogen 
stream. The water preload was usually 7 to 9 percent, corresponding to an 
equilibrium water vapor pressure of -80°F frost point. Following the water pre- 
loading, a stream of nitrogen and CO 2 was passed through the bed and its C 02 “ 
removal performance monitored. Section 6 contains the data and the discussions 
of these rather unique tests. 
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Performance-Prediction Computer Programs 

Section 7 contains discussions of the computer programs produced under 
this contract. Two major C02-remova 1 -system performance-prediction computer 
programs were produced. The first program concerns composite-bed water-dump/ 
C02-dump systems. The second and more general computer program deals with 
four-bed systems that involve any combination of water-save, water-dump, CO 2 - 
save, or C02~dump provisions. Further, it has features to allow the prediction 
of the gradual buildup of water on the CC^-removal bed over long periods of 
operation, and the prediction of the CO 2 -removal performance of the bed in the 
face of the poisoning effect of water. Additional programs useful in CC^- 
removal- system design and evaluation efforts are also discussed and presented 
in Section 7. 

Pressure Drop and Packed-Bed Density 

Section 8 contains design information on packed-bed density and pressure 
drop for molecular-sieve and silica-gel beds. 
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SECTION 3 
EQUILIBRIUM DATA 


INTRODUCTION 

This section presents the results of tests performed oh various molecular- 
sieve materials to determine their equilibrium sorption capacity for specific 
gases at various pressures and temperatures. The equilibrium sorption capacity 
refers to the ultimate capacity of a sorbent to adsorb, another material under 
given, constant conditions of temperature and sorbate partial pressure. The 
sorbate capacity, or loading, is usually expressed in terms of the mass of 
sorbate adsorbed on a unit mass of sorbent, or alternatively, as a percentage 
loading, i.e., mass of sorbate per 100 mass units of sorbent. 

This section is divided into four subsections: 

• Equilibrium test equipment and procedures - Describes the test 
setup and the test equipment used; defines test techniques used 
to obtain specific data; and discusses some of the problem areas 
encountered. 

o S i ngl e-sorbate equilibrium test data - Includes all data taken 

pertaining to single gas or vapor adsorption by various molecular 
sieve and silica gel sorbents. Where applicable, this data is 
compared with previously existing data. 

• Coadsorption equilibrium test data - Presents all data taken where 
a sorbent was preloaded with gas or vapor to some predetermined 
level and then subjected to the partial pressure of another gas 
sorbate. The . equi librium coadsorption behavior of several sorbate 
pairs on Linde 5A and on Linde I3X have been studied and are 
presented in this section. 

• Heat of adsorption - Discusses related heat of adsorption associated 
with the si ngl e-sorbate adsorptions. 


EQUILIBRIUM TEST EQUIPMENT AND PROCEDURES 


Test Equipment 

The equilibrium test setup schematic is shown in Figure 3-1; a photograph 
of the apparatus is shown in Figure 3-2. The basic purposes of the apparatus 
are (l) to provide the proper, controlled adsorption environment for the sor- 
bent test sample, and (2) to measure this environment and, by weight, to 
obtain the adsorption capacity of the sample. The test schematic is divided 
into six subsystems for functional presentation: 
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Figure 3-2. 


Photo of Equilibrium Test Setup 


23 


(1) Pressure Control Section — Provides the vacuum source to maintain 
the system pressure and to allow drying of the system and 
sorbent sample. 

(2) Gas Supply System — Provides gas supply for controlled gas 
distribution to the sorbent under test. 

(3) Water Supply Sys tern -- P rov? des water vapor for controlled vapor 
distribution to the sorbent under test. 


( 4 ) Pressure Measurement Secti on- - Provides accurate system pressure 
measurement. 

(5) Temperature Control Section --Provi des accurate temperature control 
and monitoring of the sorbent temperature. 

(6) Electrobalance System -- Provi des accurate real-time weight 
measurement of the sorbent sample. 

The heart of the system is the electrobalance module, consisting of a 
Cahn RG electrobalance and its auxiliary equipment. This device is a beam 
balance with. a feed-back servo system designed to keep the beam level. The 
sensing of rotation of the balance beam is done optically. The optical 
sensor produces a signal which, after amplification, is directed to an electro 
magnet located at the beam fulcrum. The electromagnet produces a torque which 
restores the balance to a level position. The electromagnet current is 
accurately measured and translated into a weight reading. 

At the start of a test, the balance is zeroed by using calibrated weights 
on the reference weigh pan. After this, all weight readings are obtained 
remotely. Thus, there is no need for mechanical feed-throughs and the balance 
can be placed enti rely within the test envelope. The balance is insensitive 
to pressure levels within the apparatus. Because of its excellent sensitivity 
and the ability to produce accurate wei ght readings when subjected to the 
adsorption environment, the Cahn electrobalance is an excellent instrument for 
equilibrium adsorpti on- capaci ty measurements. Throughout this study, it has 
been very reliable, and repeatability of its readings has been excellent. 

The measurement of pressures in the range of interest in this study, 
generally 760 mm Hg down to less than I p, Hg, has traditionally been difficult 
and has hindered the production of accurate adsorption equilibrium data. For 
example, McLeod gages cannot be used with water vapor in certain pressure 
regions due to the condensation of the entrapped vapor in the gage capillaries 
Recently, new wi de- range pressure measurement instruments have been devised. 
For this study, a MKS Baratron pressure meter was used. Its sensing trans- 
ducers involving variable capacitors are linear with pressure and are quite 
stable. 
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Sample temperature was measured by a chromel-al umel thermocouple sus- 
pended very close to the sample weigh pan. The thermocouple did not actually 
touch the sorbent or the weigh pan as this would have caused erroneous weight 
readings. Sorbent temperature was controlled by immersing the glass sample 
tube in a temperature controlled bath, or by affixing an electric heater 
around the sample tube. Originally, the thermocouple reference junction was 
located within the upper portion of the gas-tight glassware of the test setup. 
The junction was considered to be at room temperature. Later, a feed through 
arrangement was made and the reference junction was moved outside of the 
glassware into an ice bath. 

Table 3-1 provides a list of the test equipment (Referenced by item 
number to Figure 3-1) used and briefly describes what function the equipment 
satisfied, and, where applicable, its tolerance and range. All active lines 
and valves of the setup are Pyrex glass. Equipment which exhibits off-gassing 
characteristics at low pressures or high temperatures was not used. The test 
setup (Figure 3-1) will maintain a system pressure of 0.02 p, Hg. 


Test Procedure 

The test procedure for producing equilibrium sorption capacity data is 
relatively simple and straight forward, but the preparation and maintenance of 
the equipment prior to actually taking equilibrium data is more subtle in 
nature. To produce accurate data, the system must be thoroughly dry and leak 
tight, and must be maintained in this condition. Many laboratory hours were 
spent in producing and verifying a leak tight system. 

The first step in the equilibrium procedure involves zeroing of the 
electrobalance. This involves the use of a calibration weight (60 mg) on the 
weigh pan. Then a sorbent sample of approximately 65 mg is weighed out on 
a laboratory balance and placed onto the electrobalance weigh pan. To avoid 
picking up of any considerable amount of water from the atmosphere, this 
operation is done as quickly as possible. The system is then closed off, and 
all portions of the system, except the gas supply system, are pumped down, in 
order to dry out the system. At the same time, water and other sorbates are 
removed from the sample by means of a vacuum-bakeout. The bakeout, usually at 
600°F, is carried out by placing the electric heater around the glass tube 
which surrounds the sample weigh pan. Some bakeouts were conducted at 400°F, 
however 600° F was adopted as' the standard bakeout temperature; this will be 
discussed at greater length later. 

Bakeout of the sorbent sample is continued until a stable system pressure 
is obtained, and the Cahn electrobalance shows no weight change taking place.. 
The stabilized pressure is usually less than I pHg. Once the bakeout is suc- 
cessfully completed, the sorbent is allowed to cool to the desired temperature, 
and provisions are made to maintain this temperature with either a heater or a 
bath. When the desired sorbent temperature is established, still at vacuum 
conditions, the weight of sample as given by the Cahn electrobalance is read; 
this weight is adopted as the base or dry weight of sorbent. Usually there is 
little difference between this weight and that measured at the bakeout tempera- 
ture. However, because of the possibility of thermal currents or temperature 
gradients which might produce inaccurate readings, the hot weights are not used. 
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TABLE 3-1 


TEST EQUIPMENT 


General Description 
and Model Identification 

Pressure Control Section 

Weigh Vacuum ump, Model 1402 

Mercury iffusion Pump 

Fab co Dewar 

Fabco Dewar 

Pyrex Valves, V- type 

Gas Supply System 

Fabco Dewar 

Wallace and Tiernan Pressure 
Gage, R/N 44S015 

Hoke Qas Plenum 
Valves 

Water Supply System 

Walch Vacuum Pump, Model 1405 

Fabco Dewar 
Water Reservoir 
Pyrex Valves, V-type 
Pressure Measurement Section 

MKS Baratron Pressuremeter, 

Type 77 R/N 44S946 

MKS Baratron Pressure Heads, Type 77 

McLeod Gauge Type GM1O0-A 

Ion Tube IG 100/C 

Fabco Dewar 

Valves 

Pyrex Valves, V-type 

Temperature Control Section 

Leeds Northrop Temperature 
Bridge S/N 8690, R/N 43BI79 

Reference Thermocouple 

Constant Temperature Bath 

Electrobalance System 

Cahn Electrobalance S/N 10307 

Cahn Electrobalance Control R/N 46B009 

Leeds Northrop Speedomax R/N 430296 


Rema r ks 

Controls test system from 0.02 to 760 mm Hg 

Controls test system from 3 to 760 mm Hg 

Controls test system from 0.02 to 3 mm Hg 

Thermal sink 

Thermal sink 

Selects system route 

Provides gas to sorbent 

Thermal sink 

Monitors gas supply pressure 

Stores gas for test use 

Fill and system supply valves 

Provides water vapor for the test 

Provides low pressure source for 
vaporizing water 

Thermal sink 

Water supply for system test 

Activates water distribution 

Provides accurate system 
pressure measurement capability 

System pressure readout 

Transducers for item 4.1 
Reference system pressure readout 

Thermal sink 

Used for calibration of MKS heads and 
system moni tori ng 

System valves 

Maintains and monitors sorbent temperature 
Readout for temperature 

Reference for temperature measurement 
Temperature control for sorbent 
Provides weight monitoring 
Electrobalance devise 
Controls for weight output 
Chart recorder for weight change 







Then the test gas, from the gas supply system is introduced to the system 
at the pressure required for the first datum point. The system conditions are 
maintained constant until the pressure stabilizes and the Cahn electrobalance 
indicates no weight change taking place. The sorbent gas loading is then 
determined by the change from the bakeout dry weight. 

Additional sorbate gas is introduced and the procedure is repeated for a 
higher pressure. In this manner, any number of data points can be obtained, 
thereby producing one adsorption isotherm. The bakeout procedure is then 
repeated and another isotherm can be generated. The developed isotherms are 
then plotted and, where applicable, certain tests of consistency are made. 

The test procedure for coadsorption equilibrium tests is similar, except 
that a sorbent pre-adsorption or preloading procedure is performed. For tests 
involving water and another sorbate water vapor, produced by the water supply 
system, is first distributed to the sorbent at a specific temperature and 
pressure. When the system reaches pressure and weight stability, the percent 
of water loading is noted. Then the gas supply system is opened to the system 
and the second sorbate gas is distributed to the sorbent. The partial pressure 
of the second sorbate gas is assumed to be equal to the total pressure minus 
the initial water vapor partial pressure. This assumption is valid as long as 
none of the initially adsorbed water is desorbed by the gas and no leakage 
into the system occurs. The condition of no water unloading was verified on 
a sample test by making a complete mass balance of the gases in the equilibrium 
apparatus; this involved the measurement of internal volumes and system pres- 
sures. Once the system reaches pressure and weight equilibrium, the co- 
adsorption equilibrium data can be obtained and plotted. The developed iso- 
therms are presented later in this section. 

Experimental Problem Areas 

There are several areas which, because of procedures or limitations of 
test equipment, may have caused inaccuracy in the data. The major problem 
areas are discussed below. 

Sorbent Bakeout Temperature . -Previous to this contract, 600° F was used as 
the vacuum-bakeout temperature prior to adsorption testing. This temperature 
also seems to be widely used in industry. However, at the start of this pro- 
gram, it became apparent that spacecraft systems would not be able to produce 
bakeouts at 600°F. Instead, 400°F was considered more likely. Therefore, in 
an attempt to produce more-di recti y-useable engineering data, 400°F was adopted 
for vacuum bakeouts. After some testing, it became apparent that 400° F bakeouts, 
even with CO 2 gave equilibrium data that was lower than previous AiResearch 
data, or that was available in the literature. Some special comparison tests 
were performed; the results of these tests indicated that 600°F bakeouts should 
be used. Therefore, for all subsequent tests 600 °F bakeouts were performed. 
However, not all 400°F-bakeout tests were rerun and some of these data are 
presented in this report. As such, these data are clearly labelled as to 
bakeout temperature. 
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Sorbent Temperature Measurement . ^Earl y in the program, measurement of 
sorbent temperature caused considerable problems. The problem area developed 
because there was no way provided with the Cahn electrobalance electrical feed- 
through to pass the chromel-al umel thermocouple wire out of the apparatus with- 
out creating addition thermocouple junctions. Therefore, a reference junction 
was setup inside the apparatus. It was located near the Cahn electrobalance, 
and was supposed to be at room ambient temperature. However, due to its proxi- 
mity to the Cahn-balance light source, its temperature was actually unknown and 
varied with system pressure. Later, this problem was completely solved when 
the junction was removed and a bulkhead seal was provided which allowed the 
thermocouple wires to pass directly out of the apparatus. A 32°F reference 
junction was then provided exterior to the system. 

Some of the early data has scatter in it which is undoubtedly due to 
sample temperature-measurement inaccuracy. Although calibrations were made in 
all cases with the early thermocouple system, it is believed that the calibra- 
tions were not sufficient to account for all modes of operation. Therefore, 
some measurements are suspected of being incorrect. As much as possible, the 
data was edited and only reliable data have been presented. 

Leakage .- Leakage of atmospheric gases into the apparatus was always a 
prime concern. Depending upon the test conditions and the sorbent loading at 
the time a leak develops, the measurable effects of a leak are not readily 
predictable. Certainly if a sorbent is nearly saturated with water, a leak 
would be noticed with increase in system pressure, without apparent weight 
increase. However, at low pressure levels and low loadings with sorbates other 
than water, leakage may not be noticeable from system pressure readings, but 
unusually high loadings (due to water vapor adsorption) might result. Therefore, 
great care and diligence on the part of test personnel was exercised in noting 
trends of the readings, and in frequent leak checks. 

Slowness of Adsorption in Low Pressure Ranges . -Under some conditions 
of very low pressure (in the range of 50 jj, Hg and below), the driving force 
for adsorption or desorption is very low and relatively long periods of time 
are requi red to achieve equilibrium. Sometimes it was difficult to determine 
whether equilibrium had been achieved or if system conditions were varying 
slightly. Therefore, there is some question concerning all data in the very 
low partial pressure region. One effect of slow adsorption is that data seem 
to be slightly low when data points are taken for increasing pressures; that is, 
when gas is added to the system to produce the next adsorption condition. On 
the other hand, data seems to be slightly high when tests were taken with 
decreasing pressure levels; that is, by pumping gas out of the system. Due to 
the time involved, not all systems received as much investigation in the low 
pressure region as would be desired. 

Reversibi 1 i ty of Adsorption .- It is noted that earlier during equilibrium 
testing there was concern of reversibility of adsorption and degradation of the 
sorbent itself due to prolonged exposure to the testing environment. Special 
tests as well as continued data comparisons indicated that adsorptions were 
reversible and the stability of the sorbents (molecular sieves and silica gel) 
was excellent. This was concluded from repeatability of dry-weight measurements 
after bakeouts, and from equilibrium data repeatability. 
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SINGLE- SORB ATE TEST DATA 


Table 3-2 presents the conditions under which mono-gas si ngl e-sorbate 
equilibrium sorbent capacity data were taken. The gas/sorbent combinations 
studied are shown in the table along with sorbent temperature and pressure 
ranges. Figures 3-3 through 3-22 present the s i ngl e-sorbate equilibrium 
data. 


TABLE 3-2 

TEST MATRIX FOR SINGLE-SORBATE EQUILIBRIUM DATA ON MOLECULAR SIEVES 


Mo l ecu 1 ar 
S i eve 
Sorbent 

Manufacturer 

H 2° 

C0 2 

°2- 

N 2 

5A 

Linde 

77, 167, 120, 
392° F, 

0- 1 0 mm Hg 

Testing from 
AAP program 
used 

32, 70° F 
0-760 mm Hg 

32, 68, 1 20° F 
0-760 mm Hg 

5A 

Davi son 

1 20°F 

0-10 mm Hg 

70, 1 20° F 
0-760 mm Hg 



I3X 

L i nde 

70° F, 

0-10. mm Hg 

32, 70, 120 
1 20° F, 

0-760 mm Hg 

32, 72, 

1 20°F, 

0-760 mm Hg 

32, 72, 1 20° F 
0-760 mm Hg 

1 3X 

Davi son 

70° F, 

0- 1 0 mm Hg 




I0X 

Linde 

70°F, 

0-10 mm Hg 




4A30 

Linde 

70° F, 

0-10 mm Hg 

70° F, 

0-760 mm Hg 



5A30 

i 

Li nde 


70° F, 

0-760 mm Hg 



4A 

Li nde 



70° F, 

0-760 mm Hg 

70°F, 

0-760 mm Hg 
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Carbon Dioxide Equilibrium on Molecular Sieves 

Figures 3-3 through 3-9 present carbon dioxide isotherms obtained with 
Davison 5A, Linde I0X, Linde I3X, Linde 4A30, and Linde 5A30 molecular-sieve 
sorbents. Linde 5A data is also presented. However it has been included from 
previous Apollo Application Program (AAP) test data. The 4A30 and 5A30 
sorbents, also identified as 4AXW and 5AXW respectively, have been included 
because the manufacturers performance predictions indicated an approximate 
40 percent increase in capacity over the standard 4A and 5A sorbents. 

The basic difference in the two materials is that 4A30 and 5A30 use a 
different, supposedly improved binder, with respect to C0^ adsorption 

capacity. The test data does not support this contention; the capacities are 
essentially those of standard 4A and 5A. 

Figure 3-3 presents 70° and I20°F isotherms for the Davison 5A sorbent. 

Also included (phantom lines) is the curve for Linde 5A at 70°F (AAP data). 

A variance can be seen between the two curves; therefore the 7mm Hg point was 
re-run for the Davison 5A sorbent and is shown with a large diamond. This 
point coincides with the curve for Linde 5A, This might be expected due to 
the fact that both sorbents are basically the same crystal structure. There 
could, of course, be differences in the amount of binder used. This re-run 
casts some doubt on the original AiResearch isotherms for this material. 

Figure 3-4 presents the data on Linde 5A from the AAP program. Isotherms 
of 0°, 40°, 60°, 100°, and I20°F are included. Actual data points are indicated 
by circles; other isotherms were i nterpolated. 

Figure 3-5 presents a 70°F isotherm for the Linde I OX sorbent. Figure 
3-0 presents a tabulation of various isotherms obtained for Linde I3X sorbent. 
Reference isotherms obtained after various bakeout temperatures of 400,500, and 
600°F are also included. Figure 3-7 presents a summary of the Linde I 3X 
sorbent test data for CO^ gas . Included on the figure are 32°, 0°, 70°, 90°, 

and I20°F isotherms. This summary figure includes data taken in the AAP 
program. Figure 3-8 presents the 70°F isotherms obtained for the Linde 4A30 
and 5A30 sorbents. 

An interesting observation of the 4A30 and 5A30 absorption characteristics 
is that a cross over point exists at approximately 12.5 mm Hg and 8.8 percent 
loading. (This is also true with the standard 4A and 5A sorbents.) The 4A20 
sorbent has a greater potential load capacity at low pressures and the 5A30 
has a greater capacity in the higher pressure range. 

Obtaining equilibrium data for 4A30 was extremely tedious due to the slow 
adsorption and desorption rates. Figure 3-9 presents a transient desorption 
curve plotted for both 4A30 and I3X. The desorption was obtained by simply 
opening the apparatus to the vacuum system. 
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Figure 3-3. Carbon Dioxide Equilibrium Data for Davison 5A Sorbent 
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Figure 3-5. Carbon Dioxide Equilibrium Data for Linde I OX Sorbent 










34 


Figure 3-6. Carbon Dioxide Equilibrium Data for Linde I3X Sorbent 
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Figure 3-7. Carbon Dioxide Equilibrium Data for Linde I3X Sorbent (Summary) 
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Figure 3-8. Carbon Dioxide Equi fi bri um Data for Linde 4A30 and 5A30 Sorbents 
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Figure 3-9. Desorption of Carbon Dioxide from Linde 4A30 
and I3X Sorbents 
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Water Vapor Equilibrium on Molecular Sieves 

The following figures present the water vapor isotherms obtained with 
Davison and Linde 5A, Linde I OX, Davison and Linde I3X, Linde 4A30 and Linde 
5A30. Figure 3-10 presents a comparison of a 70°F isotherm taken from the 
Ai Research AAP program and Linde data sheet 5A- I , a least square fit of the 
Linde isotherm at 70°F, and a 70°F isotherm with Davison 5A. Figure 3-11 
presents a comparison of the Davison and Linde 5A sorbent at l;20°F. The 
Davison curve was run by Ai Research and the Linde curve was taken from Li ride 
data sheet 5A-I. Figure 5-12 presents a 70°F isotherm with Linde I OX sorbent. 
Figure 3-13 presents a comparison of the davison and Linde I3X sorbent at 
70°F. The Linde data was again taken from the data sheet and the Davison 
data produced here. Figure 3-14 presents a family of isotherms at 33°, 60°, 
70°, 90°, 120°, 208°/ 302°, and 400°F for Linde I3X sorbent. Figure 3-15 
presents 70°F isotherms for Linde 4A30 and 5A30 sorbents. 

Comparison of all the previous listed data indicate that the I OX sorbent 
has a slightly higher loading capacity at higher pressures, i.e., from 10 mm 
Hg down to 0.1 mm Hg. At approximately 0.1 mm Hg the two capacities’ are the 
same and for lower pressures the Davison 5A indicates the greater load 
capacity. Low pressure loading was extremely difficult to obtain due to the 
time required to stabilize. 

Oxygen Equilibrium on Molecular Sieves 

The following figures present the oxygen isotherms obtained with Linde 
4A, 5A, and I3X sorbents. Figure 3-L6 presents a 70°F isotherm for the 
Linde 4A sorbent. Figure 3-17 presents 32° and 70°F isotherms for the 
Linde 5A sorbent, and Figure 3- 18 presents 32°, 72°, and 120°F isotherms for 
the Linde I3X sorbent. The I20°F isotherm for the Linde I3X sorbent was 
difficult to obtain in both the high pressure and low pressure regions; 
however sufficient mid-range data was obtained to demonstrate slope charac- 
teristics compatible with the 32° and 72°F isotherms. Oxygen loading 
generally was quite small making data difficult to obtain. 

Comparison of the 70°F isotherms indicate very little to choose between 
in the three sorbents. However 4A and 5A do appear to possess a slightly 
higher potential load capacity at high pressure. 

Nitrogen Equilibrium on Molecular Sieves 

The following figures present the nitrogen isotherms obtained with Linde 
4A, 5A, and I3X sorbents. Figure 3-19 presents a 70°F isotherm for the 
Linde 4A sorbent and Figure 3-20 presents three isotherms for 32°F, 68°F, and 
I20°F, for the Linde 5A sorbdnt. Also included on this figure is a referenced 
nitrogen isotherm obtained from ESSO (Reference 3) The basic slope charac- 
teristics of the Ai Research data and the referenced isotherm are in agreement. 
Figure 3-21 also presents three isotherms, for 32°F, 72°F, and I20°F, for the 
Linde I3X sorbent. The Linde 5A appears to have the greatest potential load 
capacity of the three materials at both high and low pressure levels. This 
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Figure 3-10. Water Vapor Equilibrium Data for Linde and Davison 5A Sorbents 
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Figure 3-12. Water Vapor Equilibrium Data for Linde lOX Sorbent 
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Figure 3-15. Water Vapor Equilibrium Data for Linde 4A30 and 5A30 Sorbents 
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Figure 3-1 6. Oxygen Equilibrium Data for Linde 4A Sorbent 
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Figure 3-19. Nitrogen Equilibrium Data for Linde 4A Sorbent 
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Figure 3-20. Nitrogen Equilibrium Data for Linde 5A Sorbent 
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Figure 3-21. Nitrogen Equilibrium Data for Linde 5A Sorbent 



observation is based on an approximated 70°F isotherm. The 4A and I3X are 
essentially identical, within experimental repeat! bi 1 i ty, at low pressure 
levels; however the 4A has an increasingly higher load capacity from approxi- 
mately 200 mm Hg on up to 760 mm Hg. 

Water Vapor Equilibrium on Silica Gel 

A test series was performed on Silica gel to determine the water vapor 
adsorption characteristics. Two isotherms were plotted and are included in 
Figure 3-22. 

The silica gel used on the test was Davison Grade 05 in the form of a 
6-15 mesh irregularly-shaped particles. The isotherms were plotted from 
test data at 60°F and I20°F. A 90°F isotherm was then interpolated from 
the test data and is also included on the figure. Prior Davison data for 
60°F and I20°F are also included in phantom lines in the figure. Major 
differences are obvious in the low-pressure-range isotherm slopes, as well 
as in the high pressure termination slopes. Because of this discrepancy 
the Davison Chemical Division of W. R. Grace and Company was contacted. 
Apparently the Davison data are almost 15 years old and the silica gel 
manufacturing process has been continuously improved. An important charac- 
teristic of the silica ge 1 adsorpt ion characteristics can be observed 
on the 60°F isotherm where the slope of the loading curve decreases 
rapidly above 6 ran Hg pressure. At II mm Hg the loading is 41.6 percent. 

This compares with the calculated maximum loading for silica gel of 43 per- 
cent assuming that the maximum loading occurs when intra-crystalline voids 
are full of liquid water. 

The 60°F and 120°F isotherms generated for this program provide very 
close correlation with the associated dynamic test data. 

COADSORPTION TEST DATA FOR MOLECULAR SIEVES 

It is common in the mass transfer field to consider adsorption processes 
on the basis of one primary sorbate, and to ignore the simultaneous sorption 
of other gases. Primarily this is done for simplicity; to do otherwise might 
possibly involve cons i derable. compl icat i ons . In many instances the actual 
process can be characterized well on the basis of one sorbate; i.e., the 
drying of a gas stream by molecular sieves. Here, the water equilibrium data, 
which was taken without the carrier gas being present, is considered to be 
valid, wi th no effect due to the presence of the large amount of carrier gas. 
The mass-transfer dymanics are similarly characterized without considering 
the sorption of the carrier gas. If the carrier gas sorption is considered, 
the assumption is usually made that it has no effect on the sorption of water. 

There are, however, many cases where there is appreciable adsorption of 
more than one material, and as a consequence, interactions take place which 
alter both the dynamic and the equilibrium character of the process. Space- 
craft CO removal systems using molecular sieves certainly fall into this 
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category. The poisoning effect of coadsorbed water on the C0 2 capacity of 

molecular sieves has been well known for some time. However, until lately, 
little quantitative information was available. 

In this program it was desired to produce quantitative data on the 
equilibrium capacity of molecular sieves for a particular sorbate in the face 
of the coadsorption of another sorbate. Obviously, the most important coadsorp- 
tion sorbate pair for which data were sought is CO^ and water, i.e., the poi- 
soning effect due to water. This was not the only coadsorption pair considered. 
Other pairs of interest can be seen in Table 3-3. 

The coadsorption pairs involving either nitrogen or oxygen were of interest 
for two reasons. First, in most instances where these gases are adsorbed, 
they will eventually be lost from the vehicle. Therefore they represent a 
definite weight penalty to the CO^ removal system; and of course provisions 

would have to be made within the spacecraft to store makeup gases. It was 
desired, in this program, to ascertain any possible coadsorption effect, 
especially with water, which might appreciably alter the loss of atmosphere 
gases. The second reason for the oxygen and nitrogen equilibrium data is to 
observe if the large volumetric preponderance of these gases has an effect on 
the equilibrium capacity for CO^ and water. To a great extent^ the answer to 

this latter inquiry can be obtained from the s i ng le-sorbate dynamic tests 
contained in Section 4. From these tests it is concluded that the equilibrium 
capacity data as taken with only the sorbate (either CO^ or water) in contact 

with the molecular sieve is not noticeably affected by the presence of the 
nitrogen carrier gas. 

The following figures present the results of equilibrium sorbent capacity 
testing for conditions where a preload of one sorbate is produced followed 
by adsorption of another gas. The sorbate pairs for which tests were conducted, 
are shown in Table 3-3. Also included in the table are the sorbant temperature 
and pressure range. 


Coadsorption of CO^ with H^O Preload 

The figures 3-23 through 3-26 present the coadsorption test data for carbon 
dioxide and water vapor. It can readily be observed that for coadsorption of 
C0,j with a H^O preload the CO^ capacity is significantly reduced, however, 

even for a water load of 11.9 percent, some carbon dioxide was still adsorbed. 
Figure 3-23 presents a family of 70°F isotherms at various water vapor pre- 
load percentages, 0, 1.89, 4.02, 6.43, and 12 percent. The 12 percent curve 
was run twice to verify repeati bi 1 i ty ; excellent agreement was obtained. Figure 
3-24 presents a family of curves for constant water preload (approximately 
6.1 percent) at different temperatures, 32°, 70°, 120°, and 200°F. It is 
obvious that as the sorbent temperature is increased less and les!s adsorption 
sites for CO adsorption are available. Figure 3-25 presents the coadsorption 
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COADSORPTION TEST MATRIX 
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Figure 3-23. Coadsorption Equilibrium Data for Water Vapor and 
Carbon Dioxide on Linde 5A Sorbent 
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Figure 3-24* Coadsorption Equilibrium Data for Water Vapor 
Carbon Dioxide on Linde 5A Sorbent 
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Figure 3-25. Coadsorption Equilibrium Data for Water Vapor and 
Carbon Dioxide on Linde I3X Sorbent 
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Figure 3-26. Coadsorption Equi librium Data for Water Vapor and 
Carbon Dioxide on Linde I3X Sorbent 


data obtained with Linde I3X sorbent. This figure contains coadsorption data 
at different water vapor preloads of 0, 2.07, 3.99, 5.9, 11.9, and 17.0 
percent for the 70 or 72°F temperatures. 

Figure 3-26 presents a family of isotherms at approximately the same water 
vapor preload for the Linde I3X sorbent. The water vapor preload is approxi- 
mately 6 percent and the isotherms are 32°, 72°, 120°, 200°F, and a single point 
at 300°F. 

It appears that even though the Linde 5A sorbent has a greater C0 2 capacity 

than I3X, with no water vapor present the converse is true when a water vapor 
preload is incorporated. 

Coadsorption of Nitrogen With Carbon Dioxide Preload 

Figure 3-27 and 3-28 present the test data for the Linde 5A and I3X 
sorbent, respectively, when subjected to nitrogen gas after preloading with 
carbon dioxide. Isotherms at 70°F was taken at different C0 2 preloadings. The 

I3X data taken was not acceptably repeatable, however the basic curves and the 
re-runs are presented. 

For the Linde 5A sorbent, with a 3 percent C0 2 preload, the nitrogen 

adsorption was slightly higher than that obtained without a preload. The 
tests results indicate that the reduction in nitrogen adsorption is negligible 
up to 3 percent C0 2 preload. Above 3 percent C0^ preload, negligible adsorption 

occurs until the nitrogen partial pressure exceeds 100 mm Hg. 

In previous tests, the introduction of C0 2 to a water preload molecular 

sieve caused the weight to continuously increase until an equilibrium weight 
was achieved. In these tests, the introduction of nitrogen on a CO^ preload 

sieve caused the weight to quickly increase above the final equilibrium value 
and then slowly decrease to the ultimate value. The initial increase was greater 
than that observed for pure nitrogen adsorption on the sorbent. This initial 
increase is due to local increase in C0 2 partial pressure in the vicinity of 

the sorbent, due to the imperfect mixing of the two gases as the nitrogen is 
introduced. Since the partial pressure of CC> 2 for these tests is several mm 

of Hg, the rate of adsorption is quite high, and the initial additional CO 

adsorption is greater than the final nitrogen adsorption. The weight indication 
drops slowly as the gases mix and finally reach an equilibrium value. 

Coadsorption of Oxygen With Water Vapor Preload 

Figure 3-29 presents the data obtained for Linde 5A molecular sieve pre- 
load with water and subjected to oxygen at various pressures. A family of 70°F 
isotherms presented at various water vapor preloads, 0, 3.1, 5, and 8 percent. 
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Figure 3-27. Coadsorption Equilibrium Data for Carbon Dioxide and 
Nitrogen on Linde 5A Sorbent 








(luaqjos q | 


00 l/psqjospe usBcmiu qi) 


6ui peo[ u96oj^ i n 


61 


Figure 3-28. Coadsorption Equilibrium Data for Carbon Dioxide and 
Nitrogen on Linde I3X Sorbent 
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Figure 3-29. Coadsorption Equilibrium Data for Water Vapor and 
Oxygen on Li nde 5A Sorbent 



Figure 3-30. Coadsorption Equilibrium Data for Water Vapor and 
Nitrogen on Linde 5A Sorbent 





Coadsorption of Nitrogen With Water Vapor Preload 

Figure 3-30 presents the data obtained on the Linde 5A molecular sieve 
preload with water vapor and subjected to nitrogen at various pressures. A 
family of 70°F isotherms are presented at various water vapor preloads, 0, 
3, 5.1, and 8.1 percent. 


HEAT OF ADSORPTION 


The adsorption of gases on solids is accomplished with the evolution of 
heat, which varies with the quantity of gas adsorbed, the characteristics of 
the solid, and the adsorption temperature. The heat given off is called the 
heat of adsorption. This heat varies with the amount of material adsorbed; 
vvhen evaluated for an infinitesimally small sorbent loading change, it is 
called the differential heat of adsoprtion. As a lower limit, as the amount 
of gas adsorbed increases the differential heat of adsorption approaches the 
heat of liquefaction. Differential heats of adsorption can be calculated from 
the Claus ius-C lapeyron equation. Usual ly this involves obtaining isotherms 
(equilibrium pressure versus temperature for constant adsorption weight load- 
ings) by plotting the natural logarithm of the pressure against the reciprocal 
of the sorbent temperature. The slope of the isotherms is the heat of 
adsorption divided by the universal gas constant. 

Figure 3-31 represents several isotherms plotted on semi-log coordinates 
for COg on Linde I3X sieve (the data for this plot are from Figure 3-8). The 

differential heats of adsorption obtained by this procedure are plotted in 
Figure 3-32. It should be noted that at zero loading the maximum exothermic 
energy release occurs and as adsorption takes place the energy decreases. The 
equation used for this calculation is obtained from the- Claus i us-Clapeyron 
equati on : 


where 



Ah = differential heat of adsorption, Btu/lb 

P| = equilibrium pressure at specific weight loading for 
1 low temperature isotherm, mm Hg 

P^ = equilibrium pressure at specific weight loading for high 
temperature isotherm, mm Hg 

R = Universal gas constant, 1.987 Btu/lb-mole x °R 

W = Molecular weight 

T| = Temperature at P , °R 

T^ = Temperature at P^, °R 
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Figure 3-32. Differential Heat of Adsorption for Carbon Dioxide 
on Linde I3X Sorbent as a Function of Loading 
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Table 3-4 presents some typical heat adsorption values calculated at nominal 
loadings of the sieve. The expanded data points shown for CO^ at various 

loadings indicate the physical trend of the heat of adsorption. Where only 
two isotherms are available calculations were made because the slope's could 
not be sufficiently guaranteed. 


TABLE 3-4 

DIFFERENTIAL HEATS OF ADSORPTION 


Gas 

Molecular 
S i eve 
Type 

Weight 
Loading, f> 

Differential Heat of Adsorption, 
Btu/lb 

C °2 

LI3X 

1 

500 



2 

476 



2.5 

445 



3 

403 



4 

375 



6 

358 



8 

344 



12 

330 

CM 

O 

O 

L5A 

4 

416 

H 2° 

LI3X 

10 

1 1 I 1 

°2 

LI3X 

. 1 

245 
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SECTION 4 

SINGLE-SORBATE DYNAMIC MASS-TRANSFER DATA 


OBJECTIVES AND SCOPE 

The general objective of the s i ng 1 e-sorba te dynamic mass-transfer test 
series was to produce information that would allow a better understanding and 
characterization of molecular-sieve sorbents under dynamic sorption conditions 
typical of spacecraft systems. The method used to produce this information was 
to conduct si ngl e-sorbate breakthrough runs under controlled, constant condi- 
tions. Such testing is a standard tool in evaluating the performance of packed 
beds. The testing is carried out in two steps. The first step is the estab- 
lishment of initial bed conditions of temperature and sorbent loading which 
are uniform and can be accurately determined numerically. For example, the 
usual bed conditioning prior to an adsorption run is a vacuum bakeout to 
remove all traces of sorbates, followed by cooling to the desired adsorption 
temperature. In the second step, a gas stream containing the sorbate of 
interest is passed through the bed until the bed has adsorbed as much of the 
sorbate as is possible under the test conditions; that is, until the bed is 
saturated or equilibrated. During the transient sorption period, all inlet 
and boundary conditions (flow rate, pressure levels, inlet concentrations, 
bed temperature) are held constant and all parameters are monitored, especially 
inlet and outlet sorbate concentrations. 

Several specific objectives and goals were defined for the breakthrough 
tests : 

o The most important goal was the evaluation of mass-transfer coeffi- 
cients and intraparticle di ffusi vi ties for various representative 
operating conditions. These basic parameters were obtained from the 
breakthrough data with the help of the performance-predi cti on pro- 
grams. With this data properly correlated, the programs can be 
used with confidence for a wide variety of bed geometries and process 
condi tions. 

© Qualitatively, from the shape of the plotted breakthrough curve 

(the outlet concentration plotted vs time), it was desired to com- 
pare the performance of various sorbents under similar operating 
conditions, and to compare the performance of a single sorbent under 
various adsorption conditions. 

o Also in a qualitative sense, it was desired to compare the thermal 
performance of the different heat-transfer surfaces which were 
incorporated into the test beds. This was to be accomplished by 
observing the ability of the heat-transfer surfaces to maintain 
an isothermal bed during the exothermic adsorption period. 
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© By integration techniques to yield the amount of sorbate left on 

the bed, an equilibrium datum point, corresponding to the final bed 
temperature and inlet partial pressure, can be obtained. This was 
desired as a check on the data produced from the equilibrium test- 
ing, or obtained from the literature. 

o From the breakthrough data, removal efficiency of the sorbent bed 
(cumulative fraction of entering sorbate which is retained by the 
bed) as a function of time can be obtained. This is desired to 
properly size spacecraft beds and to set cycle times. 

Table 4-1 presents the sorbent/ sorbate combinations for which dynamic 
adsorption breakthrough tests were made. For the potential predrier-bed 
sorbents, Linde I 3X and silica gel, both adsorption breakthrough tests and 
gas strips were performed. Table 4-2 presents conditions and ranges for 
which breakthrough tests were performed. Dry nitrogen was used as the carrier 
gas in all dynamic tests. Obviously, tests could not be run for all possible 
test conditions; however, it is believed that sufficient ranges were included 
to identify trends and to permit mass-transfer parameter evaluations. 


TABLE 4-1 

SORBENT/SORBATE COMBINATIONS USED IN SINGLE-SORBATE DYNAMIC 

MASS-TRANSFER TESTS 


Sorbent 

Type 

Sorbent 

Manufacturer 

Form 

and Size 

Sorbate in 
N i trogen 
Carrier Gas 

Silica gel 
(Grade 05) 

Davi son 

Granul es, 

6 to 1 6 mesh 

Water 

I3X 

Li nde 

Pel 1 ets, 

1 /I 6-i nch 

Water 

co 2 

5A 

Linde 

Pel 1 ets, 

1 /I 6-i nch 

co 2 

5A 

Davi son 

Spher i cal 
1 2 mesh 

beads, 8 to 

C°2 
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TABLE 4-2 


CONDITIONS FOR WHICH SINGLE/SORBATE BREAKTHROUGH 
TESTS WERE PERFORMED 


Sorbate 

Inlet 
Partial 
Pressure 
(and Dew 
Point) 

Total 
Press u re 

Inlet 

Concentration 

Bed 

Temperature 

Flowrate 

Water 

9. 2 mm Hg 
(50°F dew 
poi nt) 

5 psi a 
( 258 mm Hg) 
for si 1 i ca gel 

7 ps i a 
(363 mm Hg) 

0.024 gm water/ 
gm n 2 

0. 01 8 gm water/ 
gm N 2 

60° F 
70°F 
90° F 
1 20°F 

1.3 to 5. 4 
lb/hr 


7 mm Hg 


0.031 gm C0 2 / 
gm n 2 

50° F 


co 2 

1 , 75 mm Hg 

7 ps i a 
(363 mm Hg) 

(1,9 vo 1 $ C0 2 ) 

0.0075 gm C0 ? / 
gm n 2 

70° F 

1.3 to 5. 4 
lb/hr 




(0.45 vol % C0 2 ) 




EQUIPMENT AND PROCEDURES 


Bed Designs and Heat-Transfer Surfaces 

One of the objectives of this program has been to compare the effective- 
ness of various size heat-transfer surfaces placed within packed sorbent beds. 
This information is desired for those cases where bed cooling and heating is 
to be provided by the spacecraft heat-transport system. It is not obvious 
what type of surface will produce the optimum bed performance, due to a number 
of counteracting effects. For example, from only a heat-transfer point of 
view, it is desired to produce as nearly isothermal a bed as possible; thus, 
small, close-packed heat-transfer surfaces would be desired. However, with 
very dense heat-transfer matrices, a great deal of "wall effect" is to be 
expected. The wall effect is due to the fact that sorbent pellets do not 
pack closely in the region of the walls created by the plethora of fin sur- 
faces. This results in lower overall packing density (greater void fraction 
in the wall region), and lesser contact between gas and sorbent, often des- 
cribed by the names "channelling", "bypassing", or "fingering". 

Prior to this contract, AiResearch had developed a number of research 
and development sorbent systems which involved high-density heat-transfer; 
surfaces placed within packed beds. Basically, these beds were cross-flow 
plate-fin heat exchangers, with sorbent material packed into one set of 
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passages. Host of the previous expert ence was with rectangular straight fins, 
1/4 in. in height and spaced 1/4 in. apart. The results of tests with such 
heat-transfer surfaces showed less-than-desi red packed bed densities, and a 
tendency of premature breakthrough due to channeling. For this contract, 
three specially designed beds were fabricated, each with a different size 
rectangular, offset plate-fin surface. The sizes were 1/4 in., 1/2 in., and 
I in. In any bed, the fin height and spacing were the same. In addition, 
for these special beds the fin surfaces were all offset (or interrupted) in 
the flow direction by a dimension equal to the fin height. Thus, in the 1/2- 
in. -fin bed, the continuous fin length in the flow direction is only 1/2 in. 

At this point, the downstream fins are placed midway between the previous fin 
positions. Figure 4-1 gives a comparison of plain and offset rectangular fin 
surfaces. The fins are shown in their original form, before stacking into a 
heat-exchanger configuration and brazing. The offset fin was adopted because 
of the expectation of reduced wall effects (channeling) afforded by the 
interruptions of the fins. No effect on packed-bed density was expected. 

Figure 4-2 is an outline drawing typical of the beds used in this program. 

Figure 4-3 is a photograph of the 1/2- i n . -fin bed. To promote uniform distri- 

bution of coolant in the liquid passages of the beds, the fin surfaces were 
turned sideways to the normal flow direction. This forces the coolant to 

flow in and out around the fins. This coolant fin arrangement is shown in 

Fi gure 4-2. 

A few points should be noted concerning the choice of plate-fin surfaces 
for use in' sorbent beds when bed cooling and heating are desired. First, 
the manufacturing technology for the plate-fin heat-exchanger sorbent beds is 
well developed, and is essentially qualified for use in space, Such sorbent 
bed designs are only slightly different from many heat exchangers used in 
Mercury, Gemini, and Apollo environmental control systems. Second, as in 
many heat-exchanger applications, plate-fin designs offer tremendous flex- 
ibility in tailoring the ratio of heat-transfer surface area on one side to 
that on the other side. With heat exchangers, optimum designs are usually 
characterized by having the product, hA, of heat-transfer coefficient and 

area on one side approximately equal to that on the other side. With a fluid 

stream, such as a low-velocity gas, where low heat-transfer coefficients are 
developed, it is easily possible to provide conslderabl y more plate-fin surface 

area for the gas and to tailor the ratio of areas so that the hA products of 

the exchanger are nearly equal. With packed sorbent beds, since the sorbents 
are very poor conductors of heat, it is almost mandatory to provide much more 
heat-transfer area on the gas/sorbent side than on the coolant side. With 
plate-fin designs, it is quite easy to accomplish this, no matter what the 
desired bed area and length values may be. In addition, the plate-fin 
concept provides both very low coolant-side pressure drop (when desired), 
and very low liquid holdup within the bed. Further, the large amount of sur- 
face area can be incorporated uniformly throughout the packed sorbent particles, 
and the presence of these surfaces has very little effect on the gas-side 
pressure drop. Other means of providing bed cooling, such as cooling coils 
with, external fins, cannot compete with plate-fin surfaces in the several 
areas of concern as discussed above. For, small beds like those for space- 
craft applications, the additional cost of plate-fin sorbent beds is far 
outweighed by their many advantages. 
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a. Plain Or Straight Rectangular Plate-Fin Surface 



b. Rectangular Offset Plate-Fin Surface As Used In Test Beds 

S-41922 -A 

Figure 4-1. Rectangular Plate-Fin Heat-Transfer Surfaces for Packed Beds 




Figure 4-2. Test Bed for Dynamic Adsorption Tests 
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Figure 4-3. I / 2 «» in- *“Fin Adsorption Bed 
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A set of seven thermocouples was placed within each of the beds. One 
thermocouple was placed on the fin surface near the bed inlet and the other 
six were positioned longitudinally along the bed length, about I in. apart, 
midway between fins. 

Table 4-3 gives the dimensions of each bed and pertinent physical char- 
acteristics. Each bed had a 3- by 3-in. frontal face with a 6-in. flow 
length. Depending upon the sorbent and the method of packing, from 1.0 to 
1.25 lb of sorbent can be contained by each bed. Except for the packed-bed 
fins, which are copper, and the coolant fins, which are nickel, the beds are 
constructed of stainless steel. The bed heat exchanger core was vacuum 
brazed; the end flanges and pans were welded to the core. 

The bed packing procedure was initiated by a vacuum -bakeout of the 
sorbent contained in an approximately 2-qt sample bomb. Bakeout conditions 
were 500° to 600°F at low pressure, usually approaching 10 p, Hg. Packing of 
the bed was accomplished in a nitrogen dry box which maintained a dew point 
of -40°F or lower. During packing, the bed was clamped in a vertical posi- 
tion to an electric vibrator table. Sorbent pellets were slowly poured into 
the bed while it was being vibrated. The sorbent was held within the bed by 
means of 24-mesh screens. The bottom screen was permanently welded in place; 
after packing, the top screen was secured in place by small tack welds. 

Prior to packing, the bed and the bomb were weighed. After the bed was 
fully packed, it was again weighed while still in the dry box. Thi s gave an 
approximation of the mass of dry sorbent contained in the bed. Then, as much 
as possible of the excess sorbent which was spilled during the packing operation 
was swept up and placed back in the sample bomb. The bomb was then weighed to 
obtain another approximation of the dry sorbent mass. The two values for 
bed mass were usually in close agreement, but not in every case. The reason 
for the deviation between readings, and the reason these weighings were con- 
sidered to yield only approximate sorbate weights, is that molecular sieves have 
equilibrium capacities for nitrogen at 760 mm Hg of about I percent. Also, even 
at the dry box humidity conditions, molecular sieves can hold well over 6 pei — 
cent water. It is doubtful that there was enough water vapor available for any 
appreciable adsorption; however, nitrogen adsorption could not be neglected. 

The packed bed mass was more accurately measured at the end of a sequence 
of tests. At that time, the bed was carefully unpacked and its contents 
placed in a sample bomb. After a considerable vacuum bakeout, the bomb was 
closed off and, after cooling, was weighed. This weight was then used to 
obtain the dry sorbent weight. In most instances, this weight was not sub- 
stantially different from the values obtained at the time of packing. When 
there was a difference, the latter weight was used. 

During the program in some cases where a bed was packed two or more times 
with the same type of Linde molecular-sieve pellets, it was noticed that bed 
weights could be substantially different. It is considered that the primary 
reason for the differences lies in the method of packing. During packing, 
moderate vibration of the bed was employed. High levels of vibration were 
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TABLE 4-3 


TEST BED DIMENSIONS 



; Nominal Packed-Bed Fin Height, (in.) 
Bed Designation 

Physical Item 

I 

1/2 

1/4 

PACKED-BED/GAS SIDE: 
Bed width (in.) 

3.00 

2.98 

3.00 

Bed height (in.) 

3.25 

2 . 93 

2.95 

Bed length (in.) 

6.10 

6.02 

6.02 

Number of passages 

3 

5 

5 

Number of sandwiches per passage 

1 

1 

2 

Fin thickness (in.) 

0.008 

0.008 

0.008 

Fins per in. (in. 1 ) 

1 

2 

4 

Splitter thickness (in.) 

— 

— 

0.01 0 

Plate thickness (in.) 

0.010 

0,010 

0.010 

Fin cross-section area to gas (in. 2 ) 

0.155 

0.237 

0.382 

Coolant passage cross-section area to gas (in. 2 ) 

0.780 

1.192 

1.200 

Free flow area (in. ) 

8.82 

7.29 

7.27 

Plate area (in. 2 ) 

145.5 

207.8 

205.5 

Fin area (in. ) 

89.6 

176.8 

524.5 

Total heat- transfer area (in. 2 ) 

235.1 

384.6 

730.0 

Plate volume (in. 3 ) 

5.67 

5.59 

5.64 

Fin volume (in. ) 

0.95 

1 .43 

2.30 

Free bed volume (in. 3 ) 

53.95 

43.66 

43.75 

COOLANT SIDE: 




Passage width (in.) 

5.80 

5.80 

5. 80 

Passage length (in.) 

3.25 

3.25 

3 . 25 

Number of passages 

4 

6 

6 

Fin height (in.) 

0.050 

0.050 

0.050 

Fins per in. (in. *) 

20 

20 

20 

Fin thickness (in.) 

0.002 

0.002 

0.002 

Plate area (in. 2 ) 

144.8 

217.2 1 

217.2 

Fin area (in.) 

144.8 

217.2 

217.2 

Total heat transfer area (in. 2 ) 

289.6 

434.4 

434.4 

Coolant volume (in. 3 ) 

. . . i 

3.47 

5.21 

5.21 











avoided as this tended to eject a large fraction of the sorbent bellets that 
were being poured onto the face of the bed. As the program continued, it was 
found that after the initial filling, some additional settling could be 
accomplished by a sustained vibration at very high levels and more sorbent 
could be added. Thus, later packings tended to yield higher bed weights. 

There are other factors involved in the final bed weight. The length of 
the Linde pellets is of importance. Some batches of material seemed to have 
larger amounts of fairly long pellets (about 1/4 in.). It is suspected that 
these longer pellets do not pack as wel 1 as shorter pellets; at least not 
without considerable vibration. There can also be variation in the density 
of the basic materia] produced by the manufacturer. 

Packed bed density is discussed separately in Section 8. Packed 
densities of the test beds are presented, along with densities obtained in 
special containers constructed to determine the effect of fin size on packing 
densi ty. 


Dynamic Mass-Transfer Apparatus 

Figure 4-4 is a simplified schematic of the dynamic mass-transfer 
apparatus. Figure 4-5 is a photograph of the test equipment. The heart of 
the apparatus, the dynamic test bed, is located in the lower left side of the 
photograph of Figure 4-5. The apparatus was so constructed as to allow flow 
in either direction through the bdd, and also to be able to vacuum desorb 
from either or both ends of the bed. Consequently, there are dual sets of 
inlet and outlet valves on both sides of the bed. The two large vacuum gate 
valves flanking the bed are clearly visible in the photograph. Three valves 
are attached to the inlet manifold above the bed. Two of the valves are used 
to direct flow to the desired inlet of the bed. The other valve, mounted 
between the inlet valves, is a bypass valve. The back-pressure valves on the 
bed exits, and their connections to the large vacuum duct, are obscured in the 
photograph. 

Referring to the simplified schematic of Figure 4-4, the flow path in 
the apparatus is as follows. First, laboratory dry nitrogen, derived from 
large-scale cryogenic storage (<-80°F dew point), passes through a pressure 
regulator. This regulator outlet is maintained at 20 psia. The flow is 
directed through a bank of rotameters, each being calibrated at 20 psia. 

To introduce water vapor into the stream, all or a portion of the nitrogen 
flow can be directed through a water-filled bubbler. To enhance the humidifi- 
cation process, the bubbler is contained in a hot water bath. Next, the 
humidified nitrogen steam is directed to a condenser. In some of the experi- 
ments, the condenser was used to adjust the dew point of the gas stream. In 
other cases, the condenser was dormant. 

The gas emerging from the humidity-control condenser enters the bed 
inlet manifold, where flow can be diverted around the bed directly through a 
back-pressure valve to the vacuum system. This was the case while experimen- 
tal conditions were being established. For adsorption tests, the flow is 
directed through the bed, then through a back-pressure valve to the vacuum 
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Figure 4-4. Schematic of Dynamic Mass-Transfer Apparatus 
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Figure 4-5. Dynamic Mass-Transfer Apparatus 

system. Vacuum is provided by two large pumps. One pump is used when high 
flow rates and moderate vacuums are required, as in adsorption runs. The 
other vacuum pump is used when very low pressures are desired; for example, 
during vacuum bakeouts. 

Thermal control of the bed and of the humidity-control condenser is 
provided by a heat-transport circuit using 62. 5 -percent ethyl ene-gl ycol /water 
mixture. The water-glycol emerging from the reservoir-pump system is split 
into two separate streams. One stream Is directed through temperature- 
controlled electric heaters. The other stream is cooled by a Freon refrigera- 
tion system. By selective mixing of these streams at the bed, and also at the 
condenser, the desi red temperatures can be obtained. Thus, a bed can be 
heated to 120°F for an adsorption run, while the condenser can be held at 
50°F for humidity control. The water-glycol system is capable of temperatures 
from 20°F up to approximately 200°F. Water-glycol effluents from the bed and 
the condenser are returned to the pump reservoir. 

Above 200°F, bed heating can be accomplished by use of an open-loop 
nitrogen heating circuit. Separate electric heaters were installed for this 
purpose. To use the nitrogen heating circuit, the bed is first drained of 
water-glycol and then isolated from the liquid circuit. By use of appropriate 
valving and lines, the hot nitrogen can be directed into the heat-transfer 
passages of the bed. This system was actually used in only a few adsorption 
runs. Its main use is for vacuum bakeouts of the bed. In these instances, 
the nitrogen is heated to approximately 500°F. For bakeouts, both of the large 
gate valves are opened to the vacuum system. 
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The source for carbon dioxide used in the experiments is a high pressure 
bottle of very dry CO 2 gas. Through a dual stage regulator, CO 2 is supplied 
to a precision rotameter operating at 20 psia. CO 2 flow rate is controlled 
by a valve downstream of this rotameter. The C0 2 is mixed into the main 
nitrogen stream downstream of the humidity-control condenser. 

Humidification is provided by the heated bubbler. Two methods of 
humidity control are available. One of these methods has already been 
mentioned: the condensation of excess moisture in the humidity-control 

condenser. As experience was gained with the test apparatus, it was found 
that the inlet dew point was sometimes hard to control in this manner. This 
is because the condenser and the bed both use the same hot and cold water- 
glycol streams. During the period of time when the adsorption rate in the 
bed was nearly constant, conditions could be maintained fairly steadily. 

However, as the desorption process subsided, the water-glycol stream returning 
to the reservoir from the bed became cooler and cooler. This had the effect 
of gradually cooling the water-glycol flow to the condenser. Without contin- 
ual adjustment, bed inlet dew point tended to decline. In addition, over 
humidification at the bubbler caused the rate of extraction of water from it 
to be relatively large, requiring frequent adjustments to the water level. 

A second more stable method of humidification was devised. In this method, 
the flow rate of gas flowing to the bubbler is carefully controlled so that 
precisely the desired amount of humidification will be obtained at the bubbler. 
The condenser is isolated from the coolant system. In this method, the means 
of humidity control is independent of the bed coolant system and does not 
drift with the adsorption rate. The bubbler water level is also easier to 
control. For low levels of the humidification, for example, humidifying the 
laboratory dry nitrogen to -40°F dew point, the apparatus can operate unattended 
for a period of several days without the dew point drifting more than 1° or 2°F. 
Whenever humidification to a higher level (usually 50°F dew point) is desired, 
some supervision and control of the bubbler is required; that is, small ad- 
justments might be required as often as every 30 minutes. 


Instrumentation 

Pressures are measured at all critical points in the flow system with 
precision Wallace and Tiernan gages. For unattended operation, critical 
pressures can be measured and recorded by Baraton pressure transducers. 
Temperatures are measured throughout the system and in the beds by copper 
constantan thermocouples. All temperatures are continuously monitored and 
recorded during test periods. As noted earlier, all flow rates are measured 
by rotameters calibrated at 20 psia. 

Water concentrations are measured by means of dew point hygrometers.* 
Sample lines to the hygrometer system are provided at the bed inlet, at the 


* In such instruments, dew points are measured above 32°F and frost points 
below 32°F. For simplicity, the term dew point is used throughout this 
report. Below 32°F proper conversion to water vapor pressures was made 
by use of tables based on water vapor equilibrium with ice. 
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outlet, and in the bypasss system. For those runs dealing with dew points above 
room temperature, the entire sample line system can be heated with electric 
strip heaters. During the first series of tests, dew point measurements were 
taken with a Cambridge Model 992 automatic hygrometer and with an AiResearch 
manual hygrometer. The Cambridge hygrometer, at that time, had the capability 
to measure and continuously record dew points from about +10° to I20°F. The 
manual AiResearch instrument can be used from room temperature down to - 1 00 0 F ; 
but as a manual instrument, it could produce only i nstantaneous readings. 
Generally such readings were taken about every 15 to 20 minutes. 

Later in the program, the Cambridge hygrometer was modified to allow it 
to read over an extended range, from -100° to +1 20°F. During the early part 
of adsorption breakthroughs, the Cambridge hygrometer was set to read the 
inlet dew point to the test bed. Once it was confirmed that conditions had 
settled down, the instrument was switched to the outlet stream. In this manner 
almost the entire breakthrough curve was recorded by this instrument. During 
the rest of the breakthrough, the inlet dew point was monitored with the manual 
AiResearch hygrometer. Later on in the test series, a second Cambridge hygro- 
meter was purchased and installed in the apparatus. In the final conf i gurat ion, 
inlet and outlet humidities are continuously monitored and recorded. Finer 
control of the adsorption conditions is afforded by this setup. Also, with no 
need for manual readings, the test procedure is greatly simplified. 

The hygrometers were subjected to a rigorous calibration procedure. It 
was found that the accuracy of the Cambridge instruments was good, usually 
within ±1 °F, and no corrections were employed,. However, the AiResearch 
hygrometers were found to read considerably low in the low dew-point range. 

Where possible, corrections were made to these readings. Because the region 
where these corrections were made represents very low water concentrations, 
the shape of breakthrough curves is not greatly affected by the inaccuracy 
of the manual unit. 

In early experiments, CO 2 concentrations were measured by a Beckman IR 
I 5A infrared analyzer. Later a wider range, more stable, instrument, the 
Beckman IR 3I5B, was procured. This infrared analyzer was capable of much 
greater accuracy in lower CO 2 concentration regions, and is thus more valuable 
for those experiments dealing with CO 2 partial pressures below 7 mm Hg. The 
infrared analyzers were calibrated by means of pre-mi xed and analyzed gases; 
five different gas mixtures ranging from 0.5 to 3.0 volume percent CO 2 in 
nitrogen were available. At all times when taking readings with the infrared 
CO 2 analyzers, the pressure at the instrument's detecting cell was held 
precisely at 7.0 psia. This pressure was read on a precision Wallace and 
Tiernan gage. To facilitate the sampling of the various test and calibration 
gases by the one CO 2 analyzer, a special sample manifold was installed in 
the apparatus. The flow rate of all sample flows was measured by rotameters. 

For those experiments in which appreciable quantities of both water vapor 
and CO 2 were mixed into the nitrogen stream, water traps containing magnesium 
perchlorate (trade name Anhydrone) were included in the sample lines leading 
to the C0p analyzer. Calibrations verified that these dryers produced accurate 
CO^ readings. 
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Test Procedure 


Each test in this program was preceded by a prolonged vacuum bakeout 
of the sorbent bed. Usually 450° to 500°F was attained with? n the bed at a 
pressure of approximately 10 to 20 p, Hg. The objective of the vacuum bake- 
out was to reduce the water content of the bed to about 0. I percent or less. 

The actual adsorption run was preceded by a period of usually about 
2 hr in which the proper test conditions were established. These included 
bed temperature, bed pressure, flow rate, inlet dew point, and inlet COo 
partial pressure. During this period, the bed was closed off and the fTow- 
stream was directed around the bed through the bypass system. When all con- 
ditions were satisfactory, the timer clock was initiated, and the flow was 
directed through the bed. During the first minute or two of the adsorption 
run, flow rates, pressures, temperatures, and inlet concentrations were given 
final adjustments to compensate for the slightly different conditions of the 
bypass system and the bed system. Throughout the adsorption breakthrough, 
inlet and bed conditions were continually monitored and, when necessary, small 
adjustments were made. At regulator intervals, all items of data were re- 
corded on data sheets. This was in addition to the temperatures, dew points, 
and COg concentrations which were continually recorded. 

Norma 1 1 y, breakthroughs were carried out in one continuous run, until 
inlet and outlet conditions were verified to be the same. For some of the 
runs dealing with the adsorptions of water on molecular sieve I3X or silica 
gel, it was not possible to man the equipment during the full period of the 
breakthrough. In these instances, with the I /2-in. -fin bed, it was decided to 
suspend the test, and shut off the bed and cool it below its adsorption temper- 
ature. The adsorption run was then restarted on the next day. The results 
of these long tests seem to be quite reasonable, and very little effect was 
evident from the interruption in the test procedure. With the 1/2- in. -fin 
bed, the temperature profiles were usually fairly uniform; that is, temper- 
ature excursions at the adsorbing front were not large. Thus, the restart 
of the bed at nearly uniform temperature was not a large departure from the 
conditions at the time of shut down. The continual cooling of the bed during 
the idle period tended to maintain the sorbate loading distribution on the 
bed. On the other hand, the I -in. -fin bed tests were characterized by large 
temperature excursions at the adsorption front. For these tests, it was 
undesirable to interrupt the adsorption procedure. Therefore, as much as 
possible, continual coverage of the adsorption runs was afforded. Even then, 
there were some condi ti ons where test suspension and restart were necessary. 
These runs did show appreciable effects of the interruption. 

For CO 2 adsorption runs, the infrared analyzers were calibrated at 
regular intervals by means of the special bottled gases. In some instances, 
several calibrations were made during a run. For shorter runs, the calibra- 
tions of the infrared analyzer were done at the beginning and at the end 
of the test. It was usually found that the calibration drift was very slight 
during the runs. This was especially true for the newer infrared analyzer. 
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One goal of the test program was to obtain information on gas-stripping 
of water from silica gel and I3X molecular-sieve beds. Thus, after some of 
the water breakthroughs were completed, a gas-stripping test was conducted. 
These tests usually involved just the turning off and isolating of both the 
bubbler and the humidity-control condenser, and allowing dry nitrogen to pass 
through the bed. Other than dispensing with the humidity-control functions, 
the test procedures were the same as during absorption breakthroughs.? 

"i 

In some instances, however, it was desired to run the gas-stripping 
operation at a different temperature than the adsorption run. As discussed 
earlier, one goal of the dynamic tests was to produce data on mass-transfer 
parameters while eliminating heat-transfer effects from the evaluations. 
Accordingly, for these particular gas-strip operations, a transition desorp- 
tion sequence was introduced, altering bed temperature to the conditions at 
which the gas strip was to be run. Usual conditions for this type of test 
were a 60°F adsorption to be followed by a 90°F gas strip* To accomplish 
this, after complete breakthrough at 60°F had been ascertained, the bed tem- 
perature was brought up to 90°F without changing any of the other test param- 
eters. Normal data monitoring and recording continued. During the transition 
period, because of the increase in bed temperature, the outlet humidity rose 
above that of the inlet and stayed above it for a period of time. Finally, 
enough water was removed from the bed and, again, the inlet and outlet condi- 
tions became equal. Then the bubbler and humidity control condensers were 
turned off and isolated, and dry nitrogen was allowed to pass through the bed. 

The data from the transition period was of no real use in the evaluation 
of mass-transfer parameters because of the changing bed temperature. However, 
the recording of inlet and outlet concentrations allowed calculation of the 
amount of water removed from the bed. Similarly, the calculated amounts of 
water left on or removed from the bed during the adsorption period and the 
following gas strip allows a complete material balance to be obtained. This 
material balance compares total sorbate influx to total efflux. It does not 
relate to equilibrium conditions, and is thus independent of bed temperature 
measurements. In effect, it is a direct check on the concentration measure- 
ments. 

In some rare cases, the laboratory nitrogen contained much more water 
vapor than usual (normally -80°F dew point or less). To properly run gas 
strips and other low humidity tests when this occurred, an auxiliary drying 
bed was installed in the system. Molecular sieve I3X was the sorbent used 
in this bed. 
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CHARACTERISTICS OF BREAKTHROUGH CURVES 


Data Reduction and Presentation 

The raw data taken from the mass-transfer apparatus was subjected to a 
data reduction procedure to smooth out experimental upsets and to provide 
information in more meaningful terms. The data reduction was initiated by 
checks of instrument calibrations. Then, On the strip charts of the hygro- 
meter and CO 2 analyzer outputs, short-term discrepancies due to pressure and. 
temperature fluctuations, adjustments in flow rate, and adjustments in inlet 
conditions were smoothed, (in some instances, where the fluctuations were 
large, the data was judged to be of such poor quality that the run was 
rejected and set aside for a rerun. ) This was followed by the transformation 
of dew-point and COg-vol ume-pereent readings ipto concentration readings. 

From publ i shed tables, the dew-point readings were transformed into water vapor 
pressures; water vapor concentrations were then obtained 



This concentration unit has the units of mass of water per mass of nitrogen. 

From the calibration curve produced at the time of the test, the CO 2 
analyzer output trace was transformed into CO 2 concentration values (in the 
units of volume-percent CO 2 ). The$e values were then used to obtain mass 
concentration units 



(Vol f 0 C0 2 ) 



44 

28 


(4-2) 


Since water traps were included in the sample lines to the CO 2 analyzer, 
this formula was valid regardless of whether water was involved in the test. 
The units of this concentration are mass of CO 2 per mass of nitrogen. 


The inlet and outlet concentrations were plotted against test time to 
produce the desired breakthrough curve. The unit of concentration (mass of 
sorbate material per unit mass of nitrogen) selected for the presentation 
of breakthrough curves may not allow one to easily obtain a feel for the 
mass-transfer operation by mere inspection of the breakthrough curves. That 
is, one might have a greater inherent feeling for the adsorption process if 
partial pressures of water or COp were plotted versus time. However, the 
mass concentration units are much more useful when detailed design and per- 
formance predictions are considered. Fi rst of all, by integrating either 
graphically or numerically between the inlet and the outlet traces, one is 
able to obtain the amount of sorbate left on the bed at any particular time. 
This calculation can then be used to obtain removal -efficiency vs time. In 
system designs, this information is used to set cycle times. 
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When the integration is carried out for the entire breakthrough/ the 
result is the ultimate sorption capacity of the sorbent bed. Dividing by the 
dry sorbent mass yields, by definition, the equilibrium loading of the sorbent 
with respect to bed temperature and inlet sorbate partial pressure. A direct 
comparison with published equilibrium data, or that data taken from a static 
test, is then possible. 1 

As a goal, it was desired that the experimentally deHVed bed loading 
would be within ±5 percent of that obtained from equilibrium data. In those 
cases where the comparison (usually called a material balance check) was 
considerably larger than ±5 percent, there was usually something wrong with 
the experimental run in question. Thus, the derived loading and the resulting 
material balance check were used as indicators of the quality of the adsorption 
runs. ‘ ' 

There are several 'different ways in which the material balance check was 
used to improve, or at 1 east to certi fy, the; quality of the data taken in the 
program, For example, in the f i rst runs conducted in the program concerning 
water adsorption on silica gel, it was consistently found that the bed adsorbed 
considerably more water than predicted by equilibrium. Detailed checks of the 
experimental equipment and instrumentation revealed no reason for this devia- 
tion. Data reduction procedures were carefully checked without finding any 
errors. Even selected test reruns did not : shed light on the problem. Finally, 
suspicion was cast upon the silica gel equilibria as published by the manufac- 
turer. Up to that time no silica gel equilibrium data had been taken in this 
program; therefore, it was decided to obtain 60° and i 20 0 F of isotherms. 

(These have already been discussed and presented in Section 3.) The effect 
was to show that the silica' gel as presently manufactured by the Davison 
Company is much better than would be expected from published information. 

In fact, in a personal communication, it was related that the manufacturing 
process had been altered in the hope of improving the sorbent capacity. 

However, no new equilibrium capacity data had been available. 

} 

Whenever a material balance check was decidedly negative • (that is, the 
derived bed loading was much less than the predicted equilibrium loading), 
either a leak into the bed system from the ambient, or a pOor vacuum bakeout 
was suspected. Invariably, negative material balances of as much as 20 to 
30 percent were narrowed down to one of these two causes. 

Large positive discrepancies in the material balance checks (that is, 
the bed picking up considerably more than equi 1 i bri urn), 'were invariably con- 
nected to an undesired and unknown bypass of the gas stream around the adsorp- 
tion bed. In one situation concerning several runs, upon checking for just 
this sort of a bypass around the bed, a failed (leaking) valve was discovered. 
Other cases involving single runs usually involved the inadvertent opening 
of a valve; in one instance, the failure of a differential pressure gage was 
suspected. 

In some cases, there were material balance checks of about 10 percent 
which could not be attributed to poor bakeouts, to leaks into the system, 
etc. It is suspected that the reason these material balance checks were not 
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better is a shift in the output of the temperature recorder. That is, the 
actual bed temperature was somewhat higher or lower than that desired and 
indicated. This situation is discussed later with the presentation of the 
pertinent breakthrough curves. 

There was always some uncertainty as to what the effective bed tempera- 
ture actually was during a run. There were seven thermocouples within the 
packed bed, one thermocouple each on the inlet and outlet process gas streams, 
and one thermocouple each on the inlet and outlet water/glycol coolant streams. 
During adsorption conditions, there was always a variation in the readings. 

It was considered that the inlet and. outlet temperatures of the gas stream 
or the water/glycol stream were not representative of bed temperature. This 
left the bed interior thermocouples for the determination. One of these 
seven thermocouples was attached to a copper fin, and therefore would not be 
representative of the sorbent temperature. The other thermocouples were 
distributed along the bed length. As the mass-transfer zone proceded along 
the bed, the thermocouples within this zone would read higher than the others. 
The difference could be 3° to I0°F, sometimes more. It was decided that 
effective bed temperature was most properly determined by the temperature 
readings outside of the mass-transfer zone. 

There ; were still variations in these readings, due to transverse tempera- 
ture gradients in the bed packing between fins and uncertainty in thermocouple 
location. The thermocouples were attached to a thin wire which was stretched 
longitudinally down the bed length. Even though thermocouple position could 
be determined accurately in the empty bed, there; was no guarantee that during 
packing, the thermocouples had not been displaced laterally. Moreover, due 
to the offset nature of the fins, and the normal waviness of the thin copper 
fins, not all thermocouples could be placed initially in the center of a 
passage. , Therefore those thermocouples located closer to fin surfaces would 
be cooler than mean bed temperature. The best that could be done vyas to 
average those readings for thermocouples not in the mass-transfer zone. 
Temperature recorder calibration shifts, evidenced in a few instances, 
also caused some uncertainty in bed temperature readings. 

In the foregoing discussions concerning material balance checks, it has 
been assumed that the equilibrium loading of CO 2 or water is not affected 
by the presence of the nitrogen carrier gas. It is noted that the CO 2 ,ancj 
water equilibrium data as presented in Section 3 is taken with only the sorbate 
in contact with molecular sieve. That is, the total pressure and the sorbate 
partial pressure are synonomous. In the dynamic tests, the calculated sorbate 
loadings are those developed in the presence of considerable nitrogen gas. 

From the wel 1 -control 1 ed s i ngl e-sorbate dynamic tests, it is concluded that 
the above assumption is correct; there is negligible effect of the carrier 
gas. • 



Characteristics of Dimensional and Dimensionless Breakthrough Curves 

Figure 4-6 shows a number of possible breakthrough curves. Curve A is 
a hypothetically perfect breakthrough curve. It represents conditions whereby 
all sorbent material that enters the bed is held on i t up until the time that 
the equilibrium loading capacity has been reached. Then the breakthrough is 
instantaneous, and the inlet and outlet streams become equal. Curve B repre- 
sents a realistically attainable, good breakthrough curve. With this type of 
curve, very little sorbate passes through the bed until quite late in the 
breakthrough period. Then the breakthrough is reasonably steep and the outlet 
concentration reaches that of the inlet fairly quickly. Curve B Is represen- 
tative of almost all of the si ngl e-sorbate breakthrough curves obtained under 
this program. Curve C represents a much poorer breakthrough curve. It does 
show a period of time where very little sorbate passes through the bed; then 
an initial period of steep breakthrough is followed by a rather slow approach 
of the outlet concentration to that of the inlet. Curve D represents condi- 
tions of adsorption so poor that the breakthrough of sorbate is immediate. 
Subsequently, a very long time is required for the outlet concentration to 
reach that of the inlet. 

Curve B represents conditions in which the mass-transfer resistances 
in the convective gas stream outside of the pellets, and in the diffusive flow 
of sorbate within the pellets pores, are low. Curve C represents a situation 
in which there is considerable resistance within the sorbent particle, or the 
equilibrium curve shape is somewhat unfavorable to adsorption. That is, high 



Figure 4-6. Typical Breakthrough Curves 
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loadings are only attained at higher partial pressures, and the equilibrium 
curve is pronouncedly convex upward. Curve D represents a condition of high 
convective mass-transfer resistance, possibly coupled with high internal 
diffusive resistance, or conditions where the equilibrium is definitely not 
favorable. 


As discussed previously, the breakthrough curve can furnish data on the 
equilibrium sorption capacity, and the shape of the curve tells something of 
the mass-transfer process. However, fpom the dimensional breakthrough curves 
it is often difficult to make critical comparisons of one sorbent with 
another, or of one sorption condition with another. One of the reasons for 
th-i s i s that the breakthrough curves will lie in differing time positions on 
the dimensional breakthrough graph. Without some sort of time scaling, the 
cu /e shapes will be distorted and comparisons will be difficult. To facili- 
tate critical comparisons of sorbents and sorbent conditions, it is useful 
to express the dimensional breakthrough curves in terms of dimensionless 
quantities. From the concept of the perfect sorbent bed, a representative 
dimensional tijne parameter may be designated. This quantity is the time 
required to introduce just enough sorbent into the bed to saturate it. For 
the hypothetical perfect bed, this is the time at which breakthrough occurs. 
This time, 9, is expressed mathematically as (written for C0 2 as the sorbate) 


9 = 


bed Capacity 


w M_ , 
“ Bed 


Rate of Sorbent Introduction 


W N 2 C C0 2 , inlet 


(hr) 


(4-3) 


where w ro = the bed loading at infinite time; or in other terms, the 
equilibrium loading capacity with respect to the bed 
temperature and inlet partial pressure (mass of CC^/ 
mass of dry sorbent) 

M n . = bed mass (mass of dry sorbent) 

Bed 


W^| = flow rate of the carrier gas (mass of N^/hr) 

C = inlet concentration of C0 o (mass of C0 o /mass of N 0 ) 

CO 2 , inlet * 22 

If one wants. to consider that the bed may have had a preload of sorbate at 
time zero, the definition can be restructured (shifted in concentration) 

( w - w. „ ) M D , 
a ■ , V 00 t=0 / Bed 

w N 2 I u 0 2 , inlet C0 2 , t=0 

where w t=0 = t * ie P re ^ oa ^ at time t=0 (mass of C0 2 /mass of 

dry sorbent) 


(4-4) 
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c h 2 , t=o 


the equilibrium concentration over the bed at time t=0 
and corresponding to the loading w t _Q (mass of CO^/ 


mass of N ? ) 

p 4 

tota 1 


can be converted to 



t-0 


by knowing 


Gas stripping may also be considered in the same manner as adsorption. 
The perfect-bed breakthrough time, 0, is given by Equation (.4-4) in the 
rearranged form • 

k n - w ^ j 
e -i-M > .,/ Ji t 

W N 2 ( C C0 2 , t=0 " C C0 2 , inlet) 


The dimensionless time is simply obtained by dividing the actual test 
jtime by parameter 0 


t W N 2 ( C C0 2 , inlet ■ ^COg, t=oj 

t = © = ' m d . u - v j L t 

; ; Bed 00 t=0 


(4-5) 


A dimensionless concentration ^ is chosen in a simple, classical 
manner 

r - r-iv 

co 2 ,t U CO ,t=0 

T = 7— ~~a*~ (4-6) 

L C0 2 , inlet L C0 2 ,t=0 


In the discussions that follow, it was found that the dimensionless breakthrough 
curves provided the best means of visual comparison of sorption conditions and 
the performance of beds. It should be stressed again, that the primary 
purpose of the s i ngl e-sorbate breakthrough tests was for the evaluations of 
mass-transfer coefficients and intraparticle di ffusi vi ti,es. Rigoro,us compari- 
sons of sorbents and sorption conditions can then be made using the performance- 
prediction programs applied to actual bed designs and systems of interest. 


PRESENTATION OF BREAKTHROUGH CURVES 


Range of Parameters Covered 

Table 4-2 presented earlier has given the range of test conditions 
covered in this program. Table 4-4 relates these test conditions to the bed 
size and Volume, and further relates the test condition? to those of use in 
a spacecraft system. For the breakthrough curves presented in the balance 
of this section, a run number designation has been devised as follows: 
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TABLE 4-4 


FLOW PARAMETERS RELATED TO TEST CONDITIONS 

Based upon 1.25 1b of Linde 5A in the I /2-in. -F? n Bed. 

P T = 363 mm Hg, p CQ = 7 mm Hg, T^ = 70°F 



Flow Rate (ib/hr) 

_ . 

Parameter Related to Bed Size and Flow Rate 

■non. 

2.8 

5*6 

Superficial Velocity (ft/sec) 

0.22 

0.45 

0.89 

Residence Time* (sec) 

2.3 

1 . 1 

0.56 

Space Velocity (Bed Volumes/sec) 

0.45 

0.89 

1 .8 

Reynolds No. 

3.3 

6.7 

13.4 

CO2 Removal Rate, at 100 percent CO2 
Removal Efficiency (lb C02/hr) 

'0.043 

0.087 

0.173 

"Perfect Bed" Breakthrough Time, 0, at 
100 percent CO2 Removal Efficiency (hr) 

1 .94 

0.97 

0.48 

Man Rating at 100 percent CO2 Removal 
Efficiency (No. of Men) 

0.47 

0.94 

1.89 


Individual run designation 

Sorbent/Sorbates (Fin Size) - Xs,p 

i I 

Sorbent abbrevi ati ons ' are 

SG - Silica gel., Davison (Grade 05), 6 to 16 mesh 

5A - Linde 5A molecular sieve, i / 1 6— in- pellets.' 

5AD - Davison 5A molecular sieve, 8 to 1 2 mesh spheres 

: 1 3X - Linde : I3X molecular sieve, l/16-in. pellets 

Sorbate ' desigiiati ons were W for water and CO2 for carbon dioxide. The fin 
sizes used i'fv the tests were l/2-in. and l-in. The X in the run designation 
is the serial number of test run. The s or the p which follows the serial 
number denotes the test was either a gas strip, or a water preloading sequence. 
An example of the run designation is SG/W ( I /2 ) - I . This represents run number 
one with water adsorbing on silica gel in the 1/2— in.— fin bed. 
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Adsorption of Water by Silica Gel 

Figures 4-7* through 4-16 are breakthrough curves for the adsorption 
of water by and the gas stripping of water from silica gel contained in the 
I /2-in. -fin bed. Figures 4-17 and 4-18 are water adsorption breakthroughs 
for silica gel packed in the I — i n. — fin bed. 

For the l/2-in.-fin bed, the calculated sorbent loadings compare 
acceptably with equ i 1 i br i urn da ta — genera 1 1 y within 10 percent. Originally, 
when the data was- first reduced, the tests showed 20 to 30 percent more 
sorption capaci ty than the equilibrium data published by the manufacturers. 

As discussed earlier, new data was then taken which yielded much higher 
capacities. This brought the dynamic and equilibrium capacities into general 
agreement. The remaining deviations are considered to result from errors in 
temperature measurement and control — possibly as much as 3°F. Later in the 
program, refinements were made in the bed temperature instrumentation. 

The basic dew-point and concentration measurements are considered to be 
of good quality. This is evidenced by the closeness of the material balances 
between gas stripping runs and the preceding adsorption runs: -2.5 percent 
between Runs 3 and 4s; +6.4 percent for Runs 6 and 7s; +0.8 percent for Runs 
8 (plus the transition adsorption) and 9s. These balances are not dependent 
on equi 1 i bri urn data or bed temperature measurements; they relate to only the 
concentration (dew point) and flow rate measurements. They are even independent 
of flow rate if it is held truly constant. 

The shapes of the breakthrough curves for the I /2-in. -fin bed vary from 
fair to good. The best breakthrough curve shape is given by Run SG/W ( I /2 ) - I 
(Figure 4-7) for the low flow rate of 1.37 lb/hr at 60°F. Silica gel performs 
quite satisfactorily under the conditions of this run, which are the most 
representative of those that the sorbent would undergo in a spacecraft system. 
There is a fairly long period of time where removal efficiency is over 99 
percent--5 to 7 hr. For higher flow rates at 60°F, the breakthrough begins 
disproportionately earlier. This is shown for 2.87 lb/hr in Run SG/W(l/2)-5 
Figure 4-11), where 99 percent removal efficiency is attained for only I to 2 
hr.#* For 5.37 lb/hr, Run SG/W ( I /2) -8 (Figure 4-14), the breakthrough is 
almost immediate. 

All of the adsorptions at I 20 0 F evidence early breakthroughs, which would 
not be acceptable for predrier applications. However, there was never any 
intent to have predrier beds function i isothermal ly at 1 20 0 F ; the purpose of 
these runs was to obtain information on the temperature dependence of Kg and 


* Due to the large number of pages required for Figures 4-7 through 4-52, 
these figures are presented at the end of this section. 

■iBt-Due to equipment malfunction, inlet dewpoint was nearly uncontrollable during 
this run. The earlier portion of breakthrough should not have been affected 
and this portion of the data is considered satisfactory. However, the entire 
breakthrough was not used in the subsequent evaluation of mass-transfer 
parameters. 
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This is desired as gas strips could be conducted in a spacecraft system 



at I20°F; also, in adiabatic beds, 1 oca 1; temperatures in the adsorption zone 
could approach 100° to I20°F. 

Two gas stripping runs at 90°F are included; Runs SG/W(l/2)-2s and 
SG/W(l/2)-9s. These runs followed adsorptions at 60°F, which, in turn, were 
followed by transition adsorption periods whereby bed temperature was raised 
from 60° to 90°F. For illustration. Figure 4-15 shows the transition adsorp- 
tion following Run SG/W(l/2)-8. The I 20°F gas strips (SG/W( I -2)-4s and 7s) 
were all conducted immediately after the preceding adsorption breakthrough 
was completed. 

From the dimensional breakthrough curves, some tentative conclusions can 
be drawn. First, at low enough flow rates, good performance with silica gel 
can be attained at 60°F. At higher flows, breakthroughs are not as good-- 
possibly this is an effect of channeling within the bed becoming more impor- 
tant. Second, silica gel desorption by gas stripping at 90° and I20°F proceeds 
reasonably well, but not in the same time period as adsorption. With -80°F 
dew-point nitrogen, the desorption is essentially quantitative. 

With respect to use as a predrier desiccant, it should be noted that 
breakthrough curves like that of Figure 4-7 represent the best possible 
performance. This is because the adsorption followed a vacuum bakeout, which 
left the sorbent completely dry. In actual use in a spacecraft, after a large 
number of cycles, a repeating steady-state water loading would be retained 
by the bed. In the face of this resident loading, the bed would have to 
achieve 99.5 percent removal efficiency. Certainly the bed of Run SG/W ( I /2 ) — I 
would not operate nearly as long (5 hr) under cyclic steady-state conditions. 
Detailed computer computations would be needed to obtain the exact allowable 
cycle time. 

The adsorptions carried out in the I -in. -fin beds are not as good as those 
for the I /2-in. -fin bed. This is due to the poorer heat-transfer performance 
of the I -in. fins. These tests and comparisons of fins are discussed in more 
detail later in this section. 

Figures 4-19 and 4-20 present the dimensionless breakthrough curves for 
the l/2-in.-fin bed s i 1 i ca-ge 1 /water runs. Figure 4-19 contains breakthrough 
curves for adsorption runs. Here, the effect of flow rate and sorbent temper- 
ature on the quality of the breakthrough is demonstrated. At 60°F, the 1.37 
Ib/hr breakthrough, is far better than that for 5.37 Ib/hr. Although not shown, 
the early portion of the curve for 2.87 Ib/hr (from Run SG/W ( I /2 ) -5 ; Figure 
4-11 is very close to that for the 5.37 lb/hr, showing a quite early break- 
through. The data of Curve A, 1.37 lb/hr, is well substantiated as two 
additional runs were taken at this condition, each showing nearly the same 
behavi or. 

All of the I20°F breakthroughs show performance similar to s Curve C of 
Figure 4-19. That i s ; much earlier breakthrough than the low flow run at 60°F. 
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The dimensionless representations of the gas stripping runs are contained 
in Figure 4-20. These curves all show that the stripping desorption is carried 
out with some difficulty. There seems to be little flow-rate dependence. To 
emphasize this, the high flow run (5.37 lb/hr) at I 20°F was not plotted, as 
it fell almost precisely on Curve B ( I .37 lb/hr). 


Adsorption of Water by Molecular Sieve I3X 

Adsorption and gas-stripping breakthroughs for water vapor on Linde 
molecular sieve I3X are shown in Figures 4-21 through 4-29. Four adsorptions 
and two gas strips are presented for the I /2-in. -fin bed; and two adsorptions 
are presented for the I -in. -fin bed. For comparison, the four adsorptions 
for the I /2-in. -fin bed are shown on one graph, Figure 4-27. 

The first general observation Is that all adsorption breakthroughs are 
reasonably well shaped, regardless of flow-rate, temperature, or fin size. 

For example. Run l3X/W(l/2)-5 for 2.36 lb/hr at I20°F shows an excellent 
breakthrough (Figure 4-25). This can be contrasted to that of silica gel 
for the same conditions, Figure 4-12. From these tests it appears that 
molecular sieve I3X has much more latitude in adsorption than silica gel. 

This is substantiated by comparing equilibrium data for the two sorbents 
(as contained in Section 3). It is noticed that the sorption capacity of I3X 
is greater for low vapor pressures and does not vary nearly as much with 
temperature. Thus, from the basic equilibrium data, one would conclude the 
same as these tests show, that I3X has more favorable adsorption characteris- 
tics than silica gel . 

On the other hand, the I3X equilibria is not as favorable for gas strip- 
ping. This is substantiated by the gas strips presented here. For the tests 
of Figures 4-24 and 4-26, only a small fraction of the water originally pre- 
sent is stripped off during the rather lengthy tests. There are two reasons 
for this. First, the sorbent has a substantial equilibrium water capacity 
with respect to the inlet vapor pressure of the stripping gas. For -I00°F 
dew-point nitrogen (water vapor pressure of approximately 0.0007 mm Hg), 
molecular sieve I3X has water capacities of about 5.6 and 3.3 percent at 
90° and I20°F, respectively. That is, a substantial fraction of the original 
water could not be removed at all, even at infinite time. Still, additional 
water could have been removed with more test time. For the 90°F gas strip, 
it was estimated that 0.120 lb of water could still have been removed from 
the bed. For the I20°F gas strip the estimate is 0.115 lb. 

The second reason for the poor gas strips is due to a rate phenomenon. 
The rate is low because of a very low driving force for mass transfer. This 
driving force is the difference between the equilibrium partial pressure and 
the actual free stream partial pressure. It is believed that convective or 
diffusional resistances are not appreciable enough to cause the slow desorp- 
tion rate. If they were, the adsorptions would be similarly affected, which 
of course they are not. 
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The calculated sorbent capacities for the adsorption breakthroughs are 
all very close to those obtained from the equilibrium data presented in 
Section 5. These tests indicate capability of the apparatus to produce 
quality data. The I -in. -f i n bed tests show calculated capacities less than 
equilibrium values. This is due to non-i sotherma 1 sorbent condi tions, as a 
result of the poor heat transfer of the larger fins. This situation is 
discussed in detail later in this section. 

Figure 4-30 presents typical dimensionless breakthrough curves for the 
adsorptions. The curve shapes are all quite similar and there is not clear 
cut dependence of temperature or flow rate. It should be noted that Run 
l3X/W(l/2)-5, plotted as curve C, was done on a different packing of the 
bed than the other two runs shown on this figure. Also, Run I3X/W( I /2)-3, for 
2.26 Ib/hr and 60°F, was not plotted as its breakthrough curve was very close 
to curves A and B. 

Figure 4-31 shows dimensionless breakthrough curves for the gas strips. 

It is again evident from these curves that such gas strips are carried out 
with difficulty. There is little difference in the curve shapes for the two 
temperatures. 


Adsorption of CC^ by Molecular Sieve 5A 

Adsorption breakthroughs were made for CC >2 on both Linde 5A molecular 
sieve ( I / 1 6— in. pellets) and Davison 5A molecular sieve (8 to 12 mesh 
spherical beads). A good deal of experience had been obtained in previous 
programs with the pelleted Linde 5A sorbent. Some of this experience was 
in l/4-in.-fin beds, which were suspected of excessive channeling. That is, 
the pelleted material did not seem to pack well into the |/4-in. fins. To 
gain some experience with Davidson 5A molecular sieve to evaluate channeling 
effects with the beaded material, the test program included both 5A sorbents. 

Figures 4-32 through 4-39 present adsorption breakthroughs for CO 2 on 
molecular sieve 5A. Figure 4-40 is a gas strip of CO 2 from 5A. Gas stripping 
of CO 2 from 5A has little practical application, but one such gas strip was 
made for comparative purposes. 

Figures 4-32 and 4-33 show comparisons of the Linde and Davison materials 
for essentially similar conditions. The first observation from these two 
figures is that the Linde material provides better breakthroughs at both 
flow rates. Even though both materials show good breakthroughs, Davison 5A 
might be rejected on the basis of these tests. This decision is tempered 
however when one considers the somewhat poorer Linde 5A breakthroughs which 
were obtained later in the experimental program with another packing of the 
bed, Figures 4-38 and 4-39. There is no apparent reason for the difference 
in quality between the breakthroughs of the two packings. In fact, a rerun 
of the 2.4 Ib/hr condition (Figure 4-39) was made in an attempt to determine 
a reason for the difference. The rerun was virtually the same. 
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The only explanation at present is that the sorbent used in the two 
packings (from different production batches) were manufactured with different 
quality. Possibly, there are variations in the blending and extruding process 
which result in pellets with more constricted pore structure. There are a 
number of instances in the literature where dynamic performance variations 
are thought to be caused by differences in the manufacturing process. It 
is possible that moisture, temperature, pressure, and other process variables 
could have pronounced effects on the developed pore structure. On the other 
hand, it is not expected that equilibrium values woul d be altered significantly 
by the manufacturing process. The basic molecular-sieve crystal structure 
should be easily reproduced; the equilibrium behavior would then depend upon 
only the proportion and the inertness of the clay binder. It seems these too 
could be easily controlled. Thus, differences in batches would not be detect- 
able from equilibrium tests. 

As a result of these tests, no conclusion can be made with regard to 
which 5A material, Davison or Linde, is preferred. In fact, with either, 
it is possible that dynamic performance differences between batches might 
exist. It might be desirable to run small-scale dynamic tests on samples 
before using a particular sorbent batch. 

The gas strip shown in Figure 4-40 indicates a process carried out with 
some difficulty. This is probably due to the unfavorable equilibrium charac- 
teristics for gas stripping. The stripping of CO 2 is much more complete than 
that observed for water from molecular sieve I 3X. 

The dimensionless breakthrough curves shown in Figure 4-41 show varia- 
tions due to process conditions and due to the different sorbent batches. 

Curves A, B, and D show a slight dependence on flow rate. This variation 
with flow rate is overshadowed by the variation between sorbent batches, 
curves B and E. The variation with bed temperature, curves B and C, is 
slight and may be due to experimental inaccuracies. The curve for Run 
SA/COg (l /2)-4 for 50°F sorbent temperature was not plotted as it was quite 
close to the 70°F curves, A and B. 


Adsorption of CO- by Molecular Sieve I3X 

The primary interest in molecular sieve I3X is as a desiccant or sorbent 
for the predrier beds of a spacecraft C0 2 removal system. In this use, its 
adsorption of CO 2 is negligible due to the poisoning effect of the resident 
water loading. However, there is the possibility that I 3X might be consid- 
ered as a candidate for the CO 2 removal sorbent. Above approximately 2 mm Hg 

p n , it does not have the equilibrium capacity of 5A. For example, at 7 
uu 2 

mm Hg and 70°F, 5A has 33 percent greater capacity. However, below 2 mm Hg, 
I3X is nearly equivalent to 5A. Since spacecraft beds operate predominantly 
in the low loading range, I3X might not compare too unfavorably to 5A. In 
addition, its much larger pore (or window) structure potentially can produce 
better dynamic performance than 5A. 
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Therefore, a series of CC >2 adsorption tests on Linde 1 3X were included 
in the dynamic test series. The breakthrough curves for both the I /2-in. -fin 
and the ! -in. -fin beds are shown in Figures 4-42 through 4-48. Except for the 
high flowrate run with the I -in. -fin bed (Figure 4-48), all of the break- 
throughs are of good shape. As expected from equilibrium data the developed 
loadings are not as large as those obtained with 5A. 

The high flow rate run (5.31 lb/hr) with the I -in. -fin bed, represents 
conditions of bed geometry and sorption conditions which are generally un- 
acceptable. The heat-transfer aspects of the I -in. fin are discussed in detail 
in the following paragraphs. It should be noticed that breakthrough is rela- 
tively soon and requires a relatively long time to complete. 


COMPARISON OF I /2-INCH AND I -INCH FINS 

Of the three beds constructed for the dynamic tests, the bed containing 
the I /2-inch rectangular offset fin was selected as the primary test article, 
to which the other beds would be compared. In earlier development programs 
considerable data had been obtained with l/4-inch fins. Heat-transfer per- 
formance was shown to be good, but packing density was not as high as desired, 
and there seemed to be a tendency to channeling. It was reasoned that the 
1/2- and I -inch fins, especially of the offset variety, would allow less 
channeling, but their heat transfer adequacy was questioned. To obtain at 
least a rough idea of the heat-transfer performance of the larger fins, 
a hand, relaxation finite-difference computation was performed for the 
I -inch fin. Silica gel was chosen as the sorbent; the calculation was made 
for steady-state conditions representative of vacuum desorption. It was 
expected these conditions would produce the largest temperature drop in the 
packed bed. 

Figure 4-49 shows the results of the finite-difference computation. Due 
to symmetry, the calculation was performed, for a quarter of one channel 
formed by the I -in. fin matrix. The maximum temperature difference of 46° F is 
much higher than desired, and woul d not allow efficient vacuum desorption. It 
should be noted that the fin effectiveness is not limiting; that is, the fin 
is nearly isothermal. The limiting factor is the low bed thermal conductivity. 
Under adsorption conditions with a flowing gas, a slightly lower temperature 
drop would be expected. However the limiting factor would still be the low 
bed conductance. 

It was concluded that the I -inch fin would probably have poor heat- 
transfer performance. Therefore, the I /2-inch-fin test bed was used for the 
bulk of the experimental runs. Several runs were conducted with the I -inch-fin 
bed for comparative purposes. The dimensional breakthrough curves for these 
runs have already been presented. The following paragraphs present compari- 
sons of the performance of the two beds, and some observations of their heat- 
transfer performance. 
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Figure 4-50 shows the comparison of dimensionless breakthrough curves 
obtained for water on Linde I3X molecular-sieve pellets with the 1/2- inch-fin 
bed and the I -inch-fin bed. The poorer breakthrough performance of the I -inch- 
fin bed is due strictly to heat-transfer deficiencies of this bed. In both 
cases water/glycol heat-transport fluid was fed to the bed coolant passages 
at 59° to 60°F, at such a rate that the temperature rise of the coolant was 
less than I °F. From the thermocouples placed within the bed it was easy to 
•conclude which bed had the better heat-transfer characteristics. For the 
! /2-inch-fin bed the traveling mass-transfer zone was evidenced by temperature 
readings of about I 8 0 F above the coolant temperature. For the I -inch-fin bed 
the temperature excursions were 35° to 40°F above the coolant temperature. 

The resulting shift with temperature in the sorbent's equilibrium capacity 
causes the breakthrough curve for the I -inch-fin bed to be less steep and more 
prolonged than that for the bed with the closer-spaced heat-transfer matrix. 

Figure 4-51 shows the same general behavior for the adsorption of water on 
silica gel. Here the temperature excursions from 60°F are not as noticeable-- 
4° to 5°F for the 1 /2-inch-fin bed and 8° to I0°F for the I -inch-fin bed. 
However, since silica gel's sorption capacity for water is much more tempera- 
ture dependent than molecular sieve I3X, the early portion of the breakthrough 
curve (which is the important consideration for actual spacecraft systems) 
is much poorer for the I -inch-fin bed. 

One other point can be made to show the poorer performance of the larger 
fin. All runs with the I -inch-fin bed were made with no control on the bed 
inlet gas temperature. Thus, the temperature at the bed inlet was generally 
the same as the laboratory ambient, about 75°F. Because of poor bed heat 
transfer, the first 10 to 15 percent of the bed operated above 60°F. In 
contrast, with no inlet temperature control with the l/2-inch-fin bed, usually 
a negligibly small portion of the bed was above 60°F.# The result of this 
higher temperature region in the entrance of the 1 -inch-fin bed is a lower 
overall bed sorption capacity. Indeed, most water adsorption runs with the 
I -inch-fin bed show material balance checks some 10 to 15 percent below that 
predicted by the 60°F sorbent equilibrium capacity. 

It should be noted that it is quite conceivable that some spacecraft 
systems would not include precoolers for the airstream flowing to a sorption 
bed. In such cases, if isothermal bed operation were desired, the l-inch fin 
would be downrated with respect to cooling the airstream in the entrance 
region of the bed. 

Figure 4-52 for CO 2 adsorption by Linde I 3X shows again the same general 
comparison between the two beds. This was not entirely expected as the heat 
effects of CO 2 adsorption are not as much as that for water vapor. However, 
excursions of temperature in the mass-transfer zone were noted as high as 
1 6 0 F for the I -inch-fin bed. 


* For very high flow rates, for example, 5.37 Ib/hr for Run SG/W(l/2)-8, Figure 
4-14, even in the l/2-inch-fin bed an appreciable entrance region was above 
the 60°F coolant temperature. This led to less sorbent loading than expectfed 
if the bed were truly isothermal. 
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The conclusion can be made that when internal fins are desired, the I -inch 
size fin should not be used. In the test series, it showed only 1.7 pefceHt 
greater packing density than the I /2-i nch-f i n. bed, and the adsorption perfor- 
mance was significantly and repeatably poorer. Based upon these tests, it Is 
considered that 3/8-inch offset fins would also produce very good thermal 
performance. The packed bed density would be a little less than l/Z-inch f i n$j 
but temperature excursions in the mass-transfer zone would be even less* 
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4-7. Water Vapor Adsorption Breakthrough on Silica Gel, 
Run No. SG/W ( I /2 ) — I 
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Gas Strip of Water from Silica Gel, 
Run No. SG/W(l/2)-2s 
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Bed fin surface: 1/2- In. 

Flow rate: 1.37 Ib/br 

Bed temperature: 121 °F 

Inlet p H Q : 9.21 mm Hg (50°F O.P.) 

Total pressure: 264 mm Hg 

Water Influx: 0.122 lb 

Water efflux: 0.050 lb 

Removal off. at bt: 59.0# 

Water left on bed: 0.072 lb 

Bed mass : I .086 lb 

Sorbent loading: 0.066 

Loading fromequtl. data: 0.061 

Deviation fromequtl: +8.2# 


Figure 4-9. Water Vapor Adsorption Breakthrough on 
Silica Gel, Run No. $G/W(l/2>-3 
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3£Bod fin surface: l/2-ln. 

Flow rate: 1.37 Ib/hr 

Bed temperature : I 20®F 

Total pressure: 260 nan Hg 

Calc. Initial sorbent loading: 0.061 

Bed mass: 1.086 lb 

Initial sorbate on bed: 0.072 lb 

Net sorbete removed: 0.070 lb* 

Deviation of mat’l balance: -2.5# 
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Figure 4-10. Gas Strip of Water from Silica Gel, 
Run No. SG/W(l/2)-4s 
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Figure 4-11. Water Vapor Adsorption Breakthrough on 
Silica Gel, Run No. SG/W(l/2)-5 
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Figure 4-12. Water Vapor Adsorption Breakthrough on 
Silica Gel, Run No. SG/W( 1/2) -6 
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Figure 4-15. Water Vapor Adsorption Transition on 
Silica Gel, Run No. $G/W(l/2)-8 
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Bed fin surface: l/2-inch 


ate: 5.37 Ib/hr 

■Bed temperature: 9 , '°F 

Pressure: 258 mm Hg 


Cal<;. Initial sorbent loading: 0.148 

Bed mass: I .086 lb W “ 

Initial sorbet* on bed: 0.161 lb 

■‘Net sorbate removed; 0.164 lb 




•1 J".. 


Time (hrs) 


Figure 4-16. 


Gas Strip of Water From Silica Gel, 
Run No. SG/W(l/2)-9s > 
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Figure 4-17. 
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Water Vapor Adsorption Breakthrough on 
Silica Gel, Run No. SG/W(l)-J 
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Bed fin surface: I -inch 

Flow rate: 5.3 Ib/hr 

Nominal bed temperature: 60°F 

Inlet P u 9.21 mm Hg (50°F D.P.) 
H 2° 

Total pressure: 367 mm Mg 


Water influx: 0.929 lb 

Water efflux: 0.500 lb 

Removal eff. at bt: 46.2?> 

Water left on bed: 0.429 1b 

Bed mass: 1 .350 lb 

Sorbent loading: 0.318 


Loading from equil, data: 0.376 
Deviation from equil; -15.43b 
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Figure 4-18. 


Water Vapor Adsorption Breakthrough 6h 
Silica Gel, Run No. SG/W(l)-2 
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Figure 4- 19 . D imens ionless Breakthrough Curves > Adsorption of Water by Silica Gel 
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Figure 4-20. Dimensionless Breakthrough Curves, Gas Stripping of Water From Silica Gel 
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Figure 4-21. Water Vapor Adsorption Breakthrough on Linde I3X, 
Run No. I3X/W(I 



Figure 4-22. Water Vapor Adsorption Breakthrough on Linde I3X. 
Run No. !3X/W(i/2)-2 
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Figure 4-23. Water Vapor Adsorption Breakthrough on Linde I3X, 
Run No. !3X/w(l/2)-3 
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Figure 4-25. Water Vapor Adsorption Breakthrough on Linde I3X, 
Rum No. !3X/W(l/2)-5 
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Figure 4-26. Gas Strip of Water From Linde I3X, 
Run No. I3X/W(!/2)-6 s 
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Figure 4-27. Water Vapor Adsorption Breakthrough Curves for Linde I3X 








Figure 4-28. Water Vapor Adsorption Breakthrough on Linde I3X, 
Run No. I 3X/W( I ) - I 



Figure 4-29. Water Vapor Adsorption Breakthrough on Linde I3X, 
Run No. I3X/W( I )-2 
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Figure 4-30. Dimensionless Breakthrough Curves, Adsorption of Water by Linde I3X 
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D i mens ion>l ess Breakthrough Curves, Gas Stripping of Water 
From Linde I3X 
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Figure 4-34. 


CO 2 Adsorption Breakthrough on Linde 5A, Run No. 5A/C0 2 ( I /2)-3 
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Figure 4-35. 
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Figure 4-36. CO^ Adsorption Breakthrough on Linde 5A, Run No. 5A/C0 2 ( I /2)-5 
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Figure 4-37. CO 2 Adsorption Breakthrough on Linde 5A, Run No. 5A/C0 2 ( I /2) -6 
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Figure 4-38. C0 2 Adsorption Breakthrough on Linde 5A, Run No. 5A/C0 2 (l/2)-7 
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Figure 4-39. C0 2 Adsorption Breakthrough on Linde 5A, Run No. 5A/C0 2 ( I /2)-8 


1 18 




< 3 ) 

E 


I 19 


Figure 4-40. Gas Strip of C0„ From Linde 5A, Run No. 5A/C0„( I /2)-9s 
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Figure 4-41. Dimensionless Breakthrough Curves, Adsorption of CO by Linde 5A 
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Figure 4-44. CO 2 Adsorption Breakthrough on Linde I3X, 


Run No. !3X/C0 2 (l/2)-3 
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1/2- inch 

^FIow rate: 1.23 ib/hr 

:t Bed temperature: 70°F 

Inlet ^ CQ : 7.08 mm Hg 

-!: Total pressure: 364 mm Hg 

Jnf lux: 0.070 Jb 

C0 2 efflux: 0.016 lb 

Removal eff. at BT: 76^ 

C0 2 left on bed: 0.054 lb _ 

Bed mass: 1.004 lb 

Sorbent loading: 0.0534 Ib/lb - 

Loading from equM. data: 0.05l_^ 

Deviation from equil: +4.8$ 
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Figure 4-45. CO 2 Adsorption Breakthrough on Linde I3X, 
Run No. !3X/C0 2 (l/2)-4 
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Figure 4-47. CO 2 Adsorption Breakthrough on Linde I 3X, 
S Run No. I3X/C0 2 ( I )-l 
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Figure 4-48. C0_ Adsorption Breakthrough on Linde I3X, Run No. I3X/C0 9 (I )- 




Heat transfer 

fluid channel 



120 

120 

120 120 

120 

120 


I 16.2 


III. 6 104.0 


108.7 


114.9 105.7 


III. 3 110.2 109.6 


Silica Gel 
98.1 91.0 


Bed (1/2"). 


Assumptions : 

Desorption rate, 250 Btu/hr-lb 
silica gel 

Bed effective conductivity, 0.1 

Btu/hr-ft 2 -°F/ft 

Aluminum fins, i in. high by 
I in. spacing, 0.006 in. thick 


Figure 4-49. Temperature Distribution (°F) in a Desorbing l-in 
Channel of Silica Gel 
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Water on Molecular Sieve I3X 
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SECTION 5 


EVALUATION OF MASS-TRANSFER COEFFICIENTS AND 
INTRAPARTICLE DIFFUSIVITIES 


DEFINITIONS 

In this section, values of gas-phase mass-transfer coefficient K and 
intraparticle diffusivity D are determined for the single-sorbate breakthrough 
tests presented in Section 4. 

From Reference I, the intraparticle diffusivity for sorbate k is defined 
by the diffusion equation for the interior of a sorbent pellet: 

Sw k D k i » 2 aw k , , 

at - ” iTT 3 " a? r ~ < 5 -'> 

2 

D^ has the units ft /hr; the local sorbent loading is in the units mass of 
sorbate k per mass of sorbent. Time t is in hr; r in ft is the effective 
distance from the spherical pellet center. Sorbent particle density p s is in 
the units of lb-mass/ft 5 . 


The gas-phase mass transfer coefficient Kg is defined by the following 
relation expressing the sorbate flux at the sorbent pellet boundary: 


dw k 

' p s D k M^r 


K 


g 



- P V 


(5-2) 


2 

The mass-transfer coefficient has the units of 1 b-mol es/hr-ft -mm Hg; is 
the molecular weight of sorbate k. The total pressure P is i n mm Hg; p^g is 
the sorbate partial pressure at the sorbent boundary, mm Hg; X|< is the 
molefraction of sorbate k in the free stream. Note P'X^ - p^ the sorbate 
partial pressure in the free stream, mm Hg. 

The mass-transfer coefficient is also used in the gas phase material 
balance in the x- or longitudinal flow direction: 


dp. P M 

j — = r — — a K (p. - p. ) 

dx f p u sg g *ks k 

H g 9 


(5-3) 


Here a S g is the sorbent specific external surface area, sq ft per unit volume 
of packed bed. M q is the average molecular weight of the process gas; f is 
the bed void fraction, pg is the gas density, lb-mass/ft^; Ug is the inter- 
stitial or true gas velocity, ft/hr. 



129 



It should be noted from Equation 5-3, that Kg and a S g must both be known. 
In any evaluation of Kg, the corresponding value of a S g must be given. Table 
5-1 presents several pertinent values, including a S g, used in the evaluations. 

TABLE 5-1 


PARAMETERS USED IN MASS-TRANSFER COEFFICIENT EVALUATIONS 


Parameters 

Silica Gel 
(6 to 16 
Mesh) 

Li nde 5A 
(.1/16-1 n. 
Pel lets) 

Davison 5A 
(8 to 12 
Mesh) 

Linde I3X 
(1/16-in. 
Pel 1 ets) 

Particle density, 
P s , lb/ft3 

75.0 

64.0 

66.4 

64.0 

Packed density, 

p sb , lb/f t 3 

42.96 

43.55 

49.73 

42.97 

41.06 

39.73 

Surface area per 
uni t volume, 

a sg' ft2/ft3 

575.0 

636.0 

653.7 

636.0 


METHOD OF DETERMINATION 

The method used for the determination of the coefficients Kg and D 
involves the use of the performance-prediction program, S9960, discussed in 
Section 7 of this report. If this program is considered to be a valid mathe- 
matical simulation of the actual bed mass-transfer phenomena, and if the true 
values of Kg and D are included in the program input data, the predicted per- 
formance should match that observed. Conversely, if a prediction matches the 
experimental performance, then the values of K g and D used by the program must 
be the true values. Stated in simple terms, after producing a number of per- 
formance predictions for various pairs of values of Kg and D, the adopted 
values were those which produced the best fit of the predicted performance 
to that of the experimental test. The judgment of the best fit was by a visual 
comparison of predicted and actual breakthrough curves. 

In order to aid the visual comparisons, which at times were difficult, 
a systematic computation and graphing procedure was set up. Table 5-2 shows 
the values of mass-transfer coefficient and diffusivity for which performance 
predictions were usually obtained. It should be noted from the table that 
the performance predictions were in groupings of constant diffusivity. For 
each grouping (that is, for a constant value of diffusivity) the various per- 
formance predictions were plotted on one graph wi th mass-transfer coefficient 
as a parameter. For convenience, all predictions were plotted on xerographic 
copies of the experimental breakthrough curve. For each plot, or for each 
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VALUES OF MASS-TRANSFER COEFFICIENT K AND DIFFUSIVITY D 

9 

USED IN EVALUATIONS 


Mass -Transfer Coefficient, 
2 

(ib-mole/hr-ft -mmHg ) 

0.00025 

0.0005 

0.001 

0.005 

0.000 I 
0.00025 
0.0005 
0.001 
0.005 

0.000! 

0.00025 

0.0005 

0.001 

0.005 

0.000 I 

0.00025 

0.0005 

0.001 

0.005 

0.0001 
0.00025 
0.0005 
0.00 1 
0.005 
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grouping, the best value of mass-transfer coefficient was selected and a repre- 
sentative curve drawn. Sometimes this curve was drawn precisely through the 
predicted points. In other cases, the best curve was drawn by visual inter- 
polation between predicted curves, interpolating with respect to K . Finally, . 
for all five values of diffusivity, the best breakthrough curves were compared, 
and a tentative selection of mass-transfer coefficient and diffusivity was 
made. 

For all evaluations presented here, the bed was broken into 20 finite- 
difference elements in the flow direction. For early computations, as many as . 
6 i nternal -parti cl e finite-difference elements, that is, up to six radial 
layers, were considered for calculation of diffusion in the pellet. However, 
it became evident that the best predictions occurred for large values of dif- 
fusivity which produce very small loading gradients within the sorbent particle. 
Therefore, the number of sorbent finite-difference elements was reduced to two. 
This produced distinct savings in computer time. For most of the adsorption or 
gas strips, about 4-1/2 minutes of computer time was required on a Univac 1108 
system to process the Kg/D sets given in Table 5-2. In some cases, computa- 
tions for lower or intermediate values of K were made in addition to those of 
Table 5-2. 9 


DISCUSSION OF EVALUATION RESULTS 


Sample Evaluation 

Figure 5-1 shows predicted points for Run No. 1 3X/C02 ( l/2)-2. The pre- 
dictions for this run. are typical in many aspects to all of the other evalua- 
tions for adsorptions. First, at low di ffusi vi ties (0.000004 and 0.00001 
ft^/hr) no value of K q gives a good prediction. As diffusivity increases, so 
does the accuracy of the prediction. Usually thq best prediction is obtained 
for D = 0.001 f tVhr. With this diffusivity, the computed sorbent loadings 
(part of the printed output of the program) show essentially no radial varia- 
tion; that is, the calculated resistance to diffusional. sorbate transfer within 
the pellet is so small that no concentration or loading gradient is developed. 
The same is essentially true for D = 0.0001 ft^/hr . 

One might conclude at this point that diffusion within the sorbent is 
negligible and all of the mass-transfer resistance is in the convective trans- 
fer from the free stream to the sorbent surface. This will be discussed in 
more detai 1 later. 

A second observation is that Kg does have a definite influence on curve 
shape. The dependence of Kg is not, however, linear. For example, using the 
D = 0.001 ftVhr plot of Fi gure 5-1, there is a marked difference between the 
Kg = 0.0001 curve and the Kg = 0.00025 curve. On the other hand, the differ- 
ence between Kg = 0.0005 ana Kg - 0.001 is not as great. This should be 
remembered when considering the tabled results presented later in this section. 
In effect, there is a range of Kg that will give a good prediction. For 
Figure 5-1, D = 0.001, the best curve was estimated to be for Kg = 0.00037 
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CO, concentration, C (lb CO /lb N ) CO concentration, C (lb C0 2 /lb N^) CO, concentration, C CQ (lb C0 2 /lb N 2 ) 







(midway between 0.00025 and 0.0005) . However, considering experimental 
accuracies in both the dynamic data, and in the equilibrium data which the 
programs use, the range of K_ from 0.00015 to 0.0007, or possibly higher, 
would produce satisfactory predictions. •• 

A third observation typical of almost all cases, is that the value of Kg 
that fits the early breakthrough test does not fit the entire breakthrough 
shape as well. Because of this, two pairs of K_/D values were determined 
for each run: one pair for the best curve shape and another pair for the 
best approximation of the early breakthrough performance. 

A fourth observation was usually made, although the illustrated run does 
not show this well. In several cases, no matter what value of Kg was used, 
the predicted curves were displaced in time from the experimental curve. That 
is, the predicted breakthroughs were earlier or later than the experimental. 

In most cases, the offset in time was in accord with the material balance 
error between calculated (from experiment) and equi 1 i brium sorbent capacities. 
If the calculations for a dynamic test showed greater sorbate pickup than 
equilibrium, the predicted breakthrough curves would be earlier than the 
experimental; if less capacity than equilibrium, predictions would appear 
later. Since it was usually considered that error in the measurement of the 
temperature level of the bed was the cause of the problem, giving more or 
less capacity than expected, but not appreciably altering the real values of 
Kg or D, it was decided that curve shape was the most important evaluation 
criterion for determining Kg.' 

It would be possible to 1 eliminate the offset between predicted curves 
and experimental by appropriately choosing effective bed temperature and using 
the equilibrium isotherm corresponding to it. This was not done, however, as 
it was considered unjustified.- 

The predictions for gas strips showed little variation with diffusivity 
and only slight variation with Kg. The variation with Kg occurred in the very 
early portion of the stripping procedure; later in the gas strip, curves for 
all values of Kg tended to coincide. Most of the predictions tended to show 
more ultimate desorption (higher outlet concentration) than the experiments. 

It is believed this is due to the inaccuracies in the equilibrium data for 
very low partial pressures. • : 


Accuracy of the Procedure 

There are k number of areas where this method of mass-transfer coefficient 
evaluation is vulnerable to inaccuracies in data and technique. Obviously, 
any inaccuracy in the dynamic tests will lead to scatter in resulting values 
of Kg and D. These inaccuracies have been discussed before, bed temperature 
determination being a key parameter. •' 1 1 

Accuracy of the entire equilibrium isotherm is important. The predicted 
performance depends on the entire isotherm, not just the value corresponding 
to bed inlet conditions. In this regard, several runs were subjected to the 
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evaluation procedure with differing equilibrium isotherm shapes, but with 
essenti al ly ' the same total bed capacity. The results showed that predicted 
breakthrough offset and curve shape are significantly changed by isotherm 
shape. The values of Kg and D given here- are based upon the currently avail- 
able equilibrium data, most of which is given in Section 3. 

The numerical procedure used by the program can also introduce inaccuracy 
into the evaluation procedure. Investigations showed that by varying the 
finite-difference size, that is', varying the number of finite-difference ele- 
ments into which the bed is divided, slightly different predicted breakthrough 
curves could be obtained. Furthermore, these different breakthrough curves 
did not always result in the same amount of sorbate material being left on the 
bed. In some instances, the material balance comparison between equilibrium 
data and the capacity obtained from the predicted curve was not as good as the 
comparison made with the experimental, data. As noted earlier 20 longitudi- 
nal finite-difference elements were used for all computations. The magnitude 
of the time step was also important. Small time steps tend to produce more 
accurate predictions; but consume more computer time. For the evaluations, 
the time step in the range of 0.01 hr to 0.025 hr, was chosen to be as large 
as possible while still giving good accuracy. 


Presentation of Results j I 

Tables 5-3 through 5-6 present the results and observations of the mass- 
transfer coefficient evaluations conducted for the singl e-sorbate tests on the 
I /2-inch-fin bed presented, in Section 4. Because of the large temperature 
variations noted in all of the I -inch-fin bed tests, no evaluations were per- 
formed on the dynamic tests conducted with this bed. Run conditions are noted 
in the evaluation tables to help in the assessment of trends, if any. 

It is not difficult to see that the best predictions are almost always 
obtained for a very high value of diffusivity (0.001 or 0.0001 ft /hr). With 
such diffusivity values, there is effectively no resistance to diffusional 
mass transfer in the sorbent particle, and thus no loading gradient is developed 
within the sorbent. It should be noted that the value of D = 0.001 ftVhr was 

not used for silica-gel evaluations. But as can be seen from the other tables, 

the results for D = 0.001 and 0.0001 ft 2 /hr are usually nearly identical. 

In no system does there seem to be a clear relationship of Kg upon flow 
rates. Even though the dimensionless curves seem to indicate weak flow-rate 

dependency, the scatter of the mass-transfer coefficients is large enough to 

mask the effect. . In the I3X/C02 runs there seems to be a slight flow rate 
dependence, but the number of runs and the flow rate variation is not large 
enough to substantiate this. 

For the 5A/C02 runs, the best K g values, are not consistent with flow rate 
either. Ignoring Runs 5A/CC>2( I /2 ) -7 and -8, for the second packing of the bed, 
one might conclude that Run 6, for the very low flow of 0.5 lb/hr, does show a 
lower value of K g which is due to a flow rate effect. However, noting the K g 
values for Runs 7 and 8, as well as the two Davison 5A runs, the flow rate 
dependence is uncertain. 
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EVALUATION OF K AND D FOR MOLECULAR-SIEVE I 3 X/WATER ADSORPTIONS 
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EVALUATION OF K AND D FOR MOLECULAR-SIEVE I3X/C0„ ADSORPTIONS 
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For the water adsorptions on silica gel and molecular sieve I3X, using 
diffusivity D = 0.001 ftVhr, one value of Kg can be chosen for each sorbent 
that will give good performance predictions. These values are 0.0004 lb-mole/ 
hr-ft -mm Hg for silica gel and 0.0002 for molecular sieve I 3X. 

For the COg adsorptions on Linde 5A and I 3X, the scatter of Kg values is 
large enough that it appears from the evaluation tables that no one single 
value can be used to produce acceptable predictions. To help in visualizing 
the results of the mass-transfer coefficient evaluations, Figures 5-2 and 5-3 
were prepared showing the range of K g , for D = 0.001 ft 2 /hr, which give accept 
able predictions for the COg adsorptions. For ISX/CO^ adsorptions Figure 5-3 
indicates that there is a. single value of Kg (= 0.00025 1 b-mol e/hr-f t^-mm Hg) 
which can be used to produce good predictions over the range of experimental 
conditions. However, for the SA/COg adsorptions no one value seems sufficient 

The problem lies in the second set of Linde 5A runs (Runs 7 and 8) and 
the two Davison 5A tests. The sorbent materials used in these runs may have 
had restricted pore networks, hindering internal diffusion. To test this 
assumption, a number of predictions were made for Kg = 0.0004, while varying 
di ffusivi ty over a considerable range. That is, using a mass-transfer coeffi- 
cient that in combination with a high diffusivity seemed to give good results 
for other tests, it was questioned whether a low diffusivity value existed 
which could give good predictions for the runs in questions. It was found 
that no diffusivity value coupled with Kg = 0.0004 gave, really good predic- 
tions; the curve shapes were definitely ai storted. If shape is ignored, 
the combination Kg = 0.0004, D ?= 0.000015 gives a reasonable prediction of 
the breakthrough end points for the Linde* 5A/C02 Run No. 7, 

The possibility of channeling in these tests was considered. If this 
were the cause of the poorer breakthroughs, lower effective Kg values would 
indeed have to be used in the predictions. However, there is no evidence to 
show that channeling caused the poorer breakthroughs in either of these tests 
or the substantiating, nearly identical, reruns made for each system. The 
second packing of Linde 5A material in the I /2-inch-fin bed resulted in a 
larger quantity of sorbent being held in the bed than in the first set of runs 
Thus, it would be expected there would be less channeling, rather than more. 
With regard to the runs with Davison 5A, it was observed that the sorbent 
beads packed very well in the bed. Here too channeling was expected to be 
i nsi gni f i cant. 

It was concluded that for the 5A/C02 tests no one mass-transfer coeffi- 
cient value could produce acceptable predictions over the entire range of 
experimental conditions. Therefore, a range of Kg from 0.0002 to 0.0004 lb- 
mol e/hr-f t 2 -mm Hg must be considered. Possibly, individual tests on sorbent 
samples might be required to pinpoint the value of Kg to be used. Table 5-7 
contains recommended values of mass-transfer coefficient and intraparticle 
diffusivity for all of the systems studied in this program. 
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TABLE 5-7 


RECOMMENDED VALUES OF MASS-TRANSFER COEFFICIENT K g AND 
DIFFUSIVITY D TO BE USED IN PERFORMANCE PREDICTION PROGRAMS 


System 

K 9 

p 

(]b-moles/hr-ft -mm Hg) 

D 

(ft 2 /hr) 

Si 1 ica -Gel /water 
(Davison Grade 05) 

0o0004 

0.0003 

for some conservatism 

0.001 or 0.0001 
0.001 or 0.0001 

Linde l3X/Water 
( 1 / 1 6— i nch pel lets) 

\ 

0.0002 

0.001 

Linde 5A I 

( l/16-inch pel lets) \ 

0.0002 to 0.0004 

0.001 

Davison 5A 1 

(8 to 12 mesh beads)/ 



Linde I3X/C0 2 
( l/16-inch pel lets) 

0.00025 

0.001 


Mass-Transfer Coefficients and Di ffusi vi ti es 
Predicted from Published Correlations 

In order to offer some justification of the values obtained in this study, 
mass-transfer coefficients Kg were calculated from two published correlations: 
one recommended by Vermeulen (Reference 7) and the other due to Lightfoot 
(Reference 5) as used by Lee and Cummings (Reference 6). Using a mean spheri- 
cal particle diameter of 0.008 ft for 1/16— in . Linde molecular sieve pellets, 
and conditions for 2.4 lb/hr CO 2 adsorptions in the l/2-inch-fin bed, both 
correlations gave values in the range 0.0015 to 0.0025 1 b-mol e/hr-f t -mm Hg. 
These values are somewhat higher than the recommended values given in Table 
5-7. However, as can be seen from Figure 5-1, the effect on the breakthrough 
of using the higher Kg values in this range is not as great as might be ex- 
pected. It should be noted, that the recommended values inherently include 
the effects of any channeling in the finned bed, and the non-i sothermal condi- 
tions in the mass-transfer zone. Both of these nonideal conditions tend to be 
detrimental to the breakthrough performance; it is considered that the tempera- 
ture variations in the mass-transfer zone are by far the most important. With 
the program computations, the only way to compensate for these effects on the 
breakthrough curve is to use a lower effective value of K g . 
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From a recent adsorption NTU (number of transfer units) correlation by 
Vermeulen (Reference 7), it appears that the adsorptions of this program are 
characterized as being controlled by a combination of axial dispersion and 
external transport. That is, from the flow Peclet number and the dimension- 
less NTU parameter, data points for the experiments of this program fall on 
the correlation plot generally between the regions of axial dispersion control- 
ling and external transport controlling.* 

This brings up a topic not yet discussed, that of axial dispersion. The 
mathematical model in the performance-prediction programs does not include the 
effects of axial dispersion, i . e. , the di ffusional flow of sorbate from the 
front of the bed toward the rear which is caused by the concentration gradient 
developed across the bed due to the adsorption onto the solid sorbent. For 
those familiar with heat-exchanger design, axial dispersion di ffusional effects, 
can be considered analogous to axial or longitudinal conduction effects on 
high effectiveness (high NTU) exchangers. If one were to test a heat exchanger 
controlled by longitudinal conduction effects, and were to subsequently cal- 
culate from these tests heat-transfer coefficients using simple heat-exchanger 
theory which i gnores these important effects, the effective coefficient values 
would be lower than those predicted by well-known correlations. This would 
also be the case--as evidenced here — with adsorption bed experiments where 
axial dispersion is important. Again, there seems to be justification for the 
Kg values being evaluated lower than published correlations. 

With regard to effective intraparticle di ffusi vi ty for molecular sieves, 
there is physical reasoning to believe that high D values and low loading 
gradients would result from the performance-prediction evaluations. This 
contention relates to the makeup of molecular-sieve pellets, where 3 p, crystal- 
lites (the active sorbent particles) are suspended in a'porous clay matrix. 
Ideally, the matrix will have negligible diffusion resistance; therefore, the 
only resistance would be in the very small crystallites,- Since tHfc program 
computes pore diffusion on the basis of effective spherical sorbent granules 
which have a mean diameter of the same magnitude as the actual pellet dimen- 
sions, the diffusion resistance of the crystal 1 i tes woul d not be handled pro- 
perly. It would then be expected that on the basis of the pellet diameter, 
the intraparticle diffusion resistance would appear to be negligible. There- 
fore, high di ffusi vi ty values would result. 

In actuality, both the clay binder matrix and the molecular sieve crystal- 
lites have some finite resistance to intraparticle diffusion. At the present 
time, little has been done in the adsorption field in the area of modeling 
the diffusion phenomena in sorbents, such as molecular sieves, which possess 
bi -modal pore-size distributions. Previous experimental work conducted at 
Ai Research (Reference 8), has shown an effect of gross-particle size, comparing 
the performance of l/8-inch to I / 1 6- i nch pellets. This work would then indi- 
cate that there is at least some di ffusional resistance in the clay matrix of 


* The NTU is calculated according to the relations given in Reference 7. It 
involves the mid-height slope of the dimensionless breakthrough curve. A 
typical NTU value is 52. 
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the larger particles. This is essentially the reason, that for spacecraft 
uses, the l/16-in. pellets have been chosen; from the data presented here, it 
seems the diffusional resistance within the smaller pellets is negligible. 

In summary, it appears that the I / 1 6- i n . molecular sieve pellets and the 
6 to 16 mesh silica gel particles have only slight intraparticle diffusion 
resistance, and high values of diffusivity should be used in performance 
prediction programs. The recommended mass-transfer coefficient values in 
Table 5-7 should provide good predictions for spacecraft sorbent beds. The 
recommended mass-transfer coefficient values appear somewhat low compared 
to published correlations; however, nonuni form bed temperature, and axial 
dispersion effects, as well as channeling effects, i f any, are not considered 
in the performance-prediction program's mathematical model and therefore cause 
effective K values to be lower than the correlations. 

A 
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SECTION 6 


COADSORPTION DYNAMIC MASS-TRANSFER TESTS 


GOALS 

The effect on the equilibrium capacity of molecular sieves for coadsorbed 
gases and vapors has been discussed in Section 3, Although the coadsorpti on 
equilibrium data presented in that section is certainly not all inclusive, 
this information provides a quantitative basis for more detailed analyses than 
were previously possible. However, this data does not provide insight as to 
possible dynamic effects accompanying the coadsorption of two or more sorbates. 
One of the goals of the program was to provide information on the dynamic 
behavior of molecular sieves when CO^ and water sorbates are both involved. 

For spacecraft molecular-sieve COg- removal systems, two modes of coad- 
sorption of CO 2 and water are of interest. One mode involves the predrier 
bed where appreciable quantities of both CO 2 and water are contained in the 
gas stream. The second mode concerns the CO 2 '-removal bed; it involves very 
small water-vapor concentrations in the gas stream, but with appreciable 
resident loadings of water. The first coadsorption mode is not considered 
to be very important. Usually the assumption is made that the presence of 
CO 2 has no effect on the predrier bed. There already is a reasonable amount 
of data to support this assumption, at least for design purposes. Consequently 
only a few dynamic runs of this test series involve this coadsorption mode. 

Of much greater importance, is the effect of a resident load of water 
on the CO 2 adsorption process. Therefore, a good deal of test time in this 
program was devoted to dynamic tests involving this coadsorption mode. 

SIMULTANEOUS COADSORPTION OF C0 2 AND WATER ON MOLECULAR SIEVE I3X 

Two dynamic tests were conducted which represent the first adsorption 
half-cycle of a predrier bed with both CO 2 and water vapor in the nitrogen 
carrier gas. The CO 2 partial pressure was 7 mm Hg; the water partial pres- 
sure was 9.2 mm Hg (50°F dew point); and the total pressure was essentially 
7.0 psia. The sorbent in both cases was Linde 1 3X in the form of I /1 6-in. 
pellets. Two runs were conducted; one in the l/2-inch-fin bed, and one in the 
I -inch-fin bed. Both runs were preceded by a vacuum bakeout reducing C0 2 
and water loadings on the sorbent to negligible values. Both runs were 
continued until the breakthrough of both sorbates was complete. 

Figures 6-1# through 6-4 show the water and C0 2 breakthroughs. The 
poisoning of the bed is dramatically illustrated by the CO 2 breakthroughs. 
Immediately after the C0 2 breakthrough is complete, C0 2 removal from the bed 
commences and continues until it is all stripped off. For the test with the 


* Due to the large number of figures, they are included for convenient ref- 
erence at the end of this section. 
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I /2- i nch-f i n bed, the CO2 efflux (as determined by an integration over the 
test period) matches the influx within 1.2 percent; for the l-inch-fin bed, 

CO2 efflux and influx matched within I percent. Thus, there is no net CO2 
retention by the bed. On the other hand, the water capacity is unhindered 
by the presence of CO2. In the l/2-inch-fin bed test, the calculated sorbent 
loading is 0.257 lb h^O/lb sorbent, which is within. 2.2 percent of that pre- 
dicted by equilibrium data. 

To emphasize the contention that there is little effect on the water 
adsorption due to the presence of CO2, data from similar single-sorbate tests 
(water on I3X from Section 4) have been superimposed on the water breakthrough 
curves. For the l/2-inch-fin bed test, where temperatures were relatively 
easily controlled, the two water breakthroughs are quite close together. The 
deviation between the two l-inch-fin water breakthroughs is felt to be due to 
temperature control. It is doubtful that the temperature control variations 
are due to the C02> as the C0 2 adsorption period was long past when signifi- 
cant water breakthrough was noticed. Before water breakthrough at roughly 
four hours, the exit dew points for both tests with the 1/2- inch-fin bed were 
below -87°F (0. 000006 lb h^O/lb N2), For the l-inch-fin bed runs, lowest dew 
points were -82° and -76°F. From these data, there is no discernable effect 
of the C02» 


ADSORPTION OF C0 2 BY MOLECULAR SIEVE BEDS WITH PRELOADS OF WATER 


General Description of the Poisoning Phenomena in C02-Removal Beds 

Even though predrier beds for spacecraft C02-removal systems can reduce 
the water concentration of the process gas stream to very low levels (-100° to 
-60°F dew point, approximately 2 to 45 ppm by weight), inevitably some of the 
water that does reach the CC^-removal bed will be retained there. This of 
course is due to the tremendous affinity of molecular sieves for water. In a 
spacecraft system, with its cyclic sorption/desorption operation, the buildup 
of water on the CC^-removal bed is very slow. This is as expected since much 
of the water adsorbed in one half cycle would be desorbed in the following 
half cycle. But over a. long period of continual operation, a significant water 
load would be attained. Eventually, a cyclic, steady-state condition would be 
reached and no further net water retention would occur. 

As the water load increases, the CC^-removal performance of thn system 
declines. Depending upon the actual design and performance requirements, 
a point can be reached where CC^-removal performance is unsatisfactory. If 
this does occur, provisions would have to be made to discontinue operation 
with the system, and to remove the water load by a vacuum bakeout. 

Ideally, it would be desired to produce a system where CO2 performance 
would be acceptable even with the presence of a steady-state water load. 

This may not be possible; it would probably not produce a practical system 
with respect to size and weight. In any case, it is most desirable to know 
quantitatively how long a given system design will provide satisfactory 
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performance. Up to present, there has been sufficient uncertainty in this 
regard, that all spacecraft systems have been designed very conservatively 
and with provisions for vacuum bakeouts. 

The prediction of the deviation of acceptable performance is very diffi- 
cult. This is due to two circumstances. First, the rate and equilibrium data 
required for such computations are generally insufficient in breadth and usu- 
ally questionable in accuracy. Second, the computation must involve the tran- 
sient adsorption/desorption performance of the system over many, many cycles, 
possibly hundreds of cycles. This requires an extremely accurate numerical 
procedure; and computer costs can be very large. 

Experimentally, of course, the duration of acceptable performance can be 
determined. However, the procedure is complicated and expensive, as controlled 
conditions must be maintained continuously for very long periods of time, up 
to several months. Even if an experiment does determine system performance 
lifetime, little basic information is gained, due to difficulty in relating 
the observed performance conditions to the existing bed water loading 'di stri- 
bution and to general mass-transfer rates. 

As an example of the current state-of-the-art, long duration tests indi- 
cate that the Skylab RCRS would perform satisfactorily for a period in excess 
of 120 days without interruption for a vacuum bakeout. Because of mission 
durations, and to allow a margin of safety, vacuum bakeouts for the RCRS are 
scheduled for every 28 days. The original RCRS design computations were based 
on the most conservative assumptions. One of these was that all water enter- 
ing the C02~removal bed (containing molecular sieve 5A) would be retained by 
it. Simi 1 ar 1 y, conservati ve assumptions were made as to the poisoning effect 
of the retained water. With these assumptions it was relatively easy to pro- 
duce a design that would meet the 28-day performance requirement; and the cost 
of the design process was field to a reasonable level. As indicated by tests 
mentioned above, the actual performance duration far exceeds the requirement. 
Fortunately, as can be seen by the basic bed sizes (7.5 and 9.9 lb) compared 
to total system weight ( 1 95 lb), little penalty is paid for the over design. 
However, as mission durations and crew sizes increase, it will be more impor- 
tant to produce more accurate designs. This will require considerably more 
and better information than that available for the Skylab RCRS design. 


Test Concept 


In order to produce refined data on the dynamic performance of beds with 
resident water loads, a series of controlled tests were devised. The test con- 
cept and conditions were compiled jointly between personnel of NASA Manned 
Spacecraft Center, NASA Langley Research Center, ’and AiResearch. Specific 
test procedures were implemented by AiResearch. Bas i ca 1 1y, the test concept 
was to place a known preload of water on the sorbent bed (either 5A or I 3X) 
then to run a CO 2 adsorption breakthrough. The amount of the preload in all 
cases was to be that which would be i n equi 1 i brium wi th -80°F dew point (at 
0.006 mm Hg p^ q) . Baseline tests were to be run with the preloads placed 



uniformly on the beds. Then, nonuni form preloads, but of the same total 
magnitude, were to be placed on the bed from flow streams of -40°F dew point 
and also of -20°F dew point. COg breakthroughs were then to be run on these 
preloaded beds. Table 6-1 presents summaries of the test conditions for which 
data is reported. 


CO 2 Adsorption on Linde 5A Beds Containing Water Preloads 


Procedure for Producing -80°F Dew Point Uniform Preload . - From Linde 
data, it was determined that the water loading on molecular sieve 5A corres- 
ponding to -80°F dew point, 0.006 mm Hg p u n and 70°F was 8.5 percent. At 

7.0 psia total pressure, the water concentration in nitrogen for 0.006 mm Hg 
p n is only 7.7 ppm by weight. In order to load the- I /2- inch-fin bed con- 
H 2° " : 
taining about I lb of Linde 5A to 8.5 percent water, using the maximum flow 
rate possible in the apparatus (approximately 5 Ib/hr), it would take nearly 
100 days of continuous running if the process gas stream had a -80°F dew point. 
This was clearly an impractical situation. The procedure used involved 
operation of the bed at higher temperature and higher inlet dew point, such 
that the equilibrium loading was also 8.5 percent. From the equilibrium data, 
possible operating conditions lie on the line of constant 8.5 percent loading. 
It was desired to have as high water partial pressure as possible, in order 
to increase water concentration and decrease preloading time. However, very 
high sorbent temperatures were to be avoided due to limitations in the bed 
heating equipment. The conditions of I00°C ( 21 2 0 F ) and 1.6 mm Hg p u n ( I 0 °F 

H 2 U 

dew point) were chosen as they would give the desired loading and were well 
within the capabilities of the equipment. 

The procedure used is described as follows. The bed was first heated to 
essentially 2I2°F by the hot-nitrogen thermal loop, then opened to vacuum and 
pumped down to about 27 y, Hg. From equi librium data, this operation would 
strip off all CO 2 , but would leave a 2.5 percent water load on the bed. Then, 
with the bed still at 21 2 0 F, nitrogen with a 1 3 °F dew point was passed through 
the bed at a flow rate of 4,0 Ib/hr apd pressure of 251 mm Hg. Flow was con- 
tinued for a little over 9.5 hrs. At about 5 hrs into the test it was observed 
that inlet and outlet humidities were essentially equal. The additional test 
time was provided to insure bed saturation. Accordi ng ' to equilibrium data 
with I3°F dew point ( I . 85 mm Hg p u A ) and 2i2°F isothermal bed, this opera- 

H 2° 

tion would produce 8.9 percent loading. It was intended to produce a loading 
above the desired 8.5 percent and then gas strip back to the desired value. 
Calculations from the recorded inlet and outlet dew points indicated that the 
developed loading at that time was 8.76 percent. However, during the preload- 
ing procedure, mechanical problems were experienced with the manual dew-point 
instrument. Because of this, the full water breakthrough curve was not 
recorded and thus the' 8.76 percent load value is only approximate. In view 
of the recorded temperatures in the bed, which varied from I63°F at the inlet 
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to 230°F within the bed, an average of 220°F, the 8.9 percent theoretical 
loading would not be developed. The 8.76 percent value seemed reasonable. 

The bed was then cooled to 70°F, inlet humidity adjusted to -80°F (±2°F), 
followed by several hours of gas stripping. During stripping it was evident 
that very little water was removed from the bed, and the procedure was dis- 
continued. The slowness of the gas stripping is probably due to the extreme) y 
low partial -pressure driving forces involved. 

The preload procedure was not as smooth as desired, especially with some 
data problems during breakthrough. However, the calculated 8.76 percent water 
loading is believed to be valid within about ±5 percent. This tolerance was 
considered acceptable for the test series. 

CO ^ Adsorption by Unde 5A with a -80°F Dew Point Uniform Water Preload . - 

Immediately after the preloading sequence was finished CO 2 adsorption tests 
were conducted on the bed. In all such tests the water concentration in the 
nitrogen carrier gas was adjusted to -80°F dew point ±5°F. CO 2 adsorptions 
were conducted as in the singl e-sorbate test series. After complete CO 2 
breakthrough was obtained, a gas stripping operation was performed to remove 
the CO 2 and prepare the bed for the next run. Calculations showed the gas 
strips quantitatively removed the CO 2 from the bed. 


Figure 6-5 shows the CO 2 breakthrough and gas strip for the flow rate of 
1.23 Ib/hr. The curve shapes — when the gas strip curve is i nverted--are almost 
identical. This is more clearly seen in the dimensionless curves of Figure 
6-6. These curves show a much poorer breakthrough shape than noticed in the 


si ngl e-sorbate tests with CO 2 on dry Linde 5A, 


Figure 6-7 presents the COg adsorption and gas strip for the flow rate of 
2.23 lb/hr. Altogether, 12 runs--adsorpti ons and gas strips--were run at 7 
Hg p^Q ' and 70°F. Very good reproducibility was obtained. The average C0 ? 


mm 


loading obtained for the 12 runs was 0.00406 lb CO 2 / 1 b sorbent; the standard 
deviation was 0.00026. 


This CO 2 loading capacity compares reasonably well with cdadsorption 
equilibrium data taken previously in the program, and shown in Figure 6-8. 

The loading does seem somewhat low with respect to the equilibrium data, and 
it might be speculated that the water preload on the bed was actually more 
than the calculated 8.76 percent. For example, if the water preload were 12 
percent, the CO 2 loading derived from the dynamic tests would be more consis- 
tent with the equilibrium data taken earlier. However, according to Linde 
data, the partial pressure of water in equilibrium with 12 percent water 
loading is 0.025 mm Hg, which corresponds to a dew point of -60°F. The out- 
let dew point noted throughout the dynamic coadsorption tests was always much 
lower than this value, and there was virtually no change noted in later reruns 
after many hours of nearly continuous purging with -80°F dew-point flow. It 
is felt the CO 2 loading derived from the dynamic tests may be more reliable 
than that obtained in.the equ i 1 i br i urn test apparatus. 
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CO ,, Adsorption by Linde 5A with Nonun i form- Water Preloads . - The CO^ 
adsorptions with uniform water preloads do not truly represent the phenomena 
expected in a spacecraft CO -removal bed. However, these tests can furnish 
basic data to help determine mass-transfer coefficients corresponding to the 
8.76 percent water preload. A more realistic simulation would be with the 
resident water loading distributed predominantly at the front end of the bed. 

To accomplish this, a series of tests were conducted in which the preload was 
laid down by nitrogen flows which contained water to the level of -40° and -20°F 
dew points. The total water preloads were to be the same as that in the ~80°F 
dew-point, uniform preload tests; that is, 8.8 percent. Each preload sequence 
was preceded by a vacuum-bakeout to remove all water and COg from the bed. 

After each preload only one adsorption breakthrough was made as it was expected 
further testing would distribute the water preload. During each CO^, adsorption, 
the flow stream dew point was adjusted to -80°F ±5°F. 

Figure 6-9 shows the breakthrough curves derived for the flow rate of 
1.23 lb/hr with nonuniform water preloads . Figure 6-10 shows the preload 
curves for the, -40° and -20°F dew-point tests. Also shown on Figure 6-9 are 
the data points for a duplicate prel oad/C02-breakth rough of Run No. 5A/W-C02 
(l/2)-8. This duplicate run (Run No. 5A/W-C02( I /2 ) -7 ) actually preceded the 
run that is plotted. Because of the ci rcumstances of the duplicate run, the 
comparison between the two is of interest. The preloading sequence for each 
run required 24 hours continuous operation (see Figure 6-IOa), The usual 
procedure was to bakeout the bed overnight, then to start the preload sequence 
by noon of the following day. Then in the afternoon on the day following this, 
the CO 2 breakthrough would be run. On this particular occasion, because of 
laboratory problems (-80°F dew point could not be attained due to anomalously 
high water content of the laboratory nitrogen), the CO 2 breakthrough could not 
be run. Since this was a Friday, it was decided to cool the closed bed at 
approximately 40°F during the weekend. The CO 2 breakthrough was then fun on 
Monday. Because of doubts about the condition of the preload due to the 48- 
hr delay, the entire sequence was rerun. It seems obvious, because of the 
closeness of the breakthroughs, that there was not much migration of water 
within the bed during the idle period. 

It is evident from Figure 6-9 that the distribution of the preload is 
quite important. For the ^20° and -40°F dew-point preloads, it is expected 
that the water would be accumulated near the front of the bed, completely 
poisoning this region. However, the back end of the bed would be left rela- 
tively unpoisoned, capable of much more CO 2 adsorption than if the preload 
were uniform. It is noticed that the -40°F dew-point preload case picked up 
more CO 2 than that for -20°F dew-point preload. The proposed reason for this 
is that at the lower concentration, the bed would adsorb water more toward the 
front end of the bed and allow less to penetrate deeper into the bed. The net 
poisoning effect would thereby be somewhat less, and more CO 2 capaci ty would 
be retained. 

Figure 6-11 presents COg breakthroughs for nonuniform water preloads for 
the flow rate of 2.23 Ib/hr. Figure 6-12 shows the water preload data for the 
-40° and -20°F dew-point sequences which preceded the CO 2 breakthroughs. In 
these tests there is a great di fference in the CO 2 breakthroughs between the 
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-20° and -40°F dew-point preloads. It is not known why there is such a large 
difference here, and yet there was only a small difference for the low flow 
rate tests. The water preloading sequences were nearly identical, as can be 
seen from Figures 6— 1 0 and 6-12. Possibly, wi th the longer CO^ breakthrough 
for. the low flow rate test, the preload could have migrated significantly, j 
becoming more uniform and causing the breakthrough to be earlier, with less; 
total CO 2 capacity. . 

However, the same general situation is felt to exist: that the -20° and 

-40°F dew-point preloads are concentrated at the front of the bed, leaving the 
rest of the bed relatively unpoisoned. Also, since the -40°F dew-point pre- 
load would be more concentrated toward the front of the bed, the CC >2 break- 
through would be longer, and show higher loading. 

Fi gure 6-13 presents CO 2 breakthrough curves on various water preloads 
for conditions of low CO 2 partial pressure, nominally 1.75 mm Hg. These curves 
exhibit the same effects as obtained with the 7 mm Hg p^p tests. Figure 6-14 

shows the water preload data for the -20° and -40°F dew-point sequences which 
preceded the tests shown in Figure 6-13, 

Figures 6-15 and 6-16 present dimensionless CO 2 breakthrough curves for 
the test .series with water-preloaded Linde 5A molecular sieves. Comparing 
the breakthroughs on preloaded beds to the dry breakthroughs, and also com- 
paring the preloaded breakthroughs according to the type of preload, it seems 
that there is a definite mass-transfer rate effect. due to the preadsorbed 
water. Since the flow rate and temperature in all CO 2 tests shown on each 
figure are the same, the gas-phase mass-rtransfer coefficient Kg would be 
essentially the same for all such tests. Therfore, increased 1 ntraparti cl e 
diffusional resistance due to, the water preload is expected as the cause of 
the lower rates. It seems logical that partially filling the molecular sieve 
cages with water could impede the diffusional transfer of CO 2 , or at least, 
produce a competition for adsorption sites that would impede the intraparticle 
adsorption process. On the other hand, ft is possible that the poorer break- 
through shapes could be due to adverse or unfavorable equi 1 i br i urn curve shapes. 
Data is not available to fully investigate this possibility. The possibility 
of greater intra-particle mass-transfer resistance seems much more likely. 


C0 2 Adsorption on Unde I3X Beds Containing Water Preloads 

Procedure for Producing the -80°F Dew Point Uniform Preload , - An im- 
proved procedure, over that used in the 5A series, was employed to produce the 
uniform water preload on the I3X molecular sieve beds. For -80°F dew point 
(0.006 mm Hg water partial pressure) and 70°F sorbept temperature, equilibrium 
data indicates a loading of 7.5 percent. As noted before, loading with -80°F 
dew-point flow stream is impractical; however, the, same loading can b$ 
accomplished with higher dew points and higher bed temperatures. For this 
preload, -2°F dew-point (0.86.mm Hg p H q) and 205°F bed temperature were 
chosen. 2 
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Two problem areas were encountered in the 5A procedure: nonuniform bed 
temperature and ineffectiveness of the gas strip procedure after loading. 
Nonuniform bed temperatures are due to the lack of a specifiec bed preheater 
and the rather high flow rate chosen for the preload sequence. It is noted 
that for I 20°F adsorptions at normal flow rates, experience showed that the 
bed inlet section served as an adequate preheater. In the adopted procedure 
for preloading the I 3X bed, it was decided to start with a high flow rate and 
produce the bulk of the preload, and then to reduce the flow until nearly 
uniform bed temperature was attained. Also, it was desired to load only to 
7.5 percent and forego a gas strip af tferward. 

The I3X preload was preceded by a thorough vacuum-bakeout at 500°F. This 
was to remove all traces of water; thus producing a more well-defined 
starting point than in the 5A procedure. Final calculations of the water 
preload would then be more accurate. 

The first portion of the preloading sequence used a relatively high flow 
rate of 4.0 lb/hr. The heating system was adjusted in an attempt to produce 
a bed temperature profile which averaged 205°F (the gradient actually produced 
was from 180° to 242°F). At these conditions, a water breakthrough was ob- 
tained in about 6-1/2 hours. This is shown clearly in Figure 6-17. After 
the initial breakthrough, the flow rate was reduced to 0.8 lb/hr in order to 
reduce the bed temperature gradient. The resulting temperature profile was 
from 190° to 209°F, with most of the bed at about 205°F. Immediately upon 
the adjustment of flow rate and temperature, the outlet dew point dropped to 
about -20 q F, indicating the average bed temperature for the first break- 
through was considerably above the desired 205 °F, and the bed was not yet 
saturated with respect to this sorbent temperature. Operation continued until 
the outlet dew point essentially reached the inlet at about -2°F. The total 
test time vi/a s over 17 hours. From a graphical integration of the data of 
Figure 6-17/ the water loading derived in the operation was calculated as 
7.5 percent, just as desired. 

00 ^ Adsorption by Linde 13X with Uniform and Nonuniform Water Preloads . - 
Figures^- 1 8 and 6—19 show the CO^ breakthrough curves obtained from the 
uniformly p re loaded bed. Each CO^ adsorption was followed by a gas stripping 
operation. The water preload was essentia 11 y unaffected by these operations 
as the dew point of the gas stream was maintained about -80°F. A rerun of 
the initial CO 2 breakthrough was made in order to check the condition of the 
preload. As shown in Figure 6-1 9, the rerun breakthrough was essentially 
identical to the first run. The CO 2 adsorption capacities for the three 
adsorptions and the three gas strips were all quite close, averaging 0.857 
percent. From the coadsorption equilibrium data, the COg loading for 7.5 
percent water ' preload is interpolated to be I -.15 percent.. No definite expla- 
nation for the ’ difference between dynamic capacity and the equilibrium data 
is available. It is possible that the equilibrium data’, which is rather 

difficult to obtain, is biased high. 

: .1 ' ' - 

Figures 6-20 1 and '6-22 show CO 2 breakthrough curves obtained for non- 
uniform water preloads. Figures 6-21 and 6-23 show the data from the water 
preload sequences involving -40 p and -20°F dew-point flows. Those preload 
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sequences were essential ly identical to those used in the 5A tests. There is 
one exception; for the preloading run I SX/W-COgC I /2)-7 (Figure 6-23) a lower 
carrier gas flow rate was used, 4.0 vs. 5.2 lb per hr. 

The general trend of the data for I3X is the same as that for 5A. That 
is, the beds loaded uniformly with water in equilibrium with -80°F dew point 
show considerably less CO 2 capacity than dry beds. Beds loaded to the same ' 
extent (or more), but in such a manner as to concentrate the water at the 
front end of the bed, show CO 2 capacities intermediate between those of the 
uniformly preloaded beds and dry beds. One difference seems to be apparent. 

For I3X, the -20® F dew-point preloaded bed shows greater capacity than the -40°F 
dew-point preloaded bed. The opposite was true for 5A. The suspected cause 
of this difference lies in the -20°F dew-point preload on I3X being derived 
from the lower flow of 4.0 lb per hr. It is assumed that with the lower flow, 
in spite of the higher inlet concentration, the water was adsorbed more pre- 
dominantly ; on the front end of the bed, with relatively less adsorption else- 
where in the bed. 

Dimensionless CO 2 breakthrough curves are shown for all coadsorption 
tests on I3X in Figures 6-24 and 6-25. From the figures it is apparent that 
the water preload has a lowering effect on the mass-transfer rate. As noticed 
for the 5A tests, both the magni tude of the preload and the preload distri- 
bution are factors in how much the rate is affected. 

It seems that I 3X sorbent is less effected in rate by water preloads than 
5A. This conclusion is based upon a comparison of Figure 6-24 and Figure 6-15. 
In Figure 6-15 for 5A sorbent, the curves for the preloaded beds are relatively 
less steep when compared to the dry-bed curves. As discussed earlier for the 
5A tests, the general explanation for the lower mass-transfer rate is that the 
pre.oaded water significantly restricts movement of CO 2 molecules within the 
zeolite structure. Since I 3X Has a much more open crystal structure, it is 
expected that the preloaded water would indeed have less effect than with 5A. 

Prediction of Water Preloads Derived from -20° 
and -40° F Dew-Point flows 

Performance-prediction program S9960 was used to predict the nonuniform 
water loading distributions for four of the preload sequences which involved 
-20° and -40°F dew-point flows. Two of the computations were for the Linde-5A 
series. Runs SA/W-CO^ 1/2 ) - 9p and -I2p; the other two computations were for 

the Linde-l3X series, Runs l3X/W-C0~(l/2)-7p and -9p. The best predictions 

2 2 

were obtained for K = 0.0004 lb-moles/hr- ft -°F and D = 0,001 ft /hr. Predic- 

9 

tion accuracy was judged by comparing the total bed 'loadings from the computer 
predictions with that calculated from experimental data. 

Figures 6-26 and 6-67 present the predicted water loading distributions 
(water loading vs axial position in the bed). As suggested earlier, the figures 
show that the preloads derived from -20° and -40°F dew-point flows are concen- ' 
trated in the front of the bed, with the back of the bed relatively' unloaded. 
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From Figure 6-26, showing the preloads on Linde 5A, it is difficult, to see 
why the CO^ breakthroughs conducted after the preload sequences (shown earlier 

in Figure 6-ll) are so different (the -40°F dew-point preload breakthrough 
being much better). Possibly, the additional amount of, water present migrated 
considerably during the CO^ breakthrough for the -I8°F dew-point-preload case, 

and caused considerable poisoning throughout the bed. 

On the other hand, for the preloads on Linde I3X, shown in Figure 6-27, 
it seems that the -20°F dew-point preload would yield somewhat more CO 2 capa- 

ci ty than that for -40°F dew point. Figure 6-22 does show somewhat higher 
capacity for the -20°F dew-point run. The very small amount of loading in 
the rear of the bed is considered to give the greater observed capacity. It 
should be noted that at about 5 percent water loading, the CC^ absorption 

capacity is a less than half of that of dry molecular sieve; at 7 percent it 
is about one quarter or less. Thus, it would seem that the front of the bed 
is completely poisoned and the loading in the rear of the bed determines the 
remaining capacity. 

Prediction of Mass-Transfer Parameters 
on Uniformly Preloaded Beds 

An attempt was made to evaluate mass-transfer parameters from the COg 

breakthrough curve of Run No. SA/W-CC^ I /2)-5 on the uniformly water-preloaded 

Linde 5A bed. In order to do this, the equilibrium isotherm for the adsorp- 
tion conditions had to be obtained. Experimental data was not avai 1 able ; 
therefore the isotherm was approximated. This was accomplished as follows: 
first, it was fairly certain from dynamic tests that at ebout 7.1 mm Hg p rn 

LU 2 

the bed equilibrium loading was approximately 0.00406 lb COg/lb 5A. This was 

discussed earlier in this section. Then from the data contained in Figure 3-23 
this equilibrium point was located and the slope for the isotherm was approxi- 
mated. It was then estimated that the isotherm would pass through the point 
of w = 0.001 lb C 0 g / 1 b 5A and p^ = 0.4 mm Hg. On log- log paper a straight 

line was drawn through these two points; this straight line was the isotherm 
for which tabled values were included in the performance-predi ction program. 

The evaluation yielded one combination of K /D which gave a good fit of 

9 2 

the experimental curve. The values were K = 0.0001 Ib-mole/hr-ft -mm Hg 
2 9 

and D = 0.0001 ft /hr. Figure 6-28 shows the prediction obtained with these 
values. From all of the predictions, it was fairly definite that diffusivity 

2 

values of 0.00004 and 0.001 ft /hr were too low and too high, respectively. 
There Seemed to be some latitude in the K value, possibly down to 0.00075 

2 ' 9 

1 b-mol e/hr- f t -mm Hg. 
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Higher K 


9 


va 1 ues 


than shown in Figure 6-28 


tended to predict later 


breakthrou 

K values 
9 


ghs in the early portion of the run; however, predictions for all 
for a particular diffusivity seemed to merge into one line for the 


later part of the breakthrough. 


No other evaluations were attempted due to the obvious approximate 
nature of the equilibrium data. This one evaluation does serve to point out 
that rate coefficients are reduced to some extent when the sorbent contains 
a water preload. 
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Water concentration, t- Ob H.O/Jb Nj “H C0 2 concentration C CQ (Jb/CO^/Ib Ng) 
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Bed fin surface: l/2-ln. 

Flow rate: 2.30 Ib/hr 

Bed temperature: 60°F 

Inlet dew point: 50°F 

Inlet water partial pressure: 9.1 

Inlet C0 2 partial pressure: 7.2 

Total pressure: 36B mm Hg 


Material balance error: -J.2S& 
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igure 6-1. CO 2 Concentrations for Simultaneous Adsorption of CO 2 
and Water on Linde 13X, Run No. l3X/C02“W(l/2)--l 


Bed fin surface: l/2-ln. 

“Flow rate: 2.30 lb/hr H 

Bed temperature: 60°F 

— Inlet dew point: 50°F H 

Inlet water partial pressure: 9.1 mm Hg 

__ Inlet COg partial pressure: 7.2 mm Hg ii 

Total pressure: 368 mm Hg i 

_ Water influx: 0.319 lb j 

Water efflux: 0.056 lb ii 

__ Water left on bed: 0.263 lb 

Bed mass: I .004 lb 

: Water loading: 0.262 lb/lb 1 

...1 Loading from equii. data: 0.256 j 

Deviation from equii: +2.2# j 

’# Data points from run !3X/W(l/2)-3 without^ 
the presence of C0g •] 
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Figure 6-2. Water Concentrations for Simultaneous Adsorption of CO 2 
and Water on Linde I3X, Run No. !3X/C0 2 -w(l/2)-l 
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Figure 6-3. COg Concentrations for Simultaneous Adsorption of CO 2 
and Water on Linde 13X, Run No. I3X/C02-W(I )-l 
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Figure 6-4. Water Concentrations for Simultaneous Adsorption of CO 2 
and Water on Linde I3X, Run No. 1 3X/C0 2 ~ w( I ) -.1 
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CO, concentration, C_ n (CO./Ifa N«) COg concentration 



Time (minutes) 

CO^ Adsorption Breakthrough, on Linde 5A With 8.76$> Uniform 
Water Preload, Run No« 5A/W“C0*( 1 / 2 ) — I 



Time (minutes)' 


B. Gas Trip of COg From Linde 5A with 8.76^ Uniform Water Preload 
Run No, 5AW-C0 2 (I/2)-2 s 

Figure 6-5. C0 2 Adsorption Breakthrough Run No. 5A/W-C0 2 ( 1/2)- I and 

Gas Strip Run No. 5A/W-C0 2 ( l/2)-2s 
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Figure 6-6. Dimensionless Breakthrough Curves; CO- on Linde 5A, 8.76$ Uniform 
Water Preload 
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COg Adsorption Breakthrough on Linde 5A With 9.76% Uniform Water Preload, 
Run No. 5A/W-C0 2 (l/2)-5 
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Lj Flow rate: 2.23 Ib/hr 

n Bed temperature: 70°F 

Total pressure: 365 mm Hg 

FjCalc. initial sorbent loading: 0.00393 lb/)b 

F Bed mass: I .257 1b 

\ Initial sorbate on bed: 0.00495 lb 

I Net sorbate removed: 0.00521 lb 

f Deviation of material balance: 5.3i 
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B. Gas Strip of CO^ From Linde 5A With 8.76 % Uniform Water Preload, 
Run No. 5A/W-C0 2 (l/2)-6s 

Figure 6-7. CO 2 Adsorption Breakthrough Run No. 5A/W-C02( I /2) -5 and 
Gas Strip Run No. 5A/W-C0 2 (l/2)-6s 
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Figure 6-8. Water/CO^ Coadsorption Equilibria for Linde 5A Molecular Sieve 
Pellets » 
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Figure 6-9. C0 2 Breakthrough Curves on Water P re loaded Linde 5A Molecular Sieve 

(1-23 ib/hr) 
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" Bed temperature: 70°F 

r Inlet partial pressure: 0.311. mm Hg 
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A. Water Preloading Sequence on Linde 5A, Run No. 5A/W-C0„( l/2)-8p 
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Bed fin surface: i/2**in. 

Flow rate: 5.19 Ib/hr 

Bed temperature: 70°F 

Inlet partial pressure: 0.096 mm Hg 

Inlet dew point: -40° F 

Total , pressure: 142 mm Hg 

Water influx: p. 1112 lb 
Water efflux: — 

Water left on bed: 0. M12 lb 
Bed mass: 1.257 

Water preload: 8.84# 
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B. Water Preloading Sequence on Linde 5A, Run No. 5A/W-C0„( I /2 ) - 1 Ip 


Figure 6-10. Water Preloading Sequences on Linde 5A Run Nos. 5A/W-C02( l/2)-8p 


and 5A/W-C0 2 (l/2)-l Ip 
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Figure 6-11. C0 2 Breakthrough Curves on Water Preloaded Linde 5A Molecular Sieve 

(2.23 lb/hr) 
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Bed fin surface: 1/2- in. 

Flow rate: 5.17 Ib/hr 

Bed temperature: 70° F 

Inlet partial pressure: 0.33 

Inlet dew point: -I9*5°F 

Total pressure: 228 mm Hg 

Water Influx: 0.1229 lb 
Water efflux: — 

Water left on bed: 0.1229 lb 
Bed mass: 1.257 lb 

Water preload: 9.78$ 
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Water Preloading Sequence on Linde 5A, Run No. 5A/W-C0_( I /2)-9p 
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Bed fin surface: 

l/2-ln. 


Flow rate; 5.22 

Ib/hr 


Bed temperature: 

70° F 


Inlet partial pressure: 0, 

,093 mm Hg 

Inlet dew point: 

-40. 5° F 


Tow l pressure: 

145 mm Hg 


W£ter influx: 0. 

1080 lb 


Water efflux: — 
Water left on bed 

: 0.1080 

lb 

Bed mass: I .257 

lb ' 


Water preload: 8.59$ 




Time (hrs) 

B. Water Preloading Sequence on Linde 5A, Run No. 5A/W-CC> 2 ( I /2)- 1 2p 


Figure 6-12. Water Preloading Sequences on Linde 5A, Run Nos. 5A/W-C0„( !/2)-9p 
and 5A/W-C0„( 1/2)- I 2p 
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Figure 6-13. CO 2 Breakthrough Curves on Water Preloaded Linde 5A Molecular 
Sieve at Low C0 ? Partial Pressure 
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Bed fin surface: 1/2 in. 

Flow rate : 5. 1 7 1 b/hr 

Bed temperature: 70°F 

Inlet partial pressure,: 0.31 1 mm Hg 

Inlet dew point: -20.5°F 

Tota l pressure: 228 wn Hg 

Water influx: 0.1070 lb 

Water efflux': 0.0001 lb 

Water left on bed: 0.1069 lb 

Bed mass; 1.257 lb 
Wa te r p re load : 8. 50# 
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Water Preloading Sequence on Linde 5A, Run No. 5A/W-C0- ( I /2 ) - 1 Op 
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Bed fin surface: 1/2 jin. 

Flow rate: 5.19 Ib/hr 

Bed temperature: 70°F 

Inlet partial pressure: 0.096 mm Hg 

Inlet dew point: -40°F 

Total pressure: 142 mm Hg 

Water influx: 0.1083 lb 

Water efflux: •— 

Water left on bed: 0.I08J lb 

Bed mass: 1.257 lb 

Water preload: 8.62# 
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Water Preioading Sequence on Linde 5A, Run No. SA/W-CO^C I /2) - I3p 


Figure 6-14. Water Preloading Sequences on Linde 5A, Run Nos, 
5A/W-C0- ( I / 2 ) — ,1, Op and 5A/W-C0„( l/2)-!3p 
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Figure 6-15. Dimensionless CO 2 Breakthrough Curves on Water Preloaded Linde 5A 
Molecular Sieve (1.23 lb/hr) 
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Figure 6-16. Dimensionless CO 2 Breakthrough Curves on Water Preloaded Linde 5A 
Molecular Sieve (2.23 lb/hr) 
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Figure 6-17. Water Preloading Sequence on Linde I3X, Run No. 1 3XW— COg ( I /2) — I p 
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B. Gas Strip of CO 2 From Linde I3X With 7.5 fo Uniform Water 
Preload, Run No. I3X/W-C0 2 ( I /2)-2s 


Figure 6-18. C0_ Adsorption Breakthrough Run No. I SX/W-COg C • /2') - 1 

and Gas Strip Run No. !3X/W-C0 2 (l/2)-2s 
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A. COg Adsorption Breakthrough on Linde I3X With 7.5 $ Uniform 
Water Preload, Run No. I3X/W-CC> 2 ( I /2)-3 
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Figure 6-19. C0 2 Adsorption Breakthrough Run No. I 3X/W-C0 2 ( I /2)-3 

and Gas Strip Run No. I 3X/W-C0 2 ( l/2)-4s 
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Figure 6-20. C0 2 Breakthrough Curves on Water Preloaded Linde I3X 

Molecular Sieve (1.23 Ib/hr) 
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Figure 6-21. Water Pre loading Sequence on Linde I3X, Run No. ISX/W-COgC I /2)-8p 
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Water concentration, C u n (lb H_0/ 1b N.) -»> Water concentration 
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Figure 6-24. Dimensionless COg Breakthrough Curves on Water Pre loaded 
Linde I3X Molecular Sieve (1.23 Ib/hr) 







o 

CM 



4-*|CD 


V) 

(0 

<U 


V) 

c 


^ 'uoi }ej*U9ouoo ssa^uoi suaoij a 


181 


Figure 6-25. Dimensionless CO^ Breakthrough Curves on Water Preloaded 
Linde I3X Molecular Sieve (2.33 lb/hr) 
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Figure 6-28. Evaluation of Mass-Transfer Parameters for Uniformly Water- 
Preloaded Bed of Molecular Sieve 5A, Run No. 5A/W-C0 9 ( I /2)-5 




mm BLmK 


not FILMED 


SECTION 7 

PERFORMANCE PREDICTION COMPUTER PROGRAMS 


S9960, TRANSIENT PERFORMANCE-PREDICTION 
FOR C0 2 REMOVAL WITH SINGLE OR COMPOSITE SORBENT BEDS 

Program S9960 for water- dump/COg- dump compos i te sorbent beds (originally 
produced by Ai Research under another NASA contract) was the first program 
delivered to NASA Langley Research Center under this contract. Under this 
contract, some updating of the program was made and full documentation was 
prepared (Reference I ) . It has been checked out thoroughly on CDC-6600 and 
Univac 1108 computer systems. 

This program deals with systems such as the Skylab RCRS where both water 
and COg are dumped to vacuum. The program optionally allows for internal 
coolant passages, or adiabatic beds can be simulated. Any sorbents can be 
employed in the ; desiccant and CC^-removal sect! ons of the bed. Vacuum desorp- 
tion is the main desorption mode, but gas stripping simulations can be made 
by proper input data assignments. 

The program was used exclusively in the mass-transfer coefficient evalua- 
tions as presented in Section 5. In this case, by input data assignments, 
the bed was configured to Hold only one type of sorbent and only one sorbate 
was contained in the process gas stream. Also, for the evaluation, the par- 
ticular option was specified which caused the predictions to be made for 
isothermal beds. If warranted, bed heat transfer could have been simulate^. 
However, this would have introduced a great deal of complexity into the mass- 
transfer coefficient evaluations. For the evaluations, a specific main pro- 
gram and printed output program were devised. These programs make it easier 
to run successive breakthrough predictions with different Kg and D values; 
these new programs are given later in this section. 

MAIN 4B, PERFORMANCE PREDICTION PROGRAM 
FOR C0 2 -REM0VAL SORBENT SYSTEMS WITH MULTIPLE BEDS 

MAIN 4B i s an extended version of program S9960. All of the refinements 
and additions to MAIN 4B were made under this contract. The documentation 
of this extensive program is contained in Reference I . 

Many of the numerical schemes used by MAIN 4B are essentially the same as 
S9960, and it uses many subroutines written earlier. Several features are 
included in the program that were not considered in S9960: 

® The poisoning effect of gradual water buildup on the COg sorbent 

• Estimation of overboard gas losses, oxygen and nitrogen 

o Two- and four-bed capability 
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® An option to bllow electric heater inclusion in the bed 

© The capability to analyze desorption when vacuum pumping is included 
in the system f 

By considering ( I ) the poisoning of the COg- removal bed by water and (2) 
the adsorption and subsequent loss of atmospheric gases Og and No* this program 
treats some of the very important second-order effects due to coadsorption. 

Even though the program does not handle the coadsorption in a thoroughly 
rigorous manner, it should provide good results, which will prove quite valu- 
able in system design and optimization work. The program is much more 
versatile then S9960 in that it can handle two- or four bed systems with 
water/CO^-dump or -save provisions. 


PROGRAM ADDITIONS TO S9960 


Revised Main Program S9960 and Revised Print Subroutine PRADSB 

Appendices A, B and C present computer programs which are additions to 
the program series S9960. The new programs have been written to facilitate 
the evaluation of mass-transfer coefficients from the experimental dynamic 
data as discussed in Section 5. In Appendix A, a new version of main program 
S9960 is presented. This program is used to run back-to-back adsorption half- 
cycle predictions, each starting with the same initial conditions, but each 
with a differing set of mass-transfer parameters Kg and D (Fortran variables 
GK and DIF), The program is intended to be used with the new PRADSB subroutine 
for pri nted output presented in Appendix B, It is, however, compatible with 
the original PRADSB routine (S9979), which is included in Reference I. 

The new main program requires some input by means of punched data cards. 
The first card of the data deck is a heading card. The information on this 
card is printed on the first line of each page of output during the entire 
computer run. Following the heading card in the data deck are the cards which 
specify the values of GK and DIF to be applied to the different predictions. 

One pair of GK and DIF values are contained on each card; the individual GK 
and DIF values (named locally as GKR and DIFR) are to be applied to all nodes. 
Thus, this main program can only be used where GK and DIF are desired to be 
constant throughout the bed. In order that the second and successive predic- 
tions start with the same initial conditions as given i’n the BLOCK DATA pro- 
gram (S9993), this new main program stores the initial conditions on a 
scratch device and restores them to core at the start of each successive run. 

It is also noted that the card input data for GK and DIF override the data 
contained in the BLOCK DATA program. 

Also to be' punched on each data card (in the third ten-column field) is 
the integer indicator IPRINT. This indicator when zero- or blank, instructs 
the new PRADSB program to print only outlet quantities vs sorption time, 
giving a tremendous savings in printer time and number of pages, when only 
this information is desired. When greater than zero, quantities for all 
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finite-difference elements throughout the bed are printed by PRADSB, but with 
a more compact format than in the earlier version. In addition to format 
changes, the new PRADSB gives values for gas-phase concentration for both 
water and as well as the partial pressures. The concentration values 

are those defined by Equations 4-1 and 4-2. For user convenience, values of. 
GKR and DIFR are pr i nted wi th the headi ng on each page of output. 

The new version of PRADSB can be used as is with the existing program set 
described in Reference I. However, hea.ding information and the IPRINT indi- 
cator should be added to the BLOCK DAJA program if one wants to use the full 
capabi 1 i ti es of the new output program. • 


Equilibrium Data Subroutine PKEQ 

Appendi x C lists the new, table- lookup PKEQ. subroutine 
(for sorbent/sorbate equilibrium data). This routine has been used exclusively 
in the mass-transfer coef f i c i ent eval uati ons of Secti on 4. The table-lookup 
concept was resorted, to when it became extremely di ff i cul t to obta i n 
equilibrium-data curve fits with suff i ci ent accuracy for these evaluations. 

In this regard, the program has worked very wel 1 ; and without needing to 
curve-fit each equilibrium i sotherm, it is much easier to prepare equilibrium 
data for use. 

The method for using the program* and requirements on the tabled equili- 
brium data are contained in the comment cards i ncl uded in the Fortran listing. 
As an example of the storage of a single equi 1 i bri urn i sotherm, , the listing in 
Appendix C contains data for COg on Davi son 5A molecular sieve at 70°F. Pro-' 
gram computati ons are about 30 percent s 1 owe r using the tabl e-lookup versi on 
of PKEQ, but. the ease of use, when adequate curve fits are not already avail- 
able, more than justifies the extra computation cost- 

i 

EQUILIBRIUM DATA CURVE FITTING PROGRAM 

i ) ' ; 

In order to use program set. S9960, the user. must produce a Fortran sub- 
rout i ne PKEQ whi ch upon cal l, for the sorbents of i nterest', wi 1 1 yield sorbate 
equilibrium partial pressure as a functi on of sorbent loading and temperature. 
Thi s is di scussed fully in Reference I . The program di scussed above fulfills 
this need w,i th little setup bei ng required of the program user. However, from 
the view poi n't of computation speed, a cl osed-form al gebrai c expression is 
desired--in effect a .curve fit (or a map fit) of the equilibrium data for the 
sorbent/sorbate pair. .Exper i ence has shown that it is often quite difficult 
to obtain such a curve or map fit with suff i ci ent accuracy for design purposes. 
Many of the classical i sotherm express i ons (Langm.ui r, Freundl i ch, BET, etc. ) 
have been tried; the usual results are curve fits that may have the correct 
shapes* but are not accurate enough. More success has been achieved with 
correlation equations whi.ch may resemble the theoreti cal isotherm expressions, 
but whi ch possess on.e or more arbi trary groupingsof terms and parameters . 
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In order to facilitate the curve-fitting of equilibrium isotherms and 
maps, a special program has been developed. This program and its method of 
use are presented in Appendix D. The program contains 35 different correlation- 
equation forms. The program produces coefficients and exponents for the corre- 
lation equations, and shows the deviation between the predicted values and the 
data used for input. If one is to be concerned with a large amount of 
performance-predi cti on computations for certain sorba te/sorbent pairs, the use 
of this program is suggested. If sufficiently adequate correlations result, 
computer time is saved. If correlations are not considered adequate, the 
table lookup version of PKEQ, can be used. • 
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SECTION 8 


PRESSURE DROP AND PACKED-BED DENSITIES 


PRESSURE DROP 

The knowledge of pressure drop sustained by the process gas stream flowing 
through a packed bed is an essential piece of information needed in the design 
of any fixed-bed sorbent system. All of the sorbents studied in this program 
have been in large-scale commercial use for several years, and because of this, 
considerable information is available on pressure drop. However, because of 
the plate-fin heat-transfer surfaces incorporated in the test beds, a situa- 
tion rarely if ever encountered in the chemical process industry, it was 
desired to provide accurate experimental pressure drop data. It should be 
noted that the effect of the plate-fin surfaces was expected to be relatively 
small. Generally, it would seem that the additional wall effect in the plate- 
fin beds would tend to produce slightly less pressure drop. 

In all adsorption and gas stripping runs, bed pressure drop was measured 
by a precision Wallace and Tiernan differential pressure gage. With silica 
gel (Davison, Grade 05, 6 to 16 mesh particles), the range and number of ex- 
perimental tests provided enough data to characterize pressure drop performance 
of the material. With Linde 5A and I3X molecular sieves ( I /i 6-in. pellets), 
the experimental range was not considered as great as desired, so a separate 
pressure-drop test series was conducted. With these tests, both an inclined 
oil manometer, and a Baratron electronic pressure indicator were added to the 
instrumentation. For Davison 5A molecular si eve materi al , so little adsorption 
testing was accomplished, that no pressure drop data is reported here. As 
spherical beads, this material packed well in the test beds, and it is expected 
that pressure drop can be predicted easily by using correlations for packed 
beds of spheres. One such correlation, due to Ergun, is given by Davison 
(Reference 9). 

In all of the information that follows it should be noted that the pres- 
sure drop of the empty bed, including that of the two 25-mesh sorbent-retention 
screens, was measured to be so low that it was ignored. That is, the pressure 
drop as measured is due entirely to the flow through the sorbent packing. 
Because of the inorganic nature of the sorbents, there was expected to be no 
effect on pressure drop due to the amount of sorbate held by the bed; no such 
effect was ever observed. In other words, there is no geometric change 
(swelling, shrinkage, etc) of these sorbents as they pick up water or CC^- 

The pressure drop data for silica gel was correlated so well by the 
expression given by the Davison division of W.R. Grace and Co. (Reference 10), 
that no data curves are presented here, The Davison pressure drop equation is 

KV m 

where AP/L is the pressure drop in psi per foot of bed depth, and V, the 
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superficial velocity, is in ft/sec. It should be noted that in calculating V 
for the plate-fin beds, the bed superficial free- flow area from Table 4-12 
was used. This area does not include the cross-sectional areas of the coolant 
passages, the plates, or the gas side fins. It is thus a net superficial area. 
For air, Davison gives for the coefficient K the value 0.T0; for the exponent 
m, 1.5. These values correlated well the nitrogen pressure drops in the 1/2- 
i n. -f i n bed at 5 psia. 


The test data for Linde 1/1 6-in. molecular-sieve pellets are shown in 
Figure 8-1. These data were actually taken on I3X molecular sieve; but since 
all of the Linde I / 1 6— in. pelleted materials appear dimensionally similar, 
these data should be valid for other molecular sieves of this configuration. 

The data for Figure 8-1 cover the pressure range from 5 to 14.7 psia (258 to 
760 mm Hg) and the bed flow rate of 0.5 to 5.4 Ib/hr. The ordinate of the 
figure is pressure drop (in. 1^0) per unit of flow length (ft) times the den- 
sity ratio o, which relates the actual gas density in lb/ft^ to that at stan- 
dard conditions, 0.07176 lb/ft^. The abscissa is flow rate per unit superficial 
free-flow area ( 1 b/hr-f t^) . Within the accuracy required for design, all of 
the data for both beds, at all of the test pressures, could have been repre- 
sented by one line. 


PACKED-BED DENSITY 

Because of the plate-fin heat-transfer surfaces incorporated into the 
sorbent beds, packed-bed density (sometimes called bulk sorbent density) would 
be expected to vary with fin size and would not coincide with predi cted val ues 
which are based upon the large beds as found in industry. To investigate the 
effect of fin size on packed-bed density, four small bed modules were construc- 
ted. Three of the beds contained offset fin surfaces of the type used in the 
test beds discussed in Section 4; i.e., l/4-in. l/2-in.,. and l-in. fin height 
and spacing. The fourth bed was left empty. The modules were approximately 
2 in. by 2 in. at the face, with a 4 in. bed depth. To facilitate the packed 
density measurements all of the beds were welded closed at one end. 

Using the general procedure for packing described in Section 4, the beds 
were packed with sorbent and weighed in a nitrogen dry box. After the various 
packings and weighings were completed, bed free volumes were determined by 
filling with water and weighing. 

Table 8-1 presents packed bed densities as a function of fin size, as 
derived from the bed modules and from the dynamic test beds. For reference, 
bulk density from manufacturer's literature is also given. There is some 
scatter in the density measurements given in the table, but a definite trend 
with fin size is evident. Density variations for the dynamic test beds stem 
from improved packing techniques as the program progressed, and the possibility 
that the retention screens were bulged, thus allowing more sorbent to be held. 
The densities of the bed modules have not been corrected for nitrogen adsorp- 
tion; this could be about I percent for the molecular sieves. Presumably, all 
sorbents used for the bed modules were dry due to bakeouts previous to packing. 
However, it is possible that some significant water pickup could have occurred 
during the packing test procedure. 
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Figure 8-1. Pressure Drop in the I /2-in. -Fin and I -in. -Fin Test 
Beds for Linde i / 1 6 in. Molecular Sieve Pellets 



TABLE 8-1 


PACKED-BED DENSITIES IN FINNED SORBENT BEDS 


Fin Size and Bed 

Silica Gel 
(Davison Grade 05) 
6 to 1 6 Mesh 

Linde I3X 
1 /I 6-i n. 
Pel lets 

Linde 5A 
1 / I 6— in. 
Pel 1 ets 

Davi son 5A 
8 to 12 
Mesh Beads 

l/4-in.-fin bed module 

38.8 

39.3 

41.8 

38.5 

1 /2-in. -fin bed module 

41.7 

43.1 

45.2 

4?.l 

l/2-in.-fin test bed 

43.0 

39.7 

41.1 

43.6 

49.7 

43.0 

l-in.-fin bed module 

44.0 

43.3 

46. 1 

44.4 

1 -i n. -f i n test bed 

43.2 

39.0 
43 . 2 

' 

1 

l 

-- 

Empty bed module 

44.6 

46.9 

48.6 

47.8 

Manufacturer's bulk 
density 

45 

40 

43 

44 
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SECTION 9 


CONCLUSIONS AND RECOMMENDATIONS 


The efforts under this contract have provided two large-scale computer 
program sets to predict the transient performance of spacecraft, regenerable, 
CC^-removal systems which employ molecular sieve adsorbents. The first pro- 
gram for two-bed, water- dump/C02 _ dump systems (Reference I ) has already seen 
extensive use in the design and evaluation of actual spacecraft COg-removal 
systems. The second program (Reference I ) provides a significant amplifica- 
tion and update of the first program. It can be used for two- or four-bed 
systems with either water-, or CO 2 -, save or dump provisions. Further, the 
new program treats some coadsorption phenomena, (poisoning and atmospheric 
gas losses) which, although considered second-order effects, can dominate the 
design of such systems. Several other refinements and additions are included 
in this well checked-out program. 

The equilibrium behavior of several sorbent/sorbate systems was inves- 
tigated in one of the latest, most wel 1 -equi pped, experimental setups of its 
kind. The investigations also produced data on the equilibrium behavior of 
molecular sieves when two sorbates are coadsorbed. Very little information 
of this, nature was previously available. 

The dynamic behavior of specially-constructed molecular-sieve sorbent 
beds was studied in a specially designed experimental setup. The test beds 
contain very efficient plate-fin, offset, heat-transfer surfaces, designed 
for sorbent bed applications. Differing sizes of fins were evaluated in the 
test program. In general, fin height and spacing of 1/2 in. provides good 
bed heat-transfer performance, while allowing good packed bed density. Al- 
though not tested, 3/8-inch fins should provide good heat- and mass-transfer 
performance. Fins as large as I -inch high, spaced I inch apart, should be 
avoided. To limit channeling to a sufficiently low amount, offset fins should 
be provi ded. 

On the basis of si ngl e-sorbate breakthrough tests, mass-transfer coeffi- 
cient and intraparticle di ffus i vi ties were evaluated. For all tests involving 
molecular sieves 5A and I3X, and silica gel, intraparticle diffusion resis- 
tance is quite low. Therefore, large values of effective diffusivity should 
be used in performance-prediction programs. For the systems studied, mass- 
transfer coefficient values are recommended. These values seem somewhat low 
compared to published correlations. However, it is considered that the low 
values compensate for unavoidable non isothermal effects (at the mass-transfer 
front) in the dynamic tests, and the phenomenon of axial dispersion which is 
not considered in the performance-prediction programs. These programs, when 
used with the recommended mass-transfer parameters, should provide good 
transient performance predictions, over a wide range of bed geometries and 
condi ti ons . 

The effect of coadsorbed water and CCU on the dynamic behavior of 
molecular sieve beds was studied. For predrier beds, where water pickup 
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is of interest, the presence of CO 2 in the gas stream may be ignored. For 
C02-removal beds, the presence of water, even in small quantities, has a large 
and important effect. Equilibrium capacity for C0£ is drastically lowered I 
when water is coadsorbed; being one-half of the dry value when water loading! 
is about 4 percent; and practically nil when water loading reaches 8 percent; 
Further, as water is coadsorbed, the CO 2 adsorption rate coefficients declinfe 
markedly. It is considered that diffusion in the sorbent pellets is hindered 
by the presence of water. 

At the present time, techniques for design and performance evaluation 
of regenerable C 02 -removal systems using molecular sieves are highly developed. 

One can feel confident that molecular sieve sorbent systems will meet the 
stringent objectives of space hardware, including low weight and volume pen- 
alties and reliable long life. Much of the information and many of the tools 
discussed in this report have already been used in the improvement of both 
C02~removal systems and the techniques used in their design. Still, with 
changing spacecraft conditions, there are areas where continued study and 
development should be undertaken. At the present time, with lower desired 
levels of CO 2 partial pressure (3 mm Hg) and higher cabin pressures, full 
system design studies should be implemented. It is acknowledged that 
molecular-sieve sorbents do not have large CO 2 capacities at such low partial 
pressure, but as demonstrated in this program with tests involving 1.75 mm Hg 
P^q , molecular sieves show excellent dynamic performance even in this range. 

Also, many inherent qualities of molecular sieves (their inertness, dimensional 
stability, tolerance to vacuum, etc) are not affected by the newer conditions. 

Full system design/tradeoff studies would be enlightening, especially in view 
of the new and untried schemes for CO 2 removal that are now being considered. 

Further, it is believed that basic improvements can be made in molecular-sieve 
C02~removal systems even at this stage in their development. Up to this time, 
molecular sieve 5A has been the choice for the CO 2 sorbent. However, at low 
partial pressures, (about 4 mm Hg and below) 4A has higher equilibrium capacity 
than 5A. It is possible that with smaller pore openings, the dynamic perfor- 
mance of 4A might be poorer than that of 5A. However, the comparison tests 
should be made. Also, at the present time, both Linde and Davison are market- 
ing 4A and 5A molecular sieves without binder. These materials have the 
potential of providing an immediate 20 percent increase in capacity as com- 
pared to their clay-binder counterparts. Both equilibrium and dynamic tests 
would be warranted. PredTrier sorbents that have low capacities for 0£ and N 2 
are desired. Molecular sieve 3A, or other large-cation substituted sorbents 
may be useful here. There are a number of design features that need exploring. 

One such feature, to limit the loss of atmospher i c gases (0£ and N 2 ); involves 
the pressure equalization of closed-off desorbing and adsorbing beds before 
switching. It is considered that a large fraction of atmospheric gases on 
the adsorbing bed would be quickly desorbed and saved by this operation. Other 
features such as multipoint desorption bear study. 

From the viewpoint of providing better basic data, more coadsorption 
equilibrium and dynamic testing would be fruitful. Computer program improve- 
ment could involve the effects of longitudinal or axial dispersion.- More 
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Information of rate effects with coadsorption would be welcomed. In some 
areas, rechecks of s i ngl e-sorbate equilibrium data would be idf interest. 
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PRECEDING page blank not f ilmed 


APPENDIX A 

LISTING OF REVISED MAIN PROGRAM S9950 
FOR EXECUTING BACK-TO-BACK HALF-CYCLE COMPUTATIONS 
USING THE SAME INITIAL CONDITIONS BUT DIFFERING VALUES OF GK AND DIF. 


* ELT S9950. 1.700121. 59453 , 1 


000001 

000002 

000003 

000004 

000005 

000006 

000007 

000008 

000009 

000010 

000011 

000012 

000013 

000014 

000015 

000016 

000017 

000018 

000019 

000020 

000021 

000022 

000023 

000024 

000025 

000026 

000027 

000028 

000029 

000030 

000031 

000032 

000033 


COMMON /MEAD/ HEADNGt 20 ) « l PR I NT. NLINE 

COMMON /BL0K2/ ABED( 41 ) . AUl). AVC(41),CPG(41),RH0G(41), 

1HXG( 41 ) > HXS( 41 ) . HXC( 4l ) .DlF(4l)»F(41) »C( 41 ) » VS( 2 ) >DVS(2) .DVS1( 2) i R 
2Sl<2»4J.>,RH0S<2>»UG<4l>»WM<2),UCC41>*NDXM.PSUl>.RH0SB<41>.DS<41>. 
3CSl(41},CS2(41),Cl(21.41),C2(21,41),Dl(2i,41),D2<41),PCl(41>,PC2<4 
41 ) »PC3t 41 ) < ASG( 41 ) * ASX t 41 ) , AGX( 41 ) » ClP( 41 ) » C2P( 41 ) > C3P< 41 > »D1P( 41 ) 

5, FR(2»41), RS< 2 ) .NDX1.NDR4.DX.DT. GK( 41 ) » DH( 41 ) . SK( 41 ) »P1( 41 ) > P2( 41 ) 

6. P3C 41 ) « WSt 41 ) ,CR1<2> 21 ) . CR2( 2 > 21) , CR3( 2. 21 ) . C3(21. 41 ) ,B(21> 41) . 
7Q(21. 41 ) .CPK 41 ) , CP2( 41) »X( 41 ) * VOIDFt 41 > .TIME 

8,AXC(41),RH0C<4l),CPC(41),T268.AVX(41).TKX<41>,CPX<41),RH0X(41) 

9 > NOG . PK ( 2 . 4 1 ) » PC02 1 . PH20 1 . CMR. GMW, TGI .PA .PT( 41) ,CPS( 41 > ,HSG( 41) 
COMMON /BU0K3/ W(21»4l), TG(4l). TS(41), TC(4l), TX(41). CYCLE 
DATA 111/77/ 

C 

REWIND 15 

READ (5.400D HEADNG 

4001 FORMAT (20A4) 

IF UII.EQ.77) WRITE (15) W. TG. TS. TC. TX 
REWIND 15 
C 

2 IF UIt.NE.77) READ (15) W. TG. TS. TC. TX 
REWIND 15 

III * 0 
C 

READ (5.4002) GKR. DIFR, I PRINT 

4002 FORMAT (2F10. 0.110) 

DO 3 1*1,41 

GK( I ) = GKR 

3 DIFU) a DIFR 
NLINE = 100 
CALL MADSOR 
GO TO 2 

END 
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PRECEDING PAGE BLANK NOT FILMED 

APPENDIX B 

LISTING OF REVISED PRINTED OUTPUT PROGRAM PRADSB 


9 ELT 89979,1,700121, 59451 , 1 


000001 

000002 

000003 

000004 

000005 

000006 

000007 

000008 

000009 

000010 

000011 

000012 

000013 

000014 

000015 

000016 

000017 

000018 

000019 

000020 

000021 

000022 

000023 

000024 

000025 

000026 

000027 

000028 

000029 

000030 

000031 

000032 

000033 

000034 

000035 

000036 

000037 

000038 

000039 

000040 

000041 

000042 

000043 

000044 

000045 

000046 

000047 

000048 

000049 

000050 

000051 

000052 

000053 

000054 

000055 

000056 

000057 

000058 


SUBROUTINE PRADSB 
C 
C 

DIMENSION AVL0(41) 

COMMON /BL0K2/ ABEDHl) , A( 41 > . -AVC< 41 > ,CPG< 41 ) ,RH0G( 41) , 
1HXG(41),HXS(41),HXC(41),D1F<41),F<41),C<41),VS(2),DVS<2>»DVS1<2),R 
2Sl<2,41),RH0S<2>,UG(4i>,WM(2),UC(41),NDXM,PS<41),RH0SB<41>,DS(4l), 
3CS1<41),CS2C41),C1(21,41),C2(21,41),D1(21 i41),D2<41),PC1<41),PC2(4 
41)»PC3(41),ASG(41)»ASX(41)» AGX( 41 > * C1P( 41 ) , C2P( 41 ) »C3P(41)»D1P(41) 

5. FR(2,41),RS(2),NDX1,NDR4,DX,DT,GK(41),DH<41),SK(41),P1(41),p2<41> 

6, p3(4l),WS(4i>,CRl<2.21),CR2(2,2l).CR3(2,21),C3<21,4l),B<21,4l), 
7Q(2l»41)»CPl(41) » CP2( 41 ) » X< 41 ) » V0IDF( 41 ) , TIME 

8 . AXC( 41") ,RH0CI41 ) *CPC( 41 ) »T268» AVX( 41 ) «TKX( 41 ) »CPX( 41) »RHOX( 41 ) 

9 , N0G,PK(2, 41 ) ,PC02I ,PH20I , GMR, GMW, TGI , PA , PT< 41 ) , CPSC 41 ) , HSG( 41 ) 
DOUBLE PRECISION W 

COMMON /BL0K3/ W( 21 , 4.1 ) » TG< 4l ) , T5< 41 ) , TC< 41) , TX( 41 ) , CYCLE 
COMMON /BLOK*/ POUT (10) »T1MET( 10) ,NBC0UT,NPSET(3) 

COMMON /BL0K6/ NCYCLT 

COMMON /BL0K8/ DTO» TSl< 41 > , TS2< 41 ) , TX1 ( 41) , TX2( 41 ) , TC1< 41) , TC2( 41 > 

COMMON/ BLOKIO/ NPR I NT , OTM AX , NOTCON 

COMMON /BL0K11/ T0TC02, T0TH20, SUMPTM, WTACMS, WTSG 

COMMON /BL0K12/ NTEMP 1 

COMMON /BL0K13/ -Wj , Tl 

COMMON /BL0K14/ NCYCLE 

COMMON /BL0K16/ NDXMAC, PC02C, VOLCAB, RC02C 
COMMON /BL0K177 NSTART 

DIMENSION AXS( 41 ) , AXG( 41), AGS ( 41 ) , HSX( 41 ) , HGX( 41) ,HGS( 41 ) , ACX( 41) , 
1HCX( 41 ) 

EQU I VALENCE ( ASX , A XS ) , < AGX , AXG ) , ( ASG , AGS ) , ( HXS , HSX ) , ( HXG , HGX ) ,'< HSG 
l.HGS) , ( AXC, ACX) , (HXC, HCX ) 

E8UIVALENCE(NDX,NDX1) 

COMMON /HEAD/ HEADNGC 20 ) , JPRINT, NLINE 
DIMENSION C0NCQ2U1), C0NH20C41) 

DATA NLINE/100/ 

C 

1 IF (NLINE.NE.100) GO TO 2 
NSPACE * 4 

IF ((NDX.EQ. 8 ) , OR, ( NDX«EQ. 18 ) ) NSPACE = 3 
IF ((NDX.EQ. 9) . OR. ( NDX.EQ. 19) ) NSPACE * 2 
IF ((NDX.EQ. 10). OR. (NDX.EQ. 20)) NSPACE * 1 
C 

2 AVPH20 * SUMPTM/TIME 

AVH20P =GMR«AVPH20»18 . /( PA*GMW ) 

C 

IF ((IPRINT.EQ.O). AND. (NLINE. LT.60)) GO TO 10 

IF ( ( IPRINT. GT. 0 ) . AND. ( (60”NLINE ) .GE« (NDX+9 ) ) ) GO TO 10 
WRITE (6,4000) HEADNC 
4000 FORMAT (1H114X20A4) 

NLINE 8 1 

WRITE (6,4003) GK< 1 > ,DIF( 1 ) 

4003 FORMAT ( 1H+100X’ GK( 1 )* ' F7 ,5, • DIF(1)*'E7.2) 

IF ( IPRINT. GT.O) GO TO 10 
WRITE (6,4004) 


4004 FORMAT (1H016X'PARTIAL PRESSURES CONCENTRATIONS — — TEMPE 


25H ADS4X'TIME , 9X'(MMHG)'10X , (LB/LB PR GAS) M5X' (F )' 34X* < LB/LB) ' / 

3 ' CYCLE (HRS) C02 H20 C02 H20 GAS SO 

4RBENT COOLANT HXCORE AVG 1 2 3 4 5') 
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APPENDIX B (Continued) 


000059 

000060 

000061 

000062 

000063 

000064 

000065 

000066 

000067 

000068 

000069 

000070 

000071 

000072 

000073 

000074 

000075 

000076 

000077 

000078 

000079 

000080 

000081 

000082 

000083 

000084 

000085 

000086 

000087 

000088 

000089 

000090 

000091 

000092 

000093 

000094 

000095 

000096 

000097 

000098 

000099 

000100 

000101 

000102 

000103 

000104 

000105 

000106 

000107 

000108 

000109 

000110 

000111 

000112 

000113 

000114 

000115 

000116 

000117 

000118 


NUINE * NUINE + 4 
C 

10 TIMEM~60.*TIME 
C 

IF ( I PR I NT «EQ*Q ) GO TO 15 

IF ( <NSPACE*EQ*1) , AND* (NUNE*EQ*1> ) GO TO 14 
DO 13 1*1 , NSPACE 
WRITE (6,4013) 

4013 FORMAT (1H ) 

13 NUNE * NUINE + 1 

14 WRITE (6,100) NCYCUE, TIME, TIMEM, OT 
100 FORMAT UW 1 ADSORPTION GYCUE' 13,* — - * 

1 5HT I ME * , F9 , 5 , 1X2HHR ,F12.3, IX, 3HMIN, 

2 5X , 15 HT I ME I NCREMENT#F7 , 5 , 1X2HHR ) 
NUNE s NUINE +1 

C 

15 N1 = NDXM + 1 
NDR3-NDR4-1 

DO 20 N s 1, NDXM 
SUMMS- 0 ♦ 0 

SUMMS- SUMMS+0 , 5* < W < 1 , N ) +.W < NDR4 » N ) > 

IF( NDR4 ,EQ « 2 ) AVLD(N) - SUMMS 
If ( NDR4 * EQ #2 ) GO TO 20 
Dl) 22 NR-2,NQR3 
22 SUMMS = SUMM$+W ( NR, N ) 

AVUD(N) -SUMMS/ NDR3 
20 CONTINUE 

DO 30 N a N1,NDX1 
SUMMS a 0 * 0 

SUMMS-SUMMS+O , 5* C W ( 1, N ) *W< NDR4, N ) > 

IF< NDR4 * EQ » 2 ) AVUD(N) * SUMMS 
IF(NDR4,EQ,2> GO TO 30 
DO 33 NR *2, NOR 3 
33 SUMMS -SUMM$+W< NR, N) 

AVUD(N> s SUMMS/NDR3 

30 CONTINUE 
SUM-0. 

IF( NDXMAC ,EQ* 0) GO TO 35 
DO 31 N- 1, NDXMAC 

31 SUM 5 SUM+AVUD(N)*ABE0(N>»DX*RM0^9(N> 

AVMSUD - SUM/WTACMS 

SUM s 0 . 

35 CONTINUE 

IF((NDX1-NDXM) .EQ, 0) GO TO 40 
DO 32 N s N1,NDX1 

32 SUM - SUM+AVUD(N)*ABED<N>*DX*RH0SB(N) 

AVSGUD=SUM/WTSG 

40 CONTINUE 
C 
C 

IF (I PR I NT .EG * O') GO TO 60 
C 

WRITE (6,4049) 

4049 FORMAT (1H0 


V PARTI AU PRESSURES CONCENTRATIONS — ——TEMPERATURES* 

2 — . — — ■ **»*«*«*****#*# interior SORBENT LOADINGS*********.***** 

3*'/lH 

4'AXIAL'6X'(MMHG>'1OX'(LB/L0 PR GAS >' 15XM F> * 41X' < LB/LB ) ' /1H 

5 'NODE C02 H20 C02 H20 GAS SORBENT COOLAN 

6T HXCORE A VC 1 2 3 4 5 6 7' 
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APPENDIX B (Continued) 


000119 

000120 

000121 

000122 

000123 

000124 

000125 

000126 

000127 

000120 

000129 

000130 

000131 

000132 

000133 

000134 

000135 

000136 

000137 

000138 

000139 

000140 

000141 

000142 

000143 

000144 

000145 

000146 

000147 

000148 

000149 

000150 

000151 

000152 

000153 


7) 

DO 50 N=1.NDX 

RMW = (PA - PKU»N> - PK<2iN) ) » GMW 
C0NC02(N) * PKd,N)*44./RMW 
C0NH20( N ) * PK<2.N)*18./RMW 

50 WRITE ( 6 , 4050 ) N, PKU*N), PK(2»N>* C0NC02<N>, CONH20(N>» TG(N), 

1 TS(N># TC(N>, TX(N) » AVLD(N>» <W<NR»N>.NR»1,NDR4> 

4050 FORMAT (1H 13* 3XG9 ,5, IX. G9 , 5, 1X2F7 . 5* IX, 4F8 , 1. 2XF7 , 4 , IX. 7F7 , 4 > 
WRITE (6,205) AVMSLD* AVSGLD 

205 FORMAT(30H0AVG C02 LOADING IN M.S. BED s F8,4.1X,5HLB/LB,7X, 

130H AVG H20 LOADING IN S.G,. BED * F8.4, 6H LB/LB) 

AVRC02*T0TC027TIME 
AVRH20 - TOTH20/T I ME 
WRITE (6, 204) AVRC02, AVRH20 

204 Format ( 27h tjme avg C 02 adsorp rate * 

1 F8.4, 6H LB/HR *10X, 

227H TIME AVG H2.0 ADSORP RATE » 

3 F8.4, 7H LB/HR ) 

WRITE(6, 206) AVPH20* AVH20P 

206 FORMATC21H TIME AVG EXIT PH20 - , F10.4, 3H MM ,18X, 

140HTIME AVG RATE OF M.S, POISONING BY H20 * * G12,4,10H LB H20/HR) 
NLINE ■« NLINE + 8 + NDX 
RETURN 
C 
C 

60 RMW = (PA - PK<1,1> » PK(2, 1 >') * GMW 
C0NC02( 1) • PK(1,1)«44./RMW 
C0NH20(1) * PK(2,l)*18,/RMW 

WRITE (6,4060) NCYCLE, TIME* PK(1,1>* PK<2,1), C0NC02(1>, 

1 C0NH20( 1 ) • TG ( 1 ) * TS( 1 ) * TC(1). TXd). 

2 AVLD(l) * ( W ( NR , 1 ) * NR - 1 » NDR 4 ) 

4060 FORMAT (1H 13 *F10 ,5,3x, .69 ,5, 1XG* . 5 , 1X2F7 ,5, 1X4F8 . 1 , 2XF7 ,4 , 1X5F7 . 4 > 
NLINE « NLINE + 1 
RETURN 
END 


201 




page 

APPENDIX C 


blank not 


filmed 


LISTING OF REVISED EQUILIBRIUM DATA PROGRAM PKEQ 


9 ELT PKEQ# 1 * 70X103# 55012 # 1 


000001 

000002 

000003 

000004 

000005 

000006 

000007 

000008 

000009 

000010 

oooon 

000012 

000013 

000014 

000015 

000016 

000017 

000018 

000019 

000020 

000021 

000022 

000023 

000024 

000025 

000026 

000027 

000028 

000029 

000030 

000031 

000032 

000033 

000034 

000035 

000036 

000037 

000038 

000039 

000040 

000041 

000042 

000043 

000044 

000045 

000046 

000047 

000048 

000049 

000050 

000051 

000052 

000053 

000054 

000055 

\0Q0Q56 

000057 

000058 


c 

c 

G 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

e 

c 

c 

c 

c 

c 

c 

c 

c 

0 

c 

c 

c 

c 

0 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FUNCTION PKEQ (K# W# T> 

THIS VERSION OF 'PKEQ* PROVIDES TABLE-LOOKUP AND INTERPOLATION OF 
SORBENT/SORBATE EQUILIBRIUM data. CURVE FITTING (OR MAP FITTING) 
IS NOT USED. 

EQUILIBRIUM DATA ARE STORED VIA 'DATA* STATEMENTS IN ARRAYS AS 
FOLLOWS... GIVEN FOR SORBENT/SORBATE COMBINATION NO.l# SAME FOR 2. 


•TIDATA(J)* 
*NT1 * 


fWlDATAU)* 

* NW1 * 

'P1DATAC InJ)* 


» TEMPERATURE VALUES (F) FOR WHICH P-VS-W CURVES 
(ISOTHERMS) ARE STORED, 
a NUMBER OF ISOTHERMS STORED. NOTE# WHEN THE 
COMPUTATION TO BE PERFORMED IS FOR ISOTHERMAL 
CONDITIONS. * NT! * MUST BE 6 OR LESS. OTHERWISE# 
IT MAY BE 7 OR LESS. 

« LOADING VALUES (LB/LB), COMMON FOR EACH 
ISOTHERM# MAY BE 20 OR LESS, 

« NUMBER OF LOADING VALUES, 

« PARTIAL-PRESSURE VALUES (MM/HG) CORRESPOND ING 
TO THE LOADING AND TEMPERATURE VALUES STORED IN 
' WlDATA* AND 'T1DATA*. 


actual values of all data are to be used in data statements, 

IF LOGARITHMIC STORAGE OF PRESSURE OR LOADING VALUES IS DESIRED# 
THE INDICATORS WHICH FOLLOW WILL AUTOMATICALLY ALLOW SAME.,, 

* ILOGWl * * Oi NORMAL STORAGE OF LOADING DATA# 

a 1# LOGARITHMIC STORAGE, 

1 ILOGPl * a 0# NORMAL STORAGE OF PARTIAL-PRESSURE DATA# 

« 1# LOGARITHMIC STORAGE, 


NOTE,, FOR EACH ISOTHERM# A PARTIAL-PRESSURE VALUE FOR EACH VALUE 
OF 'WlDATA* MUST BE ENTERED. BY THE NATURE OF THE METHOD BY WHICH 
EQUILIBRIUM DATA ARE TAKEN# PARTIAL-PRESSURE DATA WILL PROBABLY 
NOT BE AVAILABLE FOR ALL VALUES OF LOADING IN 'WlDATA', THAT IS# 
DATA WILL NOT BE AVAILABLE FOR ONE OR BOTH ENDS OF THE ISOTHERM 
WITH RESPECT TO THE DATA ENTERED IN * WlDAf A* • TO OBTAIN# 

MISSING VALUES ONE MAY EXTRAPOLATE GRAPHICALLY THE EXISTING 
EQUILIBRIUM MAP. ALTERNATIVELY# ONE MAY ENTER ZEROES (0.0) IN THE 
DATA STATEMENTS# IN WHICH CASE# THIS PROGRAM WILL FILL IN THE 
PRESSURE ARRAYS BY LINEAR EXTRAPOLATION, NOTE# ZEROES BETWEEN REAL 
VALUES ARE NOT ALLOWED, ZEROES ARE NOT ALLOWED AT ALL IN 'WlDATA* 


IF DATA FOR ONLY ONE SORBENT/SORBATE COMBINATION IS TO BE STORED# 
THE * NT- * AND 'NW-* VALUES FOR THE OTHER COMBINATION SHOULD BE 
STORED AS ZERO# IN THE EVENT THAT THIS PROGRAM IS ENTERED TO 
INTERPOLATE DATA WHICH HAS NOT BEEN STORED# PROGRAM EXIT WILL 
OCCUR. 


IN THIS PROGRAM INTERPOLATION OF ISOTHERM DATA IS 8Y MEANS OF 
SUBROUTINE *LAGIN2' , 

INTERPOLATION (OR EXTRAPOLATION) BETWEEN (OR BEYOND) ISOTHERM 
DATA# AT CONSTANT LOADING W# IS BY MEANS OF THE RELATION 
LN( P ) a A/(T*459,6) + B, 


WHEN ONLY ONE DATA CURVE IS STORED FOR A SORBENT/SORBATE 
COMBINATION# IT WILL BE USED DIRECTLY IN ALL CIRCUMSTANCES, 

WHEN AN ISOTHERMAL RUN IS DESIRED# AS DENOTED WHEN 'NTEMP* »0# 

AT THE FIRST CALL TO *PKEQ' # AND PROVIDED THERE ARE TWO OR MORE 
ISOTHERMS STORED# THE PROGRAM WILL GENERATE AN ISOTHERM FOR THE 
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000059 

000060 

000061 

000062 

000063 

000064 

000065 

000066 

000067 

000068 

000069 

000070 

000071 

000072 

000073 

000074 

000075 

000076 

000077 

000078 

000079 

000080 

000081 

000082 

000083 

000084 

000085 

000086 

000087 

000088 

000089 

000090 

000091 

000092 

000093 

000094 

000095 

000096 

000097 

000098 

000099 

000100 

000101 

000102 

000103 

000104 

000105 

000106 

000107 

000108 

000109 

000110 

000111 

000112 

000113 

000114 

000115 

000116 

000117 

000118 


C 

C 

C 

C 


c 


c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


GIVEN TEMPERATURE AND STORE it IN THE 'NT-*1'«AREA OF 'P-OATA’. 
THIS CURVE WILL THEN BE USED DIRECTLY FOR ALL SUBSEQUENT CALLS, 

COMMON /BL0K12/ NTEMP 

DIMENSION T1DATA< 7 ) # W1DATA<2Q), P1DATA( 20* 7 ) * 

1 T2QATA(7>, W2DATA<20># P2DATAC 20* 7 ) * 

2 TDATA< 7 , 2 > * WDATA(20» 2) • PDATA ( 20, 7, 2 > * NTD<2)i NWD<2)* 

3 IL0GW( 2 ) t IL0GPC2) 

EQUIVALENCE <TlDATA<i>*TOATA<l*l)>* ( W1DATA ( 1 ) * WD AT A ( 1 * 1 )> * 

1 ( T2DATA< 1 ) * TPATA( 1*2) > * < W20ATA< 1 > * WDATA( 1# 2 ) > , 

2 CP10ATAU*!) *PDATA< 1*1*1 )>* 

3 <P2DATA( 1* 1 ) * POATAM * 1*2 )>* 

4 (NTl.NTDU))* (NW1*NWD(1))» 

5 (NT2*NTD<2>>* <NW2*NWO< 2 > > * 

6 ( ILOGW1* JLOGWC 1 ) ) » (ILOGP1* ILOGPU))* 

7 <IL0GW2*lL0GW<2) ) # ( ILOGP2* ILOGPC 2 > > 

DATA INlTL/1/ 

DATA DELTA/0.5/ 

****##&##*#**#****«************* ******* »***«*» **#*#*• ***♦#* •••»*** 


DATA FOR DAVIS0N-5A/C02 AT 70F. 

DATA NT1/1/# NW1/12/ 

DATA T1DATA< 1 >/70 •/ 

DATA ( W1DATAI I )# I *1*12) 

1/ 0*0065, Q. 0088* 0*0130* 0»0l77* 0*02l5* 0.Q324* 0*0489, 
2 0.0615, 0 . 0671 * 0*0714* 0,0757/ 

DATA ( (PlDATAt I * J ) * I *1* 12) , J*l* 1) 

1/ 0,1* 0.2* 0.4, 0*7, 1.0* 2.0* 4.0* 

2 6,0* 7,0* 8.0* 9,0/ 

DATA ILOGW1/1/, ILOGPl/l/ 


D5A/C-70 
05A/C-70 
D5A/C*70 
D5A/C-70 
0.0553* D5A/C-7Q 
D5A/C-70 
D5A/C-70 
5.0, D5A/C-70 
D5A/C-70 
05A/C-70 


DATA NT2/0/* NW2/0/* ILOGW2/0/* ILOGP2/Q/ 


* 


IF ( ( NWD ( K > « QT . 0 ) , AND . ( NTD ( K ) • GT . 0 ) ) GO TO 1 
WRITE (6,4000) K 

4000 FORMAT MHO »*****ERROR — DATA FOR SORBENT/SORBATE COMBINATION’ 
112** HAS NOT BEEN STORED IN PKEQ.») 

CALL EXIT 
C 

c 

1 IF (INITL.EQ.O) GO TO 100 
C 


2 INI TL = 0 
C 
C 

C CONVERT ENTERED DATA TO LOGARITHMS (IF REQUIRED) AND FILL IN 
C PARTIAL-PRESSURE ARRAYS BY EXTRAPOLATION. 

C 

10 IF ( ( NW1 * EQ , 0 ) • OR , ( NTl • EQ * 0 ) ) GO TO 20 
C 

IF (ILOGWl.EQ.O) GO TO 12 
DO 11 I 2 1*NW1 
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000119 


ii 

WIDATA(I) * ALOGCWlDATAtl)) 

000120 

C 



000121 


12 

DO 19 U*l. NT1 

000122 

C 



000123 



I ZERO A * 0 

000124 



1ZEROB * 0 

000125 



II * 0 

000126 

C 



000127 



DO 15 1=1, NW1 

000128 



IF (PIOATAU.J)-GT.O.O) GO TO 14 

000129 



IF UI.GT.O) GO TO 13 

000130 



IZEROA = IZEROA + 1 

000131 



GO TO 15 

000132 


13 

IZEROB = I ZEROS + 1 

000133 



GO TO 15 

000134 

C 



000135 


14 

II * 1 

000136 



IF (ILOGPl,CTiO) FlOATAMsJ) * AtOGC PlDATAl I # J> > 

000137 


15 

CONTINUE 

000138 

C 



000139 



IF (IZEROA. EQ. 01 GO TO 17 

000140 

c 



000141 



SLOPE * ( P1D AT AC I ZERO A+2 . J ) - PlOATAt IZEROA+l.J) > t 

000142 



l ( W1D AT A ( I ZEROA+2 ) - WIDATAt IZEROA+l) ) 

000143 



DO 16 1=1, IZEROA 

000144 



PlDATA(I.J) * P1DATA( I ZEROA+1 » J ) + 

000145 



l $L0PE*<W1DATA(I) - W10ATA<iZER0A*l)> 

000146 



P = EXP(P1DATA(I,J>) 

000147 



WRITE <6,4001) WIDATA(I). P 

000148 

4001 

FORMAT ( 1H0 F5,4, E13.5) 

000149 



IF < ( ILOGPl vEQ.O) , AND, (PlOATAt I » JKLT «1«QE'7 ) > P1DATAC 

000150 


16 

CONTINUE 

000151 

C 



000152 


17 

IF ( I ZERO'S. EQ,0> GO TO 19 

000153 



11 = NW1 * 1 - IZEROB 

000154 



SLOPE = (PlOATAt 11-2, J) - PlOATA(Il-l.J)) t 

000155 



l ( WIDATAt 11—2) - WlOATA(Il-l) ) 

000156 



DO 18 1=11, NW1 

000157 



PlOATAt I, J) « PlOATAt 11*1, J) + 

000158 



1 SL0PE«< W10ATA t l) - WIDATAt 11-1 > > 

000159 



P = EXP (PlOATAt I , J) J 

000160 



WRITE (6,4001) WIDATAt 1 ) , P 

000161 


18 

CONTINUE 

000162 

C 


* 

000163 


19 

CONTINUE 

000164 

e 



000165 


20 

IF ( (NW2.EQ.0) .OR , (NT2.EQ.0) > GO TO 50 

000166 

c 



000167 



IF ( IL0GW2.EQ. 0 ) GO TO 2? 

000168 



DO 21 1=1, NW2 

000169 


21 

W2D AT A ( I ) * AL0G(W2DATA( I ) ) 

0Q017Q 

c 



000171 


22 

IZEROA = 0 

000172 



IZEROB = 0 

000173 



II =0 

000174 

c 



000175 



DO 29 J=1»NT2 

000176 

c 



000177 



DO 25 1=1, NW2 

000178 



IF t P2D AT A ( I , J ) . GT , 0 , 0 ) GO TO 24 
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000179 

000180 

000181 

000182 

000183 

000184 

000185 

000186 

000187 

000188 

000189 

0QQ19Q 

000191 

000192 

000193 

000194 

000195 

000196 

000197 

000198 

000199 

000200 

000201 

000202 

000203 

000204 

000205 

000206 

000207 

000208 

000209 

000210 

000211 

000212 

000213 

000214 

000215 

000216 

000217 

000218 

000219 

000220 

000221 

000222 

000223 

000224 

000225 

000226 

000227 

000228 

000229' 

000230 

000231 

000232 

000233 

000234 

000235 

000236 

000237 

000238 


c 

e 

s 


c 


G 

C 

c 

£ 

c 


c 


0 

e 

c 

P 

§ 

c 


IF (II.GT.O) GO TO 23 
IZEROA * IZEROA + 1 
GO TO 25 

23 IZEROB = IZEROB + 1 
GO TO 25 

II =1 

IF (ILOGP2.GT.O) P2DATAU.J) * AL0G(P2DATA< I . ) 

25 CONTINUE 

IF ( IZEROA. EQ.O) GO TO 27 

SLOPE = (P2QATAI IZER0A+2.J) - P20ATAUZER0A+1. J) > / 

1 <W20ATA( 1ZER0A+2) r W2DATA! IZER0A*1 > ) 

DO 26 1=1 . IZEROA 

P2DATAU.J) * P2DATA< IZER0A+1.J) * 

1 SL0PE*<W2DAT*(i) - W2QATA! IZERQA+1) ) 

IF < C I L0GP2 . EG • 0 ) .AND . ( P2DATA ( I » J 1 • Lt . 1 • 0E*7 > ) P2D*TAU. J)»I*l.E-7 

26 CONTINUE 

27 IF ( IZEROB. EQ.O) GO TO 29 

II = NW2 +1 - IZEROB 

SCOPE » < P2DATA 1 1 1-2 » J ) - P2D*TAH1-1*U>> / 

1 <W20ATA( 11-2) - W20ATACU-1) > 

DO 28 1=11. NW2 

28 P2DAT A ( 1 1 J > * P2DATAI U-l.J) + 

1 SLOPE*! W2DATAU) - W2DATA ( 1 1-1 ) ) 

29 CONTINUE 

IF (NTEMP.NE.O) GO TO 100 

IF ( ( NT1 • LE > 1 ) .AND . < NT2 . LE > 1 > ) GO TO 100 

GENERATE ISOTHERMS FOR TEMPERATURE *T» — TO BE USEO IN ISOTHERMAL 
PERFORMANCE PREDICATIONS. 

50 K1 =1 

K2 = 2 

IF (NTl.GT.l) GO TO 5l 
K1 * 2 

51 iF (NT2.GT.1) GO TO 52 

K2 =1 

52 DO 53 KK=K1»K2 

NWK = NWO(KK> 

NTK = NTD(KK) 

DO 53 1*1. NWK 

53 CALL EQMAP (KK, WDATA(I.KK). T. PDATAU , NTK+1 . KK ) . 0) 
»*»***#«**«..****#*****#***»*»******»»*»•••***»•***•***•**•*•*•**** 
ACTUAL INTERPOLATION. 

ioo we = w 

IF (W.LE.0.0) WC * 1.0E-5 
CALL EQMAP (K. WC. T. PKEQ. 1) 

RETURN 


SUBROUTINE EQMAPIL. WE. T. P. INDEX) 
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000239 

000240 

000241 

000242 

000243 

000244 

000245 

000246 

000247 

000248 

000249 

000250 

000251 

000252 

000253 

000254 

000255 

000256 

000257 

000258 

000259 

000260 

000261 

000262 

000263 

000264 

000265 

000266 

000267 

000268 

000269 

000270 

000271 

000272 

000273 

000274 

000275 

000276 

000277 

000278 

000279 

000280 

000281 
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200 NT 3 NTD(L> 

NW = NWD( L ) 

C 

IF (NT.GT.l) GO TO 210 
JT 3 1 

GO TO 290 
C 

210 IF ( ( !NDEX#EQ» 0) • OR# (NTEMP* NE# G > > GO TO 220 
JT 9 NT * 1 

GO TO 290 
C 

220 IF ■<T.GE.TOATA<t#L.n GO TO 230 
J1 =1 

J2 3 2 

GO TO 260 
C 

230 IF <T.Lt,TDATA<NTil>> GO TO 240 
J1 3 NT « 1 

J2 3 NT 

GO TO 260 
C 

240 DO 250 J=1»NT 

IF (T.GT.TDATA('J+1,L)> GO TO 250 
J2 - J + l 

J1 3 J 

GO TO 260 
C 

250 CONTINUE 
C 

260 IF CABS(T-TDATACJ1»U >.GT. DELTA) GO TO 270 
JT 3 J1 

GO TO 290 
C 

270 IF (ABS(T*TDATA(J2#U).GT. DELTA) GO TO 280 
JT 3 J2 

GO TO 290 
C 

280 WLOOK 3 WE 

IF ( ILOGWC L> i GT.O) WLOOK 3 ALOG< WE.) 

CALL LAGIN2 (331. WDATA( 1 #L ) # NW# 2# WLOOK# Pli PDATA< 1# J1#L > > 

CALL LAGIN2 (332# WDATA(1#L># NW, 2, WLOOK# P2, PDATA( 1# J2#L) ) 

IF ( I LOOP ( L ) * GT # 0) GO TO 282 
PI 3 ALOG(Pl) 

P2 3 AL0G(P2) 

C 

282 A 3 (P2-P1) / U./(TOATA< J2#.U + 459.6)-l./(T0ATA( J1#L)M59.6> ) 

B 3 PI - A/ ( TD AT A ( Jl# L ) +459 » 6 ) 

P 3 EXPX A/< T+459 , 6 ) + B) 

RETURN 

C 

290 WLOOK 3 WE 

IF ( ILOGW(L ) ,GT • 0 ) WLOOK 3 ALOG(WE> 

CALL LAGIN2 (333# WDATA( 1#L) # NW# 2, WLOOK# P, PDATA( 1# JT#L ) ) 

IF (ILOGP(L).GT.O) P 3 EXP(P) 

RETURN 

C 

C 

SUBROUTINE LAG JN2 (Mi X# NP# ND# XA, YAi ¥) 

C 

C INTERPOLATION SUBROUTINE# 


207 



000299 

000300 

000301 

000302 

000303 

000304 

000305 

000306 

000307 

000308 

000309 

000310 

000311 

000312 

000313 

000314 

000315 

000316 

000317 

000318 

000319 

000320 

000321 

000322 

000323 

000324 

000325 

000326 

000327 

000328 

000329 

000330 

000331 

000332 

000333 

000334 

000335 

000336 

000337 

000338 

000339 

00034Q 


APPENDIX C (Continued) 


C ADAPTED FROM THE GENERAL •LAGIN2* PROGRAM, 

C 'X* ARRAY INCREASING. LINEAR TWO-POINT INTERPOLATION, NO MESSAGE. 

C 

DIMENSION X(2> » Y<2> 

C 

300 ILO = 1 

INI » NP 

C 

C DETERMINE IF EXTRAPOLATION IS NECESSARY. 

310 IF (XA - X<1>) 320, 370. 311 

311 IP (XA - X<NP)J 330, 380, 321 


C 


C 

C 


C 

C 

C 

c 


c 


320 11 * 1 

12 * 2 

GO TO 360 

321 II « NP - 1 

12 * NP 

GO TO 360 

BINARY LOOKUP. 

330 LOOK * (ILO ♦ IHI ♦ 1) / 2 

IF (XA - X ( LOOK ) ) 331, 390, 332 

331 IHI * LOOK 

GO TO 340 

332 ILO * LOOK 

340 IF ( ( IHI-ILO) . GT . 1 ) GO TO 330 

350 11 - ILO 

12 * IHI 


LINEAR, TWO-POINT INTERPOLATION, 

360 YA » Y ( II) + < Y< 12 )*Y< 1 1 ) ) * ( XA-X< 1 1 > ) / < X< 1 2 ) -X< 1 1 ) > 
RETURN 


370 YA 8 Y( ILO) 

RETURN 

380 YA * Y< IHI ) 

RETURN 

390 YA » Y(LOOK) 

return 


END 
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LEAST SQUARES PROGRAM TO CORRELATE 
EQUILIBRIUM DATA 


INTRODUCTION 

This program generates a polynominal equation to correlate equilibrium 
data using a stepwise multiple regression least square analysis. The 
polynominal equation is in the form: 

P = C + B ( X, + B 2 X 2 + + Bn Xn 

where: 

P = equilibrium pressure, mm Hg 

B, C = constants 

X = function of T and W 

T = temperature, °F 

W = absorbent loading, lb sorbate/lb sorbent 

The program uses a least squares subroutine (CORLAT) which is an adaptation 
of Share Fortran program G2 ER MPR2 No. 477. The correlation developed by 
this program is used in the molecular sieve adsorption/desorption program 
to predict regenerative carbon dioxide bed performance. The polynominal 
developed must have derivatives which are continuous and nonzero; there- 
fore polynominals up to the third order were used to correlate the 
equi librium data. 

PROGRAM LOGIC 

Main Program 

The flow chart for the main program is shown in Figure D-l. The main 
program reads in all the control data (N,K) and the equi librium data points 
(T, W, P) . N is the total number of data sets (T, W, P) and K is the 
number of different equations to be correlated by subroutine CORLAT. The 
present program contains 35 equations; these equations are shown in Table D-l. 

The variables for each of the thi rty-f ive equations are defined within 
the main program, the constant and the coefficients for a given equation are 
determined by subroutine CORLAT. After the constant and coefficients have 
been determined, control is returned to the main program. For a given 
correlation the deviation at each data point and the average percent error 
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MAIN 

PROGRAM 


( ) 


READ 

N,K 


READ I = I , N 
T(l) , w(i), 
vP(l) 


DiS 800J=I ,K 


DEVS=0.0 

NiS>DATA=0 



1 

f 

D«! 7501=1 , N 


1 

N&DATA 

4 

= NtlDATA 
1 





Figure D-l. Flow Chart for Main Program 
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TABLE D-l 

CORRELATING EQUATIONS USED IN THE MAIN PROGRAM 


CsrohSTANT 
L'lrt.'ATURAl. LOG 
P=l ,f >t:S>SlJRF, MM-HG 
TrTFtlRERATlIRcf PEG. F 
W=LC'/\P» WT, PERCE- it >120 OR CO?. 

dQUAl I-0N5 A T rERRTEP . 

i . LfUP) r-i/ ( i 4t6<;)4t ! i ( V.) / ( T 4 460 ) + LN I a) 4 LM I V?) **? >0 

a . LM( P) --i / C T+460 ) 1 1/ i ( v; ) /|T +«60 ) +Lf. ( „■ » + l II I W ) *+.3+.l/LN.( W) **3-1 / ( T+460 ) **2+C 
LG ■(•-) 5 - * / < r ♦ •:.» r !j ) * I . i i ( 1 . ) / ( T + 4*,0 ) + 1.1 1 (,„•)**? + 1 N ( W ) * * 3 + C 

‘I . >- > ( i ' ) ~~i/ < ... ♦ 4 o'.* ) 1 1 f - < v; ) +U : ( * ) **2+C 

■i , L. ,.(.*) =.!..« i •) /v:i. V ( T 4 460 ) +L * ! ( >' ) +Lf i ( W ) ##?-4.C 

•v, ■ . ( . ) 1 1 •' ) / t 74 4. OJ +'LM ) +.LUI V) **P+C. 

?. a.')-! i< ■ ) + 1/ (T+4 U '-).)-+l/<T++60) **.?+LU<W)*>2.+C 

i . L V >--1/ 1 1 + 4',0) (■•..) ♦+?4LN(«)+i/(T+4fO)4-C 

V. <-•(!■. =t.v* ! ( l/< 7+4. 0.) 4*r+U-. (W)/1 3. 6+LI I (L0(5( too*”') ) + 1/ (T+460 ) **2 ILMI W ) /13.6“ 

L.itl 0. ( 1 00'. ) ) 4 C > 

l " . rl-Jt I t.h I • ) / ; 3 , <> +L J I LUG 1 1 0 0* V* ) ) 4 1/ ( T+460 ) ** ? I Li) I ’*> ) /t 3.6-1 N ( 1.00 ( 1 00 *.W 

) ) ) 4 C ) 

1 1 . .. f”)-!. ( 1/ ( 7 + 4i ft ) ** 2 ( LT ■ { V. ) / 1 3. G-Ll-i I LOG 1 1 00*>/ ) ) )+C) 

U. . . j < 4- ) - l <(■)/( T44LO) >0. :1 1 >')+*?/ 1 T+4u0 ) .) + I LN I W ) *#3./ 1 T+460) V+C 

I . ... .it ) -:l.i i ( • ' ) / ( T 44, ()') + ( |_; i ( V. ) **2/ ( T +4uQ ) **2 ) + (LN( V- ) **3/1 T+460 ) **3 ) +C 

it. i. 4 O' )=(.•!( . I 4 ( L (. ( ) 4*?/( r 44611) ) + (LHI ty ) * *3/ ( T+ 400 ) **2 ) +C 

r.-j. .. .0 )*4.-.+l. ( ;) /< r+4G >) **?4i ri(W)**2/(T44«so)4C 

J. .. ! iit 1 - L I ( V ) / ( 7 44 •vO+LM ' ) * + 3/ ( T + 460 ) +Lli («/)*♦? 4C 

1 / , L- i 14") il.t-O ) 4 (. !(•■, ) 4,2/ ( r 4 460) +LH Ite ) + + 3/ 1 T + 4&0 ) +F 

l . 1.. T +'*!}!) ) 4 i ! ; ( ,y ) 4*2/ ( T +460) * 4td+LM ( W) **3/ 1 T 4460 ) **3+C 

t 4 ?. • .(>' ) -L t( +4.. i ( .v ) ti, 241. 4* 3 / { T + 4 uO ) + C 

cU. !-L 1(!')4L>'(< )+*2/| T>46u.) 4LII(V)**3 + C 

Ri. L i ) -L. t( V) / ( T + Rj ‘i ) +!.!•> I 4)*+2+LMlW) 4*3 + C 

t i. . f r- ) .” I / | T 4 4»'.f > +i. ‘ ) * * 5*Li l| *)*♦?/ ( T+460) **3 +C 

4. 3 . .. ( :') -u.G t ) 4 *3+1-1 ( -i ) / | T 4460 ) +LM7 it ) *+?/ IT +460 ) **?+C 

*4. ; V i ■ ) — t_i »{••?) 4 i_. i t ;.)*+, Ml. 1(1 ) .+ *3+ 1 / ( T +460 ) +C 

» >. i (i ) -t.M ( ) * ♦ 24 i- ! ( .;) 4* .) + !,"!< .) /(T + 4o0 )4 1 /(T+460 ) 4C 

.(! >=:.,<(•■)* + 5+L (.,-)* *2/ (T+460) +t r- I »)/( t+ 460 )+l/ (T+460 )+C 
ft. U. (i )=L.;i(/ J + i 4/( : 4», >0)4i.' ( .\)+*2/(T+4»'P)+tM(X)/(T+46n)+l/(T+4f>0)+r 

c X. (I )=L.:(‘/)4 4 ‘*Ui.( . ) 4*.! + L'i( «. ) 4 1/ (T+460 ) **? + C 

*V. i. vU > -1..4 1 •“ )4 *3-+Li .( /!)** 3*1 f.u,,) /(T+4o0 ) +*2 + l/(T+460) **?+C 

6... t a! ) -L. i ('■)* 42/1 M + bO) **?+L! (¥)/( T + 460 ) **2+1/ (T+460 ) **?+I.G ( vJ) **.3+C 

jl . i_ (i ) -L. . ( ■. )* t. V ( T+460) **£+Li S(R)*+?/( T+460 >**2+LI'{W) /( T+460 ) **2+ 

j / i I +4 >0 j + 4 ,:+C 

32. 4.1(1 )=LU(i-;)/( T44, ; r. )-)/( r+46.0) **2+l/LN C * ) ♦*3+LG 1 4') +1/1 t+460 ) ++3+C 

35. Li'!(H)=UI(W)/ (T+460 ) +1 /L>i ('-.') **3-1/ ( T+46C ) **2 +LM ( W ) **3+LN ( W V + 1/ ( T+46U ) **3+C 

5*4 . I...(i'i =L t() )/( I +«r.'i) +LIJ( . ) +LM V ) *♦?+)-.* (V ) **3+1 /U 4 (W) **3-1/ ( T+460) **2+ 

\/ ( T 4'. >1! ) 4* 3 +C • 3 • 

56. 4 (. in-1/ ( r +.»6i ) .Li (i ) /(T + +6T) 4l r( ., ) +I.NIV) +*3 + l/IJJ( ./) +*3-1./( T+460 ) **2+ 

’... !( )* *;'4i~ 
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for all data points are determined. Assuming a normal error distribution, 
the average percent error for a correlation is determined from: 



( 100 ) 


where: 

TER = average percent error 
N = number of sets of data points 


DEV(N) = deviation in pressure, at data set N 
p(N) = pressure at data set N 


Subroutine CORLAT 


Subroutine CORLAT receives the values of the defined variables from 

the main program and performs the stepwise multiple regression analysis 

in the form of P = B, X, + B« X„ * ... B X + C, where n is the number 

ii c c n n 

of independent variables, 59 or less. The number of sets of data points 
that may be analyzed at one time Is 300. This restriction can be increased 
by enlarging the size of the dimension statement. Some properties of 
the stepwise regression follow: 

1. A number of intermediate regression equations are obtained 
adding one variable at a time. The variable added is that one which 
makes the greatest improvement in "goodness of fit". This provides 
valuable statistical information at each step in the calculation. 

2. A weighting factor can be assigned to each set of data. Certain 
sets of data may have higher reliability than others, and weighting factors 
can be applied in relation to their reliabilities. Weighting factors 
were not used for the equilibrium data. 

3. Only significant variables are included in the final regression. 

4. A variable which is approximately a linear combination of other 
independent variables is not entered into the regression. This eliminates 
the possibility of degeneracy and permits single precision floating 
point operations. 
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Calculation Procedure 

The setpwise multiple regression procedure consists of four parts, 
as follows: 

I. Calculation of weighted sums of squares and crpss products. The 
following values are calculated: 


a, 

b. 


c. 


Weighted number of data = £w t 
Weighted sums of variables = Ew x. 

t I • t 

Weighted sums of squares and cross products = Ew^x 


i.t x j.t 


where 


t = I to m 
i = I to n + I 
j = i to n + I 


Note: 


x 


n+l 


= y 


2. Calculation of weighted means, residual sums of squares, and 
residual sums of cross products. 

a. Mean = x. = Ex x. ./Ew 

I t I . t t 

b. Residual sums of squares and cross products = 


S-t Zw t x l ,t x j.t 


Ew^x. ^Ew.x. . 

t i ,t t i . t 


Ew. 


where 

t = I to m 
i = I to n+l 
j = i to n+l 

y = Vi 
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3. Calculation of simple correlation coefficients, r . j ? 

a. Definition . • • 

. .. * 

;X . x . = Residual sums of squares 

J ,r : * ' ' • ,i •> ‘ ? •• 

X X 

i j = Residual sums of cross products 

» I > \ 

b. Calculations 

a, = x,x. ‘ 1 

r . . 

J i j at partial correlation coefficient 

o.cr . 

• J ! 

r . . = r. . 

J • i J 

r . . == 1 . 000 

where 

i = I to n+ 1 
j = i+ 1 to n+l 

4. Stepwise calculation of regression coefficients and reliability. 

In the stepwise procedure, intermediate results are used to give 
valuable statistical information at each step in the calculation. These 
intermediate answers are also used to control the method of calculation. 

A number of intermediate regression equations are obtained by adding 
one variable at a time, thus giving the following intermediate equations: 

a. P = B j X j + C 

b. P = 6 1 X j + B 2 Xg + C, etc. 
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The coefficients for each of these Intermediate equations and the 
reliability of each coefficient are obtained by the stepwise procedure. 
The values and reliability may vary with each subsequent equation. The 
coefficients represent the best values when the equation is fitted by 
the specific variables included in the equation. The variable is added 
that makes the greatest improvement in "goodness of fit" or, stated 
another way, gives the greatest reduction in variance of the dependent 
variable. 

A variable may be indicated to be significant at an early stage and 
enter the regression equation. After several other variables are added 
to the regression equation, a variable in the equation may be indicated 
to be insignificant. Under this situation the stepwise regression 
procedure will remove the insignificant variable before adding, an 
additional variable. Thus, at the various steps in the regression 
procedure, only those variables which are significant will be included 
in the regression equation. 

USAGE 

Main Program 


All control words must be assigned in the main program and stored 
in the COMMON storage area. The sequence of storage must be the same as 
that of CORLAT. The common statements are: 

COMMON NTPIN, NTPOUT, TOL, EFIN, EFOUT, NOPROB, INVAR, NODATA 
COMMON IFLIST, IFWT, IFSTEP, IFRAW, IFAVE, IFRESD, IFCOEN, IFPRED 
COMMON IFCNST 
COMMON CP 

Tapes used by CORLAT are defined in the main program. The input tape 
is denoted by NTPIN, and is written by the main program in binary. The 
output tape is denoted by NTPOUT; it is BCD. 

For later use in the main program, the constant and the coefficients 
obtained as a result of correlation are stored in the common area. CNST 
is the constant, and the coefficients 8|, Bg, etc., are stored in the 
array COEN. If the coefficients are to be used in the main program, 

CNST and COEN must appear in the COMMON statement and, in addition, 

COEN (60) must appear in the DIMENSION statement of the main program. 

The coefficients are then obtained using COEN (i), I =1, 2 -■ — 60. 
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TOL 

EFIN 

EFOUT 

NOPROB 

INVAR 

NODATA 

The 

IFLIST 

IFWT 

IFSTEP 

IFRAW 

IFAVE 

IFRESD 

IFPRED 

IFCNST 


Tolerance (normally .001) 

F level to enter variables in the regression 

F level to remove variables in the regression 

Problem number (up to 8 integer digits) 

The total number of variables, including the dependent 
variable, (if N is the number of independent variables 
then INVAR = N+l) 

Number of sets of data. (This need not be determined before 
object time, but "counted" in the main program). 


following control words should never be negative. 


Controls listing of the input data to the subroutine. 

If IFLIST = 0, data are printed on the output tape. 

If IFLIST = I, data are not printed 

Controls weighting of sets of data. 

If IFWT = I, then all weights are 1.0 and are not to be 
included. 

If IFWT = 0, then weights are assumed to be given with the 
input data. 

Controls printing of results of intermediate steps; printed 
if = 0, not printed if == I. Usually IFSTEP = 0, as this 
gives very valuable information. 

Controls printing of raw sums, sums of squares, and sums of 
cross products; printed if = 0, not printed if = I. 

Controls printing of average values of the variables; printed 
if = 0, not printed if = I. 

Controls printing of residual sums and squares; printed 
if = 0, not printed if = I. 

Controls printing of predicted values and the comparison 
with actual values; printed if = 0, not printed if = I. 

Controls whether a constant appears in the regression 
equation, if = I, then the constant is set to zero; if = 0, 
then the constant term is included and determined by the 
analysis. 
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Input Data " 

Data which are used by CORLAT is written by the main program 
(equations 1-35) in binary form and the number of this tape assigned to 
NTPIN; e.g., NTPIN = 15, if the data is written on tape 15. One set 
of data is stored in the following form: (all floating point): the run 
number, RUNNO; the independent variables (x.> Xg, Xj — ) ; and then the 
dependent variable (XP) . " 

The first input card contains the number of data sets (N) and the 
number of equations (K) to be correlated. The following cards contain 
one set of equilibrium conditions (T, W, P) . If variable weights are < 
used, the weight for each set of data follows the independent variable. 
Sets of data are written one after the other until all data are used. 
Then the tape is rewound by the main program and control shifts to 
CORLAT by the call statement: CALL CORLAT. 


Should an error occur in the execution of CORLAT, a statement is 
written on the output tape and control is returned to the main program 
by means of CORLAT statement number 910 RETURN. By suitable coding of 
the main program, several different regression problems can be handled 
sequentially. This involves proper rewinding, writing, and rewinding 
of the input tape, and assignment of control words in the main program. 

It should be noted that control words are not changed at any time by 
CORLAT. Thus, control words, once assigned for one problem, do not need 
to be reassigned for a subsequent problem, unless they are to be changed. 
The normal return to the main program is also by subroutine statement 
910 RETURN. 

Output 


For each equation calculated in the main program the following is 
printed from: 


Subroutine CORLAT 

1. The number of the equation calculated. 

2. The number of variables used in the equation, including the 
dependent variable. 

3. The number of input data points and degrees of freedom. 

4. The statistical level to enter and remove each independent 
variable. 
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For each variable entered into the correlation computation the 
fol lowing is given. 

1. The computation step number 

2. The variable, by number, entered into the correlation computation. 

3. F level of the variable entered. This is the significance of 
the variable to the equation. The larger the F level, the more important 
the entered variable. 

4. The standard error of the polynomial with the entrance of the 
given variable. 

5. The constant C, for the polynomial. 

6. The coefficients of the polynomial, 

7. The diagonal elements of the final vector matrix used to develop 
the final polynomial. 

Main Program 

1 . Data point sequence number, 

2. Data point (T,W, P) 

3. Calculated value pf polynomial, per data point. 

4. Deviation between the data point and its calculated value. 

5. Percent error of the calculated point. 

6. Average percent error of the calculated polynomial. 

Operation Details 

This program was written in Fortran IV for the Univac 1108 system. 

The core storage for the total program is 134741 octal words. Compilation 
and execution time for all 35 equations is one minute. 

SAMPLE PROBLEM 

The input equilibrium data for water on molecular sieve type 5A is 
shown in Table D- 2, the output obtained with the two best data fits 
(equations 33 and 34) are shown in Table D-3. Table D-4 contains the 
source listing for the program. 
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TABLE D-2 

INPUT DATA FOR SAMPLE PROBLEM 
t>Al n I iPvl (H£0-5A) 
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